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I Abstract: We consider a heavy quantum particle with an internal degree of freedom moving on the d-dimensional lattice Z** (e.g., a heavy 
atom with finitely many internal states). The particle is coupled to a thermal medium (bath) consisting of free relativistic bosons (photons 

1 I or Goldstone modes) through an interaction of strength A linear in creation and annihilation operators. The mass of the quantum particle is 

assumed to be of order A^^, and we assume that the internal degree of freedom is coupled "effectively" to the thermal medium. We prove that 
O ^. the motion of the quantum particle is diffusive in d > 4 and for A small enough. 
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^ 1 Introduction 

> 

00 1.1 Diffusion 
l> . 

T— I Diffusion and Brownian motion are central phenomena in the theory of transport processes and nonequilibrium 
statistical physics in general. One can think of the diffusion of a tracer particle in interacting particle systems, the 
\^ diffusion of energy in coupled oscillator chains, and many other examples. 

From a heuristic point of view, diffusion is rather well-understood in most of these examples. It can often be 
. successfully described by some Markovian approximation, e.g. the Boltzmann equation or Fokker-Planck equation, 
depending on the example under study. In fact, this has been the strategy of Einstein in his ground breaking work 
of 1905, in which he modeled diffusion as a random walk. 

However, up to this date, there is no rigorous derivation of diffusion from classical Hamiltonian mechanics or 
unitary quantum mechanics, except for some special chaotic systems; see Section fLS-ll Such a derivation ought to 
d allow us, for example, to prove that the motion of a tracer particle that interacts with its environment is diffusive 
at large times. In other words, one would like to prove a central limit theorem for the position of such a particle. 

In recent years, some promising steps towards this goal have been taken. We provide a brief review of previous 
results in Section [131 In the present paper, we rigorously exhibit diffusion for a quantum particle weakly coupled 
to a thermal reservoir. However, our method is restricted to spatial dimension d > 4. 



o 



^Postdoctoral Fellow FWO-Flanders at K.U.Leuven, Belgium, email: wojciech.deroeck@fys. kuleuven .be 



1 



1,2 Informal description of the model and main results 

We consider a quantum particle hopping on the lattice U^, and interacting with a reservoir of bosons (photons or 
phonons) at temperature (3^^ > 0. In the present section, we describe the system in a way that is appropriate at 
zero temperature, but is formal when /3 < oo. The total Hilbert space, J^/f, of the coupled system is a tensor product 
of the particle space, J^, with a reservoir space, Jf^. Thus 

-.^ J^i <E) JfR. (1.1) 

The particle space J^s is given by S^, where the Hilbert space describes the internal degrees of freedom 

of the particle, e.g., a (pseudo-)spin or dipole moment, and the particle Hamiltonian is given by the sum of the 
kinetic energy and the energy of the internal degrees of freedom 

Hs ■■= HsMn ® 1 + 1 «) i?S,spin (1-2) 

The kinetic energy is chosen to be small in comparison with the interaction energy, and this is made manifest in its 
definition by a factor A^, where A is the coupling strength between the particle and the reservoir (to be introduced 
below). Hence we set 

HsMn = A2e(P), (1.3) 

where the function e is the dispersion law of the particle and P is the lattice-momentum operator. The most natural 
choice is to take e{P) to be (minus) the discrete lattice Laplacian, — A. The energy of states of the internal degree of 
freedom is to a large extent arbitrary 

^^s,spin ^, for some Hermitian matrix F, (1.4) 

the main requirement being that Y not be equal to a multiple of the identity. 

The reservoir is described by a free boson field; creation and annihilation operators creating /annihilating bo- 
sons with momentum q eM.'^ are written as a* , Qq, respectively. They satisfy the canonical commutation relations 

[a*,a*]=0, [aq,a;,]^S{q-q% (1.5) 

where a* stands for either a or a*. The energy of a reservoir mode q is given by the dispersion law oj{q) > 0. To 
describe the coupling of the particle to the reservoir, we introduce a Hermitian matrix W ony and we write X for 
the position operator on P{Z'^). 

The total Hamiltonian of the system is taken to be 

Hx:=Hs + J dga;(g)a*a, + Ay"dg(e'«-^«)W^(g)0(g)a, +e-'«-^«)iy(g)0(^a;) (1.6) 

acting on <g) J^. The function is a form factor and A e K is the coupling strength. We write Hs instead of 
Hs fE) 1, etc. 

We introduce three important assumptions: 

1) The kinetic energy is small w.r.t. the coupling term in the Hamiltonian, as has already been indicated by the 
inclusion of A^ in the definition of i?s,kin. Physically, this means that the particle is heavy. 

2) We require a linear dispersion law for the reservoir modes, u}{q) = \q\,in order to have good decay estimates at 
low speed. This means that the reservoir consists of photons, phonons or Goldstone modes of a Bose-Einstein 
condensate. 

3) We assume that the amplitude of the wave front of a reservoir excitation (located on the light cone) has 
integrable (in time) decay. This is satisfied if the dimension of spac^l is at least 4. 

^Since the integrability in time is only needed for reservoir excitations, we can in principle also treat models in which the particle is 3- 
dimensional, but the reservoir is effectively 4-dimensional. 
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Additional assumptions will concern the smoothness of the form factor cj) and the "effectiveness" of the coupling 
to the heat bath (e.g., the interaction between the internal degrees of freedom and the reservoir, described by the 
matrix W, should not vanish.) 

The initial state, p^, of the reservoir is chosen to be an equilibrium state at temperature {3^^ > 0. For math- 
ematical details on the construction of infinite reservoirs, see e.g. fl0l l4ll2l. The initial state of the whole system, 
consisting of the particle and the reservoir, is a product state pg p^, with pa, a density matrix for the particle that 
will be specified later. The time-evolved density matrix of the particle ('subsystem') is called ps.t and is obtained 
by "tracing out the reservoir degrees of freedom" after the time-evolution has acted on the initial state during a 
time t, i.e., formally. 



(1.7) 



where Tr is the partial trace over Mr. We warn the reader that the above formula does not make sense math- 
ematically for an infinitely extended reservoir, since the reservoir state is not a density matrix on Mr. This 
is a consequence of the fact that the reservoir is described from the start in the thermodynamic limit and, hence, 
the reservoir modes form a continuum. Nevertheless, the LHS of formula jl.Ti can be given a meaning in the 
thermodynamic limit. 

The density matrix ps,t obviously depends on the coupling strength A, but we do not indicate this explicitly. We 
also drop the subscript S and we simply write pt, instead of ps.t/ in what follows. 
We will often represent pt as a !3§{,5^)-va{\xed function on 1/^ x U^: 

PtixL.Xj,) e ^(^), Xl,Xj, € Z". (1.8) 

Although this is not necessary for many of our results, we require the initial state of the particle to be exponentially 
localized near the origin of the lattice, i.e., 

\\pt{xL,XR)\\sgi.9') < Ce-*'l^^le-'^'l^«l, for some constants C, 5' > (1.9) 
Our first result concerns the diffusion of the position of the particle. 



1.2.1 Diffusion 

We define the probability density 

pt{x) ■.^Ttj/' pt{x,x) (1.10) 

where Tr^' denotes the partial trace over the internal degrees of freedom. The number pt (x) is the probability to 
find the particle at site x after time t. 

By diffusion, we mean that, for large t, 

,„W. (_L)*\,e.I,)-"^e.ri-l i^^D-'j,)}. (Ill) 

where the diffusion tensor D = Dx is a strictly positive matrix with real entries; actually, if the particle dispersion 
law e is invariant under lattice rotations, then the tensor D is isotropic and hence a scalar. The magnitude of D 
is inferred from the following reasoning: The particle undergoes collisions with the reservoir modes. Let tm be 
the mean time between two collisions, and let u„, be the mean speed of the particle (the direction of the particle 
velocity is assumed to be random). Then the mean free path is Vm x <,„ and the central limit theorem suggests that 
the particle diffuses with diffusion constant 

D ~ ^"""^ (1.12) 

tm 

The mean time tm is of order tm ^ A^^ since the interaction with the reservoir contributes only in second order. 
The mean velocity is of order Vm ^ because of the factor A^ in the definition of the kinetic energy. Hence 
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We now move towards quantifying ((LTT|. Let us fix a time t. Since pt (x) is a probability measure, one can think 
of Xt as a random variable with 

Prob(a:t = x) :— pt{x). (1.13) 
The claim that the random variable converges in distribution, ast /• oo, to a Gaussian random variable with 
mean and variance D is called a Central Limit Theorem (CLT). It is equivalent to pointwise convergence of the 
characteristic function, i.e., 

y e'^'^ Vtla;) e-^«-°9, forall^eM'', (1.14) 

and it is this statement which is our main result, Theorem l3.1l A stronger version of the convergence in ||1.14| | (also 
included in Theorem |3.1t implies that the rescaled moments of /if converge. For example, for i,i = 1, . . . ,d, 

iTr[p,X,] = i^a;,/.,(a;) ^ (1.15) 

X 

jTr[ptXiXj] = j^XiXjHt{x) ^ Aj, (1-16) 

X 

In fact, the first line (vanishing of average drift) is expected only if one assumes that the model has space inversion 
symmetry, which is assumed throughout. 

1.2.2 Equipartition 

Our second result concerns the asymptotic expectation value of the kinetic energy of the particle and the internal 
degrees of freedom. The equipartition theorem suggests that the energy of all degrees of freedom of the particle, 
the translational and internal degrees of freedom, thermalizes at the temperature f3~^ of the heat bath. We will 
establish this property up to a correction that is small in the coupling strength A. This is acceptable, since the 
interaction effectively modifies the Gibbs state of the particle. We prove that, for all bounded functions F, 

MF{HsMn)) M dfcF(A2e(fc))e-^^'^(^)+o(|An (1.17) 

PtiFiHs,sp^n)) ^ Fie)e-^^ + oi\\\% as A \ (1.18) 

where Z, Z' are normalization constants and the sum J^eespY ranges over all eigenvalues of the Hamiltonian Y. 
We note that the factor e^^^^^*^'') can be replaced by 1 (as in Theorem l3.2|l since we anyhow allow a correction term 
that is small in A and the function e{k) is bounded. For this reason, one could say that, for very small values of A, 
the translational degrees of freedom thermalize at infinite temperature {/3 — 0). 



1.2.3 Decoherence 

By decoherence we mean that off-diagonal elements pt{x, y) of the density matrix pt in the position representation 
fall off rapidly in the distance between x and y. Of course, this property can only hold at large enough times when 
the effect of the reservoir on the particle has destroyed all initial long-distance coherence, i.e., after a time of order 
A^^. Thus, there is a decoherence length 1/^dch and a decay rate g such that 

l|pt(a;L,x^)IU(^)<Ce-^-"l"--"-l+C'e-^'s*, as t ^ oo (1.19) 

for some constants C, C . The magnitude of the inverse decoherence rate ^dch is determined as follows: The time 
the reservoir needs to destroy coherence is of the order of the mean free time tm, while the time that is needed for 
coherence to be built up over a distance 1/^dch is given by {jdch x Vm)^^, where Vm is the mean velocity of the 
particle. Equating these two times yields 

(1.20) 

and hence, recalling that t„i ^ A^^ and Vm ^ A^, as argued in Section ri.2.1[ we find that 'jdch does not scale with A. 
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1.3 Related results and discussion 



1.3.1 Classical mechanics 

Diffusion has been established for the two-dimensional finite horizon billiard in ||6J. In that setup, a point particle 
travels in a periodic, planar array of fixed hard-core scatterers. The finite-horizon condition refers to the fact that the 
particle cannot move further than a fixed distance without hitting an obstacle. 

Knauf [27J replaced the hard-core scatterers by a planar lattice of attractive Coulombic potentials, i.e., the po- 
tential is V{x) = — J2jez^ \x-j\ ■ ^ case, the motion of the particle can be mapped to the free motion on a 
manifold with strictly negative curvature, and one can again prove diffusion. 

Recently, a different approach was taken in [24j: the authors of [24.| consider ad ^ 3 lattice of confined particles 
that interact locally with chaotic maps such that the energy of the particles is preserved but their momenta are 
randomized. Neighboring particles can exchange energy via collisions and one proves diffusive behaviour of the 
energy profile. 

1.3.2 Quantum mechanics for extended systems 

The earliest result for extended quantum systems that we are aware of, II30I , treats a quantum particle interacting 
with a time-dependent random potential that has no memory (the time-correlation function is S{t)). Recently, 
this was generalized in |25| to the case of time-dependent random potentials where the time-dependence is given 
by a Markov process with a gap (hence, the free time-correlation function of the environment is exponentially 
decaying). In [32J, we treated a quantum particle interacting with independent heat reservoirs at each lattice site. 
This model also has an exponentially decaying free reservoir time-correlation function and as such, it is very similar 
to ||25l . Notice also that, in spirit, the model with independent heat baths is comparable to the model of [24J, but, 
in practice, it is easier since quantum mechanics is linear! 

The most serious shortcoming of these results is the fact that the assumption of exponential decay of the correl- 
ation function in time is imrealistic. In the model of the present paper, the space-time correlation fimction, called 
t) in what follows, is the correlation function of freely-evolving excitations in the reservoir, created by interac- 
tion with the particle. Since momentum is conserved locally, these excitations cannot decay exponentially in time 
t, imiformly in x. For example, if the dispersion law of the reservoir modes is linear, then ^(x, t) is a solution of the 
linear wave equation. In d = 3, it behaves qualitatively as 

^{x, t) ~ ■|iTi5(c|i| — with c the propagation speed of the reservoir modes (1.21) 

\x\ 

In higher dimensions, one has better dispersive estimates, namely sup^ \'>}j{x, t)\ < (under certain conditions), 
and this is the reason why, for the time being, our approach is restricted to d > 4. In the Anderson model, the ana- 
logue of the correlation function does not decay at all, since the potentials are fixed in time. Indeed, the Anderson 
model is different from our particle-reservoir model: diffusion is only expected to occur for small values of the 
coupling strength, whereas the particle gets trapped (Anderson localization) at large coupling. 

Finally, we mention a recent and exciting development: in [13J, the existence of a delocalized phase in three 
dimensions is proven for a supersymmetric model which is interpreted as a toy version of the Anderson model. 

1.3.3 Quantum mechanics for confined systems 

The theory of confined quantum systems, i.e., multi-level atoms, in contact with quasi-free thermal reservoirs has 
been intensively studied in the last decade, e.g. by |3l|23l[l2l. In this setup, one proves approach to equilibrium for 
the multi-level atom. Although at first sight, this problem is different from ours (there is no analogue of diffusion), 
the techniques are quite similar and we were mainly inspired by these results. However, an important difference 
is that, due to its confinement, the multi-level atom experiences a free reservoir correlation function with better 
decay properties than that of our model. For example, in {23\, the free reservoir correlation function is actually 
exponentially decaying. 
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1.3.4 Scaling limits 



Up to now, most of the rigorous results on diffusion starting from deterministic dynamics are formulated in a 
scaling limit. This means that one does not fix one dynamical system and study its behaviour in the long-time limit, 
but, rather, one compares a family of dynamical systems at different times. The precise definition of the scaling 
limit differs from model to model, but, in general, one scales time, space and the coupling strength (and possibly 
also the initial state) such that the Markovian approximation to the dynamics becomes exact. In our model the 
natural scaling limit is the so-called weak coupling limit: one introduces the macroscopic time t :— and one 
takes the limit X\0,t oo while keeping r fixed. In that limit, the dynamics of the particle becomes Markovian 
in r (as if the heat bath had no memory) and it is described by a Lindblad evolution. The long-time behavior of 
this Lindblad evolution is diffusive. This is explained in detail in SectionlH One may say that, in this scaling limit, 
the heuristic reasoning employed in the previous sections to deduce the A-dependence of the diffusion constant 
and the decoherence length becomes exact. The same scaling is known very well in the theory of confined open 
quantum systems as it gives rise to the Pauli master equation. This was first made precise in 1 9 1 . 

If we had set up the model with a kinetic energy of 0(1) (instead of 0{X^)), then one should also rescale space 
by introducing the macroscopic space-coordinate x '-^ ^^x- The reason for this additional rescaling is that, between 
two collisions, a particle with mass of order 1 moves during a time of order A^'^, and hence it travels a distance of 
order A^^. The resulting scaling limit 



is often called the kinetic limit. In the kinetic limit the dynamics of the particle is described by a linear Boltzmann 
equation (LBE) in the variables (x, r). The convergence of the particle dynamics to the LBE has been proven in 1 14j 
for a quantum particle coupled to a heat bath, and in |17J for a quantum particle coupled to a random potential 
(Anderson model). The long-time, large-distance limit of the Boltzmann equation is the heat equation, which 
suggests that one should be able to derive the heat equation directly in the limiting regime corresponding to 



This was accomplished in fl^.TSl for the Anderson model. An analogous result was obtained in |28J for a classical 
particle moving in a random force field. 

1.3.5 Limitations to our result 

Two striking features of our model are the large mass, of order A~^, and the internal degrees of freedom described 
by the Hamiltonian i?s,spin = Y ■ Physically speaking, these choices are of course not necessary for diffusion, they 
just make our task of proving it easier. Let us explain why this is so. First of all, once the mass is chosen to be of 
order A^^, the internal degrees of freedom are necessary to make the model diffusive in second order perturbation 
theory. Without the internal degrees of freedom, it would be ballistic. This is explained in Section fOl in particular, 
it can be deduced immediately from conservation of momentum and energy for the processes in Figure |3l Note 
also that the dependence on A is chosen such that the kinetic term i/s,kin — X^e{P) is comparable to the particle- 
reservoir interaction in second order of perturbation theory (both are of order A^). The large mass ensures that 
the position of the particle remains well-defined for a time of order A^^, which permits us to sum up Feynman 
diagrams in real space. 

Further, we note that our result requires an analyticity assumption on the form factor see Assumption 12.31 
This assumption ensures that the free reservoir correlation function ^p{x, t) is exponentially decaying for small x, 
even though it has slow decay on the lightcone, as explained in Section ri.3.21 

1.4 Outline of the paper 

The model is introduced in Section|2]and the results are stated in Section|3l In SectionlH we describe the Markovian 
approximation to our model. This approximation provides most of the intuition and it is a key ingredient of the 
proofs. Section |5] describes the main ideas of the proof, which is contained in the remaining Sections [SE] and the 
four appendices IAIIDI 



x-^\-'^x, t^X-'^t, A\0 



(1.22) 



X A-^^+^'x, t A"(2+2«)^^ A \ 0, for some k > 0. 



(L23) 
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2 The model 

After fixing conventions in Section lZTl we introduce the model. Section l2l2l describes the particle, while Section lZSl 
deals with the reservoir. In Section l234l we couple the particle to the reservoir, and we define the reduced particle 
dynamics Zt. Section |231 introduces the fiber decomposition. 

2.1 Conventions and notation 

Given a Hilbert space S , we use the standard notation 

i^p(^) js-e ^(^),Tr (S'*S')P/2 < 00} , 1 < p < 00, (2.1) 
with SSr,^{S) = ^{S") the bounded operators on S', and 

\\S\\p:^ (Ti-[{S*Sr^']y^\ \\S\\ := \\S\U (2.2) 

For bounded operators acting on ^p((f ), i.e. elements of I^{^p{(S')), we use in general the calligraphic font: 
V, W, T, .... An operator X e ^{S') determines an operator ad(X) € .^{^p{S')) by 

'Ad{X)S := [X, S]=XS- SX, S G 3§p{<ff). (2.3) 

The norm of operators in ^{^p^S")) is defined by 

We will mainly work with Hilbert-Schmidt operators (p = 2) and, unless mentioned otherwise, the notation ||>V|| 
will refer to this case. 

For vectors v G C'', we let Rev,lmv denote the vectors (Re , . . . , Re Vd) and (Im wi , . . . , Im Vd), respectively. 
The scalar product on is written as v ■ v' and the norm as \v\ y/v ■ v. 

The scalar product on a general Hilbert space S" is written as (•,•), or, occasionally, as {■,■)<?■ All scalar products 
are defined to be linear in the second argument and anti-linear in the first one. We use the physicists' notation 

\(p){(p'\ for the rank-1 operator in 3§{S') acting as (ys" i-> ((/?', (2.5) 

We write Tg^S') for the symmetric (bosonic) Fock space over the Hilbert space S" and we refer to flO] for defin- 
itions and discussion. If w is a self -adjoint operator on S', then its (self-adjoint) second quantization, drs(w), is 
defined by 

n 

drs(a;)Sym(<^i (g) . . . (g) (pn) ^ Syni(v3i (g) . . . (g) ujipi g) . . . g) (2.6) 

i=l 

where Sym projects on the symmetric subspace of and ipi, . . . , ipn G S'. 

We use C, C to denote constants whose precise value can change from equation to equation. 



7 



2.2 The particle 



We choose a finite-dimensional Hilbert space , which can be thought of as the state space of some internal 
degrees of freedom of the particle, such as spin or a dipole moment. The total Hilbert space of the particle is given 
by := P{'Zi'^, S^) ~ ^ y (the subscript S refers to 'system', as is customary in system-reservoir models). 

We define the position operators, Xj, on ^ by 



{Xjip){x) = Xjifix), 



(2.7) 



In what follows, we will almost always drop the component index j and write X = {Xj ) to denote the vector- 
valued position operator. We will often consider the space ,3% in its dual representation, i.e. as L'^(T'', S^), where 
T'' is the d-dimensional torus (momentum space), which is identified with L^([— tt, tt]'', .y). We formally define the 
'momentum' operator P as multiplication by A; € T"^, i.e.. 



PLp{k) = fc(^(fc), k e [-TT, nf, if e L^{T'^, y) 



(2.8) 



Although P is well-defined as a boimded operator, it does not correspond to a continuous function on T'', and it 
is not true that [Xj, Pj] = — i. Throughout the paper, we will only use operators F{P) where F is a function on T'' 
that is extended periodically to M'^. We choose such a periodic fimction, e, of P to determine the dispersion law of 
the particle. The kinetic energy of our particle is given by A^e(P), where A is a small parameter, i.e., the 'mass' of 
the particle is of order A^'^ 

The energy of the internal degrees of freedom is given by a self-adjoint operator Y G B{,y), acting on as 
{Yip){k) = Y{ip{k)). The Hamiltonian of the particle is 

Hs X^s{P) (g)l + l(g)Y (2.9) 

As in Sectionlll we will mostly write s(P) instead of s{P) (E) 1 and Y instead of 1 (E)Y. 

Our first assumption ensures that the Hamiltonian Hg ~ Y + \^e{P) has good regularity properties 

Assumption 2.1 (Analyticity of the particle d5mamics). The function e, defined originally on T^, extends to an analytic 
function in a neighborhood of the complex multistrip of width > 0. That is, when viewed as a periodic function on R'^, e is 



analytic (and bounded) in a neighborhood Oj 



i[—5e,5^])'^. Moreover, e is symmetric with respect to space inversion, i.e., 
e(k)=e(-k). (2.10) 



Furthermore, we assume there is no v & R'^ such that the function k ^ v ■ Ve{k) vanishes identically and that e does not 
have a smaller periodicity than that o/T'', i.e., we assume that 



e{k) = e{z + k)forall kcT"^ 



(2.11) 



The most natural choice for e is e{k) — 2(1 ~ cos(/ci)), which corresponds to —e{P) being the lattice Lapla- 

cian. As already indicated in Section ri.2.1[ the symmetry assumption l|2.101 is necessary to exclude an asymptotic 
drift of the particle. 

By a simple Paley- Wiener argument. Assumption 12.11 implies that one has exponential propagation estimates 
for the evolution generated by the operator e{P). Indeed, from the relation 



■x 



-ite{P+u) 



< Qle(\lii^v\)\t\ 



with (7e(7) := sup|j,j-jp|<^ |Ime(p)|, one obtains 

(e-"^(^))(x,.,XK) 



< Q~lWL~XR\^qc{'l)\t\ 



for llm^l < 5f 



for 7 < (5e 



(2.12) 



(2.13) 



where we write S{xliXb) for a ^(^)-valued 'matrix element' of 5 G ^(J^). 
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2.3 The reservoirs 

2.3.1 The reservoir space 

We introduce a one-particle reservoir space f) = L'^{M.'^) and a positive one-particle Hamiltonian uj > 0. The 
coordinate q E M.'^ should be thought of as a momentum coordinate, and u acts by multiplication with a function 

{u:^){q)^uj{q)p{q) (2.14) 

In other words, u is the dispersion law of the reservoir particles. The full reservoir Hilbert space, M^i, is the 
symmetric Fock space (see Section IZT] or IfTOl ) over the one-particle space \), 

^R:=r,(()) (2.15) 

The reservoir Hamiltonian, iJ^, acting on J^^ir is then the second quanitzation of w 

i?R, := drs(w) = j dquj{q)a*aq. (2.16) 

with the creation/ annihilation operators a*, Uq to be introduced below. 

2.3.2 The system-reservoir coupling 

The coupling between system and reservoir is assumed to be translation invariant. We choose a 'form factor' 
(j) e L^(R'') and a self-adjoint operator W = W* e ."^{.y) with || W|| < 1, and we define the interaction Hamiltonian 
HsR by 

i/sR:=y"dg(e'«-^(g)H/(g)0(g)a, +e^''?-^«)W^(g)0(^a;) on » ^R, (2.17) 
where a,, a* are the creation/ annihilation operators on f) satisfying the canonical commutation relations (CCR) 

[aq,a*q,] = S{q-q'), [a*,a*] = (2.18) 
with a* standing for either a or a* . We also introduce the smeared creation/ annihilation operators 



a 



{ip):= dqip{q)a* a{^) := dg^(g)a„ ip e L\R^). (2.19) 



In what follows we will specify our assumptions on iJsR/ but we already mention that we need [W, Y] ^ Q for 
the internal degrees of freedom to be coupled effectively to the field. 

2.3.3 Thermal states 

Next, we put some tools in place to describe the positive temperature state of the reservoir. We introduce the 
density operator 

Tfi = {c^'^ - l)-^ on f} = L^{R'^). (2.20) 

Let € be the *algebra consisting of polynomials in the creation and annihilation operators a{(p), a* {(p') with ip, (p' e 
f). We define as a quasi-free state defined on £. It is fully specified by the following properties: 

1) Gauge-invariance 

pi[a*M^p^^[aip)]=0 (2.21) 

2) The choice of the two-particle correlation function 

pi[a*{p)aip')] p^^[a*{^)a*{p')] \ {p^'lTp^) \ 

pi[a{<p)a{ip')] pi[a(p)a*{ip')] J { + W) j ^'^ 
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3) The state pj^ is quasifree. This means that the higher correlation functions are related to the two-particle 
correlation function via Wick's theorem 

pi[a*{^,)...a*{^2n)]^Y. n /'rK('^^)«*(^^)] (2.23) 

where stands for either a* or a, and Vn is the set of pairings tt, partitions of {1, ... , 2n} into n pairs (r, s). 
By convention, we fix the order within the pairs such that r < s. 

The reason that it suffices to specify the state on £ has been explained in many places, see e.g. l4l l20l[T0l 



2.3.4 Assumptions on the reservoir 

Next, we state our main assumption restricting the type of reservoir and the dimensionality of space. 
Assumption 2.2 (Relativistic reservoir and d > 4). We assume that the 

dimension of space rf > 4 (2.24) 
Further, we assume the dispersion law of the reservoir particles to be linear; 

coiq) \q\ (2.25) 
For simplicity, we will assume that the form factor cf) is rotationally symmetric and we write 

0(g) = (/.(|9|), qeW' (2.26) 
We define the "effective squared form factor" as 

, . f ^-^|0(|w|)P uj>0 

where we are abusing the notation by letting cj denote a variable in R. Previously, u! was the energy operator on 
the one-particle Hilbert space and as such, it could assume only positive values. Indeed, at positive temperature, 
the function ipi^) plays a similar role as |0(|cij|)p at zero-temperature: It describes the intensity of the coupling to 
the reservoir modes of frequency uj. Modes with a; < appear only at positive temperature and they correspond 

physically to "holes". One checks that = C^"^, which is Einstein's emission-absorption law (i.e. detailed bal- 

ance). This particle-hole point of view can be incorporated into the formalism by the Araki-Woods representation, 
see e.g. HHHOl. 

The next assumption restricts the "effective squared form factor" ip. 

Assumption 2.3 (Analytic form factor). Let the form factor be rotation-symmetric (j){q) = (^(|g|), as in (|2.26l l, and let ip be 
defined as in i2.27i . We assume that ?/'(0) — and that the function uj — ij{uj) has an analytic extension to a neighborhood 
of the strip R + i[(5R,, d^],for some (5r > 0, such that 

sup / duj\-4j{uj)\ < oo. (2.28) 

R+iX 

We note that Assumption |2]3] is satisfied (in d > 4) if one chooses: 

m) ■■= -^m) (2.29) 

Vm 

with a fimction on R with cj) = ^{uj) and analytic in the strip of width S^, and such that | |2.28|I holds with 
substituted for 
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The motivation for Assumptions 12.21 and 12.31 will become clear in Section 15.11 where we discuss the reservoir 
space-time correlation function ip{x,t). 

The last assumption is a Fermi Golden Rule condition that ensures that the spin degrees of freedom are effect- 
ively coupled to the reservoir. To state it, we need the following operators 

Wa:= J2 ie'iY)WleiY), a G sp(ad(r)) (2.30) 

e, e' £ spy 
e — e' = a 

Note that the variable a labels the Bohr-frequencies of the internal degrees of freedom of the particle. 
Assumption 2.4 (Fermi Golden Rule). Recall the function ijj as defined in H2.27b . The set of matrices 

■= ^^^{a)Wa,a e sp(ad(r))} C (2.31) 

generates the complete algebra .^{S^). This means that any S € ^(^) which commutes with all operators in is 
necessarily a multiple of the identity. We also require the following non-degeneracy condition 

• Every eigenvalue ofY is nondegenerate (multiplicity 1) 

• For all eigenvalues e,e',e", e'" ofY such that e ^ e! , we have 

e'-e = e"'-e" ^ e' = e"'flnrfe" = e (2.32) 

This condition implies that all eigenvalues o/ad(F) are nondegenerate, except for the eigenvalue 0, whose multiplicity is given 
by dim ,y. 

The strict nondegeneracy condition on Y, in contrast to the condition on is not crucial to our technique 
of proof, but it allows us to be more concrete in some stages of the calculation. In particular, the matrices Wa^o, 
introduced above in | |2.30|I , can be rewritten as 

Wa = {e\We)x\e'){e\, (2.33) 

where e, e' are the unique eigenvalues s.t. e — e' = a ^ 0, and we denoted the corresponding eigenvectors by 
the same symbols e, e'. The condition that .^w generates the complete algebra, can than be rephrased as follows: 
Consider an imdirected graph with vertex set spy and let the vertices e and e' be connected by an edge if and only 
if 

7/'(e'-e)|(e,iye')|V0 (2.34) 



(note that this condition is indeed symmetric in e, e', as long as f3 < oo). Then Assumption 12.41 is satisfied if and 
only if this graph is connected. 

Assumptions of the type above have their origin in a criterion for ergodicity of quantum master equations 
due to [34, 19 J. In our analysis, too. Assumption 12.41 is used to ensure that the Markovian semigroup At (to be 
introduced in SectionlU has good ergodic properties. This can be seen in Section lC.l.ll in Appendix ICl 

2.4 The dynamics of the coupled system 

Consider the Hilbert space :— ^ (g) . The Hamiltonian Hx (with coupling constant A) on ^ is (formally) 
given by 

Hx := Hs + Hn + XHgn (2.35) 

If the following condition is satisfied 

< oo, (2.36) 

then HsR is a relatively boimded perturbation of Hs + and hence Hx is a self-adjoint operator. One easily 
checks that l|2.36t is implied by Assumptions 12.2] and 12.31 
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For the purposes of our analysis, it is important to understand the dynamics of the coupled system at positive 
temperature. To this end, we introduce the reduced dynamics of the quantum particle. 

By a slight abuse of notation, we use to denote the conditional expectation SSiM'^^), given by 

pi{S®R) = SiJi{R) (2.37) 

where (i?) is defined by (|2. 21112.22112. 23b for i? e £, i.e. a polynomial in creation and aimihilation operators. 
Formally, the reduced dynamics in the Heisenberg picture is given by 

Zl{S) [e"^^ (5 ® 1) &-'*^^\ . (2.38) 

However, this definition does not make sense a priori, since e'*^^ {S ® l)e~'*^^ ^ ^(J^) ® £ in general. A 
mathematically precise definition of Zj* is the subject of the upcoming Lemma |Z51 

Since both the initial reservoir state and the Hamiltonian Ii\ are translation-invariant, we expect that the 
reduced evolution is also translation invariant in the sense that 

r.Z^T-. = Z*, where {%S){xi^,xr) := S{xl +z,xr + z) (2.39) 



By the requirement e(fc) = e(— fc) in Assumption 12.11 and the requirement that ^(q) — 4>{—q) in Assumption[ 
the Hamiltonian H\ is also invariant with respect to space-inversion x i— >^ —x, or, equivalently, k i-^ —k. Since the 
initial reservoir state is also invariant with respect to space inversion (this follows from the fact that uj{q) — uj{—q)), 
we expect that 

TeZITe^ - Z*, where iTES){xL,XR) 5(-xl, -a;^) (2.40) 
Finally, the self-adjointness of the time-evolved density matrix implies that 

TjZiTi' = Z*, where {TjS){xl,xr) S*{xr,xl) (2.41) 
where the * in S*{xr, —xl) is just the Hermitian adjoint of operators on .y. 
Lemma 2.5. Assume Assumptions 12 . 1 \\T2\ and \23[ and let 

Ho := Hs + i/R, Hsn{t) := e'*^°iJsRe-'*^° (2.42) 

The Lie-Schwinger series 

ZliS) ^(iA)" J dti...dtn (2.43) 

"^^^ 0< ii < ... < t„ < t 

(ad(ffsR(ii))ad(i/sR(i2)) • ■ • ad{Hsn{tn)) e"-'^^^") (5 ® 1)) 

is well-defined for all A, t e M, that is, the RHS is a norm convergent family of operators and has the following properties 

1) Z*(l) - 1. 

2) TzZ^T-z = with Tz as defined in (|239t . 

3) TeZ*Te^ = Z* with Te as defined in jllOll . 

4) TjZ*Ti^ = Z* with Tj as defined in (|24]) . 

5) ||Z*(5)||oo < \\S\\oo. 

6) Z*(5) > OforS> 

7) For S e ^2(=^)/ the map S Z^{S) is continuous in t in the Hilbert-Schmidt norm \\ ■ \\2- 
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These properties of should not come as a surprise, they hold true trivially if one pretends that the initial 
reservoir state is a density matrix and is obtained by taking the partial trace over the reservoir space, as in 
||1.7| |. One can prove this lemma, under much less restrictive conditions than the stated assumptions, by estimates 
on the RHS. For this purpose, the estimates given in the present paper amply suffice. However, one can also 
define the system-reservoir dynamics as a dynamical system on a von Neumann algebra through the Araki- Woods 
representation and this is the usual approach in the mathematical physics literature, see e.g. ll3l [T2l[T0ll20l l23l. 

We also define Zt : ^i(J^) 3§i{J^%), the reduced dynamics in the Schrodinger representation by duality 

Tr[5Z*(S")] =Tr[Zt(5)5'] (2.44) 
Physically, Z^ is the reduced dynamics on observables of the system and Zt is the reduced dynamics on states. 



2,5 Translation invariance and the fiber decomposition 

In this section, we introduce concepts and notation that will prove useful in the analysis of the reduced evolution 
Zt- These concepts will be used in Section IT2l However, Section |3IT1 which contains the main results, can be 
understood without the concepts introduced in the present section. 
Consider the space of Hilbert-Schmidt operators 

^2{M) ~ ^2(/^(Z'') ® ^) L^{T'^ X T'^,£§2{y),dkLdkR) (2.45) 

and define 

^(fci,fcfl,):= 5(xL,Xi^)e-'("^'=^-"^■'=«^ 5 G ^a^'l^'') ® ^). (2.46) 

XL ,xa£Z'' 

Note the asymmetric normalization of the Fourier transform, which serves to eliminate factors of 2tt in the bulk 
of the paper. In what follows, we will write S for S to keep the notation simple, since the arguments x ^ k will 
indicate whether we are dealing with S or S. To deal conveniently with the translation invariance of our model, 
we change variables, see also Figure [ij 

^^fcL+_fcfl^ P-fci-fcfl, k,peT^ (2.47) 
and, for a.e. p G T'^, we obtain a function Sp G L'^{T'^, S§2{-5^)) by putting 

(5p)(fc) :=5(fc + |,fc-|). (2.48) 

This follows from the fact that the Hilbert space ^2(^3) - iP-^f^ x T'', ^2(^), dfcid/cfl) can be represented as a 
direct integral 



.(JTs) = y dp^P, S = JdpSp, (2.49) 

turally 

by ^, ^' and the scalar product is 



where each 'fiber space' is naturally identified with J# = L^(T'', ^2(^))- Elements of ^ will often be denoted 



{^,0^-^ / dfc Tr^[r(fc)^'(fc)] (2.50) 



jd 

with Tr the trace over the space of internal degrees of freedom S^. 

Let Ti, z G Z"* be the lattice translation defined in ( |2.39t . In momentum space, 

{%S)p = e-'P'Sp, S G ^2(^3). (2.51) 

Since Hx and p^^ are translation invariant, it follows that 

T-.Zt%^Zt. (2.52) 
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Figure 1: The thick black square [—tt, tt] x [— tt, tt] is the momentum space T'* x T'' (drawn here for d = 1), with k^.kji G T**. After changing 
variables to (fc, p) G T"^ x T"^, the momentum space is transformed into the gray rectangle. On sees that the four triangles which lie inside the 
square but ouside the rectangle, are identified with the four triangles inside the rectangle but outside the square. 

Let yy e ^(^2(^)) be translation invariant in the sense that T-z'WTz = W (cf. l l232ll ). Then it follows that, in 
the representation (|2.491 , W acts diagonally in p, i.e. {WS)p depends only on Sp and we define Wp by 

{V\;S)p = yVpSp, Sp e ^, Wp G ) (2.53) 

For the sake of clarity, we give an explicit expression for Wp. Define the kernel 'Wxj^,Xj,;x'^,x'^ by 

{WS){x'^,x'j,)^ J2 W.^,.^;.^.x'„5K,x^), x'L,x'j,eZ''. (2.54) 
Translation invariance is expressed by 

y^Xj^,Xj^-x'^,x'j^ — y^x^+z,Xj^ + z;x'^+z,x'j^+zi Z G Z"*, (2.55) 

and, as an integral kernel, Wp G ^(L^ (T'*, ^2 (^))) is given by 

Wn{k',k)^ V f:M=^L-^R)--^k'(x'^-x'^)^-i^((x'^+x'„)-(x^+x^))y^^ X -x' .X' ■ (2.56) 

x^ = 
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To avoid confusion with other subscripts we will often write 

{S}p instead of and {Wj^ instead of Wp (2.57) 

We also introduce the following transformations. For i/ E T'', let C/^ be the unitary operator acting on the fiber 
spaces as 

iU,Oik)^ak + '^), i^'S (2.58) 
Next, let K = (k^, k^) e x and define the operators by 

{J.S)(x^,x^) := e'^-.-.5(x^,x^)e-'5-H-« (2.59) 

Note that is unbounded if k ^ R'' x R''. 

The relation between the operators J7k and the fiber decomposition is given by the relation 

{XW^-Jp = (2.60) 

^ 4 P 2 4 

as follows from l|2.56|l and the definition (|2.59l l. From l|2.56l l and ll2.60l l, we check that 

p is conjugate to i ((a;'^ + x^) - (xl + xb)) (2-61) 

V is conjugate to {xi^ — xr) — (x^ — x'^). (2-62) 

We state an important lemma on the fiber decomposition. 

Lemma 2.6. Let S e Bi{L^{T'^, y)). Then, Sp is well-defined, for every p, as a function in L^{r'^, ^2(^)) and 

TrJ^S = J2 e~''"'S{x,x) = (l,5p)^, withp = -IL—HR andK = {k^,Kj^) (2.63) 

where 1 stands for the constant function on T'' with value 1 e If, moreover, Jk.S is a Hilbert-Schmidt operator for 

|Im«;^ ^1 < S\ then the function 

T'^^'^ : Sp, (2.64) 

as defined in | |2.48| |, is well-defined for all p e T"* and has a bounded-analytic extension to the strip \lmp\ < 5' . 

The first statement of the lemma follows from the singular-value decomposition for trace-class operators. In 
fact, the correct statement asserts that one can choose Sp such that | |2.63|I holds. Indeed, one can change the value 
of the kernel S{kj^, kj^} on the line kj^ — kj^ = p without changing the operator 5, and hence Sp in (|2.63b can not be 
defined via | |2.48| | in general, if the only condition on S" is G . 

The second statement of Lemma 12.61 is the well-known relation between exponential decay of functions and 
analyticity of their Fourier transforms. Since we will always demand the initial density matrix po to be such that 
I1i7kPo||2 is finite for k in a complex domain, we will mainly need the second statement of Lemma |Z61 

By employing Lemma |Z61 and the properties of Z'l listed in Lemma [2. 5[ it is easy to show that the fimction 

k >^ {Ztpolo (k) e ^{y) (2.65) 
takes values in the positive matrices on and is normalized, i.e., 

J dkTv^[{Ztpo}„ (k)] = {l,{Ztpo}oh = 1 (2.66) 

Further, the space-inversion symmetry and self-adjointness of the density matrix (the third and fourth property in 
Lemma 1231 ) imply that 

E{Zt}^E^{Zt}_^^, where {EO{k) := a-k): for ^ G (2.67) 

J{Zt}pJ = {Zt}_^, where (JO (fc) (?(^))*, for ^ G ^. (2.68) 
where the * on £^{k) is the Hermitian conjugation on ^(^). 
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3 Results 



In this section, we describe our main results. In Section |3A1 we state the results in a direct way, emphasizing the 
physical phenomena. In Section lSTH we describe more general statements that imply all the results stated in Section 
El 

3.1 Diffusion, decoherence and equipartition 

We choose the initial state of the particle to be a density matrix p e 3§i {Jfs ) satisfying 

p > 0, Tr[p] = 1 WJ^ph < oo, (3.1) 

for K in some neighborhood of € x C^. The condition || Jkp||2 < oo reflects the fact that, at time t = 0, the 
particle is exponentially localized near the origin. 

Our results describe the time-evolved density matrix pt := Zfp. Note that pt depends on A, too. First, we state 
that the particle exhibits diffusive motion. 

Define the probability density pt = l^-t' depending on the initial state p E ^^i{Jf%), by 

pt{x) := Tr^ [pt{x, x)] . (3.2) 

It is easy to see that 

^^tix) > 0, /^*(^) = Tr[Pt] = 1- (3.3) 

The following theorem states that the family of probability densities pt{') converges in distribution and in the 
sense of moments to a Gaussian, after rescaling space as a; ^ ^ . 

Theorem 3.1 (Diffusion). Assume Assumptions \2.1l\T2\\2.3\ and \2.4\ Let the initial state p satisfy condition ^3.1^ and let 
Pt be as defined in l|3.2b . 

There is a positive constant Aq such that, for < | A| < Aq, 

^ — ' t /"oo 

with the diffusion matrix Da given by 

Da = A2(D™ + o(A)) (3.5) 

where D^w is the diffusion matrix of the Markovian approximation to our model, to be defined in Section^ Both D\ and Drw 
are strictly positive matrices (i.e., all eigenvalues are strictly positive) with real entries. 

The convergence of pt{-) to a Gaussian also holds in the sense of moments: For any natural number e N, we have 

\ E A**(x)e-^'-^ ) ^-^ (VJV^^-^^^ (3.6) 

In particular, for £ = 2, this means that 

i V a;,Xj/it(x) — > iD\)^,j (3.7) 

Our next result describes the asymptotic 'state' of the particle. Not all observables reach a stationary value as 
t oo. For example, as stated in Theorem 13.11 the position diffuses. The asymptotic state applies to the internal 
degrees of freedom of the particle and to functions of its momentum. Hence, we look at observables of the form 

F(P) (8) A, F = F e L°°(T'^), A = A* e (3.8) 
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with P = P (g) 1 the lattice momentum operator defined in Section I2l2l Such observables can be represented as 
elements of the Hilbert space L'^{T'^) (g) ■S§2{y) ^ i^(T'^, ^2(^)) = (recall that y is finite-dimensional) by the 
obvious mapping 

F[P) ®A^F®A (3.9) 

since L°°(T'^) c L^(T''). Consequently, the asymptotic state is not described by a density matrix on but by a 
positive functional on the Hilbert space This positive functional is called ^'^^ = ^^"^ ('eq' for equilibrium) and we 
identify it with an element of The asymptotic expectation value oi F ® Ais given by 

(F®An)^=/ dfcP(fc)Tr^[r«(fc)A] (3.10) 

We also state a result on decoherence: Equation (|3.13l l expresses that the off-diagonal elements of pt in position 
representation are exponentially damped in the distance from the diagonal. Note that this is not in contradiction 
with Theorem l3.1l as the latter speaks about diagonal elements of pt- 

Theorem 3.2 (Equipartition and decoherence). Assume Assumptions \2.1\ |Z21 \2.3\ and \2.4\ Let the same conditions on 
the coupling constant A and the initial state p he satisfied as in Theorem \3l\ Let A, F be as defined above. Then 



TT[pt{F{P)<»A)]^{F<i^A,e'')^ + 0{e'S^"'), t ^ oo (3.11) 
for some decay rate g > 0. The function f = (^-ff^^ is given by 

r(fc)-^e-'^^ + o(|A|0), forallk^T\ A\0 (3.12) 

with the normalization constant Z{(3) := {2ttY Tr(e~'^^). 

Further, there is a decoherence length (^dch)^^ > such that 

\\pt{x, y)\\m(y) < Ce-^-'-l^-^l + 0(e-f^'*), t / oo (3.13) 

In particular. Theorem 13.21 implies that the inverse decoherence length '^dch remains strictly positive as A \ 0. 
Theorems 13.11 and 13 . 2 1 are derived from more general statements in the next section. 

3.2 Asymptotic fomi of the reduced evolution 

In the following theorem, we present a more general statement about the asymptotic form of the reduced evolution 
Zt ■ The two previous results. Theorems 13.11 and 13. 2[ are in fact immediate consequences of this more general 
statement. 

As argued in Section 12.51 , the operator Zt is translation invariant and hence it can be decomposed along the 
fibers, 

e 

Zt = Jdp {Zt}^ , {Zjp G (3.14) 

The next result. Theorem 13. 3[ lists some long-time properties of the operators {Zt} and U„ {2^t}p U-i, with Ui, 
as defined in l|2.58b . To fix the domains of the parameters p and v, we define 

S)'""" {peT'^ + lT^j/eT^ + iT'' |Rep| <p*, |Imp| <<5, |Imi/| <(5} (3.15) 

S'^^s'' := {peT'^ + iT'*,z/eT'^ + iT'* |Rep| >p72, |Imp| <(5, |Im:/| <5j (3.16) 
depending on some positive constants p* , (5 > 0. 
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Theorem 3.3 (Asymptotic form of reduced evolution). Assume Assumptions 12.21 12.21 12.31 and 12.41 and let the same 
conditions on the coupling constant A and the initial state p be satisfied as in Theorem \3l\ Then there are positive constants 
p* > and (5 > 0, determining the sets ^^^igh ^f^o^g^ gy^/i that the following properties hold: 

1) For small fibers p, i.e., such that {p, 0) e '£)^ow^ there are rank-1 operators P{p, A), bounded operators R''°^{t,p, A) and 
numbers f{p, X), analytic in p on and satisfying 

sup \\U,P{p,X)U_,\\ < C (3.17) 

sup sup\\U^R'°"'{t,p,X)U-^\\ < C (3.18) 

such that 

{Ztjj, = e/(P'^)*P(p, A) + R^'"^{t,p, A)e-(^'9'°™)* (3.19) 

sup Re/(p, A) > -A^g'""' (3.20) 
(p,o)es'°™ 

for a positive rate g^°^ > 0. 

2) For large fibers p, i.e., such that {p, 0) e 'j^high^ there are bounded operators R'"s'^{t,p, X), analytic in p on ^^^9^ ^nd 
satisfying 

sup sup||C/^i?"^'''(i,p,A)?7_^|| = 0(1), A\0 (3.21) 

{p,v)e&^i3h t>o 

and 

{Zjp = R'^^<^\t,p, A)e-(^'^'''"'')*, t /• oo (3.22) 

for some positive rate g^^s^ g_ 

3) The function f{p, A) and rank-1 operator P{p, A) satisfy 

sup |/(p,A)-AV™b)| = odAp) (3.23) 

(p,o)esi'<"" 

sup \\U,P{p,X)U^„-U,P„MU-.\\ = o(|A|°), A\0 (3.24) 

(p,i^)e£)'''"' 

zf/zere the function frwip) and the projection operator Prw{p) ure defined in Section^ 

The main conclusion of this theorem is presented in Figure |2l Let 7?^(z) be the Laplace transform of the reduced 
evolution Zt and {R-iz)}^ its fiber decomposition, i.e.. 



n{z) -.^ J dtc-*'Zt and n{z) = j dp {n{z))^ 

E+ T'' 



(3.25) 



The figure shows the singular points, z = f{p, A), of {'R-{z)}p. Those singular points determine the large time 
asymptotics. If we had not integrated out the reservoirs, i.e., if Zt were the unitary dynamics, then one could 
identify /(p, A) with resonances of the generator of Zt. 

The proof of Theorem 13.31 forms the bulk of the present paper. 

3.3 Connection between Theorem 13.31 and the results in Section |3Jl 

In this section, we show how to derive Theorems 13.11 and 13.21 from Theorem |3.3l 
Since P{p, A) is a rank-l operator, we can write 

P(p,A) = |e(p,A))(|(p,A)|, forsomee(p,A),|(p,A)e^^ (3.26) 
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Figure 2: The singular points of {'R-{z)}p as a function of the fiber momentum p. Above the irregular black line, the only singular points are 
given by /(p, A), in every small fiber p. Below the irregular black lines, we have no control. 



using the notation introduced in (|2.51 . We derive a bound on the eigenvectors £,{p,X) and £_{p,X), analytically 
continued in the coordinate k. This bound follows from the analyticity and uniform boundedness of UvP{j>, X)U-u 
on and straightforward symmetry arguments; 

Lemma 3.4. The vectors UyS,{p, A) and U„^{p, A) can be chosen bounded-analytic on In other words, the operator 

P{p, A) has a kernel 

(3.27) 



p{p,x){k,k')^ apAm) {apA){k') 

which is bounded-analytic in both k and k' in the domain | Im fc | , | Im fc' | < (5 

Note that for fixed fc, k' , the RHS of | |327|| belongs to ^{.^2 (^))- 

Proof. Since P{p, A) is the dominant contribution to {Zt] for large t, the symmetry properties of Zt, stated in Lemma 
12.51 imply in particular that {JE)P{p, X){JE)^^ = P{p, A). Moreover, as JE is an involution, JE = {JE)^^, the 
eigenvectors ^(p, A) and ^(p, A) can be chosen such that JE^{p, A) = ^{p, A) and JE^{p, A) = ^(p, A). Then 



iic/,c(p,A)ii = \\u,jEap,x)\\ = \\jEu-,ap,m = \\u-,ap,m 

Since Ui, = 6^"^^*=, we have also 

||C(P,A)|| < ||2cosh(Im!.VfcK(p,A)|l - |l((7ii„,. + C/-ii„,.K(p, A)|l < 2\\U,^\\, foranyz^e 
The same relation holds for ^{p, A) and hence none of the factors on the RHS of 



|C/.P(p, A)C/-,.||.^(^) = |C/.e(p, X)){U.ap, A)| 



\u.ap,m^\\Uuap,x)y 



can become small as i' varies. The lemma now follows from the uniform boundedness of U,^P{p, A)[/_ 



(3.28) 
(3.29) 

(3.30) 
□ 
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For p — 0, the vectors £,{p, A) and (,{p, A) play a distinguished role, and we rename them as 

C = C ■■= - 0, A), r = er lb - 0, A), (3.31) 

Note that ^^"^ was already referred to in Theorem l3.2l 

By exploiting symmetry and positivity properties of the reduced evolution Zt, we can infer some further prop- 
erties of the function f{p, A) and the operator P(p, A). 

Proposition 3.5. The function f{p, A), defined for all p with {p, 0) € j)'"""^ has a negative real part, Re f{p, A) < 0, and 
satisfies the following properties 

/(p = 0,A) = 0, and V^/b, A)|^^j, = (3.32) 

The Hessian D\ := {Vp)'^f{p, A) has real entries and is strictly positive (3.33) 

The functions ^'^'^ and 1*^' can be chosen such that 

i"i = l, r'^(fc)>0, / dfcTrj^ [r'^(fc)] = (1,^'^) = 1 (3.34) 
where 1 £ ^ is the constant function on with value 1 € ^2 (-5^)- Moreover, it satisfies the space inversion symmetry 

ie") (k) = (e^) i-k). 

The fact that f{p = 0, A) = 0, (^'^'^ = 1 and (|3.341 follow in a straightforward way from (|2.66b and the asymptotic 
form | |3l9l l. The symmetry property C'^9(fc) = ^'"^{-k) and Vp/(p, A) j^^g = follow from iZSTt and | |3l9l l. The fact 

that Da has real entries follows from f{p, A) = f{—p, A) which in turn follows from the reality of the probabilities 
Utix) and the convergence (I3.4I I. 

To derive the strict positivity of D\, we use the claim (in Proposition I4.2|l that Drw, the Hessian of frw{p) at 
p = 0, is strictly positive. By the convergence Il3.23b and the analyticity of frw{p), it follows that \Dx — \^Drw\ \ 
as A \ 0. Indeed, if a sequence of analytic functions is uniformly bounded on some open set and converges 
pointwise on that set, then all derivatives converge as well. 

3.3.1 Diffusion 

We outline the derivation of Theorem l3.1l 

Let p be such that (p, 0) e S'°"'. Then we can calculate the logarithm of the characteristic fxmction: 



log^e-'f>t(x) = log^e-'f-Tr^p*(x,x) 



log(l,{pJp) 
log(l,{Zap{po},) 



= log (e^(^''^)*(l, A) {po}p) + e-^'«'°™*(l, A) {po}p)) 

= loge/(f'^)* ((1, P{p, A) {po},) + e-(^^^"'™-/(''-^))*C||l|| || {po}, ||) 
= /(p, A)t + log ((1, P(p, A) {p„}p) + e-^>^'<^'--fiv,mc\\i\\ II {po}^ 11) (3.35) 



where the scalar product (•, •) and || • || refer to the Hilbert space The second equality follows from Lemma [ 
the fourth from l|3.19l l and the fifth from | |3.18|| . The second term between brackets in the last line vanishes as i 00 
by l|3.20|l . To conclude the calculation, we need to check that the expression in log (•) does not vanish. We note that 

(1, P{p = 0, A) {polo) = (1, r > (r , {polo) = 1 (3.36) 
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as follows from the fact that f ^"^ = 1 and the normalization of f ^"^ in (|3.341 . Hence, for p in a complex neighborhood 
of 0, the expression (1, P{p, A) {po}p) is bounded away from by analyticity in p. Consequently, 

lim i log V e"'P>t(a:) - f{p, A). (3.37) 
tyoo t ^ — ' 

X 

Next, we remark that, for ip real, the LHS of | |3.37|I is a large deviation generating function for the family of measures 
(m* ('))teK+ • ^ classical result [5J in large deviation theory states that the analyticity of the large deviation generating 
function in a neighborhood of implies a central limit theorem for the variable both in distribution, see 1 13.41 , 
as in the sense of moments, see (|3.6b . 

3.3.2 Equipartition 

To derive the result on equipartition in Theorem l3.2[ we consider F, A as in H3.8t . Since pt{F{P) ® A)\sa trace-class 
operator. Lemma |Z61 implies that 

Tr [(F(P) ® A)pt] - (1, {{F{P) (g> A)pt}oh ^ {F (g> A, {ptjoh (3.38) 

where, as in l|3l0l l, F (g) A stands for the function k F{k)A in L'^{T'^, ^2(^)). 
Using Theorem |33] for the fiber p ^ 0, we obtain 

{F ® A, {pt}o) = e/(0'^)*(f^ ® A, Pip = 0, A){po}o) + e-(^'3'°™)'(F ® A, i?'°-(i,p = 0, X){po}o) 

= {F^A,^-'i)+Ce-^^"s"^^'\\F^Ay\\{po}oy (3.39) 

To obtain the second equality, we have used the uniform boundedness of the operators B}°'^{t,p — 0, A) (Statement 
1) of Theorem l3.3b , the fact that /(p = 0, A) = (Proposition l3.5|l and the identities 

P{p = 0, A){po}o - , {po}o)r^ = (1, {po}o>r'' = f " (3.40) 
Hence, from | |3.39|I , we obtain the asymptotic expression l|3.11l l by choosing g < g'"™. 

3.3.3 Decoherence 

In this section, we derive the bound ll3.13l l in Theorem l3.2l We decompose pt as follows, using Theorem l3.3l 



pt J dp {ptjj, (3.41) 

- r dpe^'/(^''^)*P(p,A){po}p+e-^'«"""* dpi?'°"'(t,p,A){po}p (3.42) 

j|p|<p* J\p\<p' 

" V ' " V. ' 

= :Ai =:A2 

+ e-^'s"'''*/® dpi?'"^''(t,p,A){po}p (3.43) 
J\p\>p' 



(3.44) 



The terms A2 and A3 are bounded by 



\A,,3U<C dp\\{po}J'^^C\\poU<C\\po\\i (3.45) 



21 



where the first inequality follows from the bounds l|3.18|l and | |3.22|I . Hence, for our purposes, it suffices to consider 
the first term Ai. To calculate the operator Ai in position representation, we use the kernel expression l|3.27ll for 
P{p, A) to obtain 

{Ai){xl,xh)= [ dpe'5 ("-+"«)e/(P'^)*(e(p,A),{po}J / dfe^b, A)(fc)e"= (^---«) (3.46) 
J\p\<P* -iv 

We now shift the path of integration (in k) into the complex plane, using that the function ^(p, A)(-) is bounded- 
analytic in a strip of width 5. This yields exponential decay in [xl — xu)- Using also that Re f{p, A) < 0, for \p\ < p* 
(see Proposition |3.5t , we obtain the bound 

||(Ai)(a:i,X;,)||.^,(^) <Ce-^l----«l, for7<5 (3.47) 

Combining the bounds on Ai and A2, A3, we obtain 

||pt(a:i,x^)|U,(^) < Ce-^l-----l +C'e-(^'s)*, for 7 < 5 (3.48) 

with q := min 

(^giow ^ ghtghy jj^g fgj,^ ^^^^^ ^^^^g bouud Is vaHd for any 7 < (5, confirms the claim that the inverse 
decoherence length ^dch can be chosen uniformly in A as A \ 0. 



4 The Markov approximation 

For small coupling strength A and times of order A~^, one can approximate the reduced evolution Zt by a "quantum 
Markov semigroup" Af which is of the form 

At = e*(-'"'*(^)+^'-^) (4.1) 

where Y = 1 (E) Y is the Hamiltonian of the internal degrees of freedom, and is a Lindblad generator, see 
e.g. |1J. Lindblad generators, and especially the semigroups they generate, have received a lot of attention lately in 
quantum information theory. The operator A4 has the additional property of being translation-invariant. Translation- 
invariant Lindbladians have been classified in [21 1 and, recently, studied in a physical context; see |35| for a review. 
In Section 14.11 we construct Ai and we state its relation with . We also describe heuristically how A4 emerges 
from time-dependent perturbation theory in A as a lowest order approximation to Zf. In Section l4!2l we discuss 
the momentum representation of Ai (the derivation of this representation is however deferred to Appendix O , 
and we recognise that the evolution equation generated by is a mixture of a linear Boltzmann equation for the 
translational degrees of freedom and a Pauli master equation for the internal degrees of freedom. In Section l4!3l 
we discuss spectral properties of Ai, which are largely proven in Appendix |C] Finally, in Section l4.3.1[ we derive 
bounds on the long-time behavior of Atp, for any density matrix p e ^i(J^). 

4.1 Construction of the semigroup 

First, we define the operator C*{t) on 3§{J^s)' 

C*it){S) = -pi (ad(i?sR) e'*^'i(^+^-) adiHsR){S ® 1)) (4.2) 

This definition makes sense since the conditional expectation p^ is applied to an element of ^^{J^ ) (?) £, see Section 
12.41 Then we consider the Laplace transform of C*{t), i.e. 

C*{z) = J dte-'' tit), Rez>0, (4.3) 

R+ 
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and, finally, we let C{z) be the dual operator to C*{z), acting on i^i(J??g), see (|2.441 . Then the operator M is 
obtained from C by "spectral averaging" and adding the "Hamiltonian" term — iad(e(P)): 

M -iad(e(P)) + ^ la(ad(r))£(-ia)la(ad(y)) (4.4) 

aesp(ad(y)) 

For now, this definition is formal, since it involves (|4.3l l with Re z = 0. 

The following proposition provides a careful definition of M and collects some basic properties of the semig- 
roup evolution At. 

Proposition 4.1. Assume Assumptions \2.1l |Z2] and |Z31 Then, the operators C{z), defined above, can be continued from 
Rez > Oto a continuous function in the region Re z > and 

sup \\J^C{z)J^^\\ < oo, for keC^xC^ (4.5) 

Rcz>0 

(In fact, J^C{z)J^, = £{z)). 

The operator M, as defined in ||4.4| |, is bounded both on ^i{J^) and ^2(=^)- Recall the constants qdl), 7 > 0, defined 
in Assumption \2.1\ Then 

WJ^MJ.^ - M\\ < (ZedlmKil) + gedlm^^l), \1^^^l,r\ < (4.6) 

where the norm \\ ■ || refers to the operator norm on ^{^2i'^))- 
The family of operators At, defined in l|4.1b , 

A^^e'^'i-d(Y)+x'M)^ ieM+ (4.7) 
is a "quantum dynamical semigroup" . This mean^ 

i) At^At^ = Ati+t2 for allti,t2>Q (semigroup property) 

a) Atp > /or flwyO < p e ^i(J^) (positivity preservation) (4.8) 

Hi) Tr Afp = Trp for any Q < p £ (trace preservation) 

We postpone the proof of this proposition to Appendix O 
4.1.1 Motivation of the semigroup Af 

The connection of the semigroup Af with the reduced evolution Zt is that, for any T < oo, 

sup ||Zf ~Af||^(,^,(^3)) -o(A"), A\0 (4.9) 

0<t<A-2T 

Results in the spirit of (|4.91 have been advocated by [22^ and first proven, for confined (i.e. with no translational 
degrees of freedom) systems, in |9| . They go under the name "weak coupling limit" and they have given rise to 
extended mathematical studies, see e.g. Il29l [TTI . In our model, (|4.91 will be implied by our proofs but we will 
not state it explicitly in the form given above. In fact, statements like (|4.91 can be proven under much weaker 
assumptions than those in our model; see 1133 1 for a proof which holds in all dimensions d> I. 

Here, we restrict ourselves to a short and heuristic sketch of the way the Lindblad generator A4 emerges from 
the full dynamics. First, we consider the Lie-Schwinger series (|2.43t in the interaction picture with respect to the 
free internal degrees of freedom, i.e. we consider 2^*e~'*'*'^^^' instead of Z^. Keeping only terms up to second order 
in A in | |2.43l l and substituting our definition for £* (t) we obtain 

2^*g-itad(y) ^ 1+ iA2ted(e(P))+ A2 J diidt2e'(*-*^^"'*(^)£*(t2-ii)e"'^*-*^)"'^^^^+0(A^) (4.10) 

0<tl<t2<t 

^Most authors include "complete positivity" as a property of quantum dynamical semigroups, see e.g. fTl. Although the operators At satisfy 
complete positivity, we do not stress this since it is not important for our analysis. 
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where we have also used [s{P), Y] = to simplify the second term on the RHS. It is useful to rewrite the third 
term by splitting ^2=^1 + (^2 - ^i) and inserting the spectral decomposition of unity corresponding to the operator 

ad(r) 

J2 f dtie'(*-*i)("-"') la(ad(r)) ( /" ' due~''"'C*{u)]la'{sid{Y)) (4.11) 

a,a' ^sp{ad{Y)) 

Next, we analyze the RHS of (|4.10l l for long times; we choose t — X^^t, and we argue that, in the limit A \ 0, 
the three terms reduce to the operator 

1 + iad(e(P))t + U(ad(y))£*(ia)l,(ad(y))^ t (4.12) 

This limit can be straightforwardly justified if the function t — > £* (t) is (norm-)integrable, which will follow from 
Assumptions 12.21 and 12.31 in our case. Indeed, if < — > £*{t) is integrable, then the integral ' *^ du . . . in (|4.11)l 
converges to du . . . for fixed ti, as A — > 0. This yields the Laplace transform C*{ia). The restriction a = a' 
appears then because 

A2 / dtie'(*-*i'("-"') ^ t6aa' (4.13) 

Jo -^^0 

and one finishes the argument by invoking dominated convergence. Comparing with (|4.4| | and using that (1^ (ad(y)))* 

l_a(ad(y)), one checks that l|4.12|l is equal to 1 + iA4*, with A4* the dual to A4. This is the beginning of a series 

defining the semigroup c'-^ (we got only the first two terms because we kept only terms of order A" and A^ in the 

original Lie-Schwinger series). 

4,2 Momentum space representation of Ai 

In this section, we give an explicit and intuitive expression for the operator Ai. As Ai is translation covariant, i.e., 
TzAiT-z = Ai, as in (|2.52l l, we have the fiber decomposition. 



M= dpAip (4.14) 



where the notation is as introduced in Section 1231 We describe A4p explicitly as an operator on It is of the form 

(XpO(fc) = m] - + f ) - e{k - |))e(fc) + e e ^ (4.15) 



where e is the dispersion law of the particle, see Section |Z21 and T is a self-adjoint matrix in whose only 

relevant property is that it commutes with Y, i.e., [Y, T] = 0. Physically, it describes the Lamb-shift of the internal 
degrees of freedom due to the coupling to the reservoir and its explicit form is given in Appendix [C] The operator 
Af is given, for ^ G C{T^, ■S§2{y)), by 

(AAO(fc)= j dk' (^ra{k\k)Waak')W:-]^Ta{k,k'){mW:Wa+W:Wam)) (4.16) 

with the (singular) jump rates 

r { i_c-V(,) '^(^(g) - a)5fd {k-k' - q) a > 

ra{k,k'):=2^ dq\<j,{q)\^{ (4.17) 

{ ^F:^^5{uj{q) + a)5^d{k + q - k') a<0 

where (f) is the form-factor, see Section |231 and 5fd{-) is the Dirac delta distribution on the torus, i.e.; 

5^d{-):= Yl S{--qo), (4.18) 

go=0+(27rZ)<' 
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Note that ra{k,k') vanishes at a = 0, due to the fact that the 'effictive squared form factor' ?/'(•) vanishes at 0, see 
Assumption 12.31 

Equation (|4.15l l is most easily checked starting from the expressions for M. in Section IC.ll In particular, the 
three terms in (|4.16b correspond to the fiber decompositions of the operators — ■^<l>*(l)/9, — ip(f>*(l) in | |C.7| |, 
and the first two terms on the RHS of | |4.15l l correspond to the commutator with T and e{P) in | |C.7| |. 

We already stated that M is translation-invariant, hence it commutes with ad(P). However, the operator M. 
also commutes with ad(F), as can be easily checked starting from the expressions (|4.151 and | |4.16|I and employing 
the definitions of Wa in ^TM and the fact that \Y, T] = 0. 

We can therefore construct the double decomposition 

® / dpTWp.a (4.19) 

aesp(ad(y)) Jfd 

where 

Mp.a la(ad(r))Xpl,(ad(y)) (4.20) 



To proceed, we make use of our strong nondegeneracy condition in Assumption 12.41 Indeed, the operators AAp^a 
act on functions ^ G that satisfy the constraint 

e(fc) = la(ad(y))e(fc) = le{Y)m^e'{Y), ^ ^ ~ (T'i, (^)) (4.21) 

e,e'(^spy,e— e'— a 

Due to the non-degeneracy assumption, the sum on the RHS contains only one non-zero term for a 7^ 0, i.e., there 
are unique eigenvalues e, e' such that a = e — e' . Let us denote the eigenvector in the space ,5^ of the operator Y 
with eigenvalue e by e as well (cf r. the discussion following Assumption 12.41 , then this unique term in (|4.211 can 
be written as 

le(r)e(fc)le.(y)-((e,^(fc)e')) |e)(e'|, e - e' = a (4.22) 
It follows that the matrix valued function ^(fc) satisfying (|4.21) l can be identified with the C-valued function 

^{k) = {e,m^').y (4.23) 

For a — Q,a function ^(fc) satisfying (|4.21) l is necessarily diagonal in the basis of eigenvectors of Y . In that case, we 
can identify ^ with 

^{k,e) = {e,p{k,k)e)y (4.24) 

Hence, we can identify A^p,a^o with an operator on L'^iT'^) and M.pfl with an operator on L'^iT'^ x spF). 

A careful analysis of these operators is performed in AppendixO Here, we discuss the operator A^o.o because 
it is crucial for understanding our model. 

4.2.1 The Markov generator A^o,o 

Let us choose ip E C{T'^ x spY). Then, by the formulas given above, the operator Alo,o acts as 

MooV{k,e):= [ dk' V {r{k' ,e';k,e)ip{k' ,e') - r{k,e;k' ,e')ip{k,e)) (4.25) 

where r{k, e; k' , e') are (singular) transition rates given explicitly by 

r{k,e;k',e') := r,,_^, {k, k')\{e' ,We)\^ (4.26) 

In formula l|4.25l l, one recognizes the structure of a Markov generator, acting on absolutely continuous probability 
measures (hence -functions) on T'' x spY. The numbers 

j(e,fc) := / dfc' r{k,e-k',e') (4.27) 
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are called escape rates in the context of Markov processes. Let = Jj^ '^f^J2eespY Ifi^^^)]- Then one easily 
checks by decomposing ip = (p+ + and using that r{k, e; k' , e')dk' is a positive measure, that 

||Xo,o¥'lli < 4 |^supj(e,fc)^ \Mu (4.28) 

which implies that A4o,o is bounded on L^(T^ x spF). In particular, this means that A^q.o is a bonafide Markov 
generator (i.e. it generates a strongly continuous (in our case even norm-continuous) semigroup) and e*^°'"(p is a 
probability measure for all t > 0. Physically speaking, the probability density ip is read off from the diagonal part 
of the density matrix p, see l|4.24|l . 

It is interesting to see that the transition rates r{k,e;k',e') satisfy detailed balance at inverse temperature /? for 
the internal energy levels e, e', and at infinite temperature for the momenta k,k'. 

r{k, e; k' , e') = e^('=-'='V(/fc/ e'; k, e) (4.29) 

Physically, we would expect overall detailed balance at inverse temperature /?, i.e. 

r{k, e; fc', e') = e'^(^('=''=)-^('='''='»r(fc/ e'; k, e) (4.30) 

where the energy E{k, e) should depend on both e and k. To understand why E does not depend on k in | |4.29l l, we 
recall that the kinetic energy of the particle is assumed to be of order A^; hence, the total energy is e + A^e(fc) which 
reduces to e in zeroth order in A. 

One can associate an intuitive picture with the operator Alo.o- It describes the stochastic evolution of a particle 
with momentum k and energy e. The state of the particle changes from (fc, e) to (fc', e') by emitting and absorbing 
reservoir particles with momentum q and energy uj{q), such that total momentum and total energy (which does 
not include any contribution from k, k') are conserved, see Figure|3] 




1 q 

(a) Emission of a boson, k = k' + q and e = uj{q) + e'. (b) Absorption of a boson, k + q = k' and uj{q) + e = e'. 

Figure 3: The processes contributing to the gain term (the first term on the RHS in )4.25) ) of the operator A4o,o- Emission corresponds to 
e > e' and absorption to e < e'. 

It is clear from the collision rules in Figure|3]that, in the absence of internal degrees of freedom, the particle can 
only emit or absorb bosons with momentum q — 0, and hence it cannot change its momentum. This means that 
without the internal degrees of freedom, the semigroup Af would not exhibit any diffusive motion. This is indeed 
the reason why we introduced these internal degrees of freedom. 

4.3 Asymptotic properties of the semigroup 

The following Proposition 14.21 states some spectral results on the Lindblad operator A4 and its restriction to mo- 
mentum fibers Aip £ ^(^). These results are stated in a way that mirrors, as closely as possible, the statements of 
Theorem l3.3l 

These results are useful for two purposes. First of all, they show that our main physical results. Theorems 
13.11 and 13.21 hold true if one replaces the reduced evolution Zt by the semigroup At (see the remark following 
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Proposition 132). Second, a bound which follows directly from Proposition 14.21 will be a crucial ingredient in the 
proof of our main result Theorem 13. 31 This bound is stated in l|4.551 in Section l4.3.1l 
We introduce the following sets (cf. l l3.15H3.16|l ) 

:= {p e T"^ + ir', e T'* + IT-^I |Rep| < p*,^,, \lmp\ < Sr^, |Imi/| < (4.31) 
^rl"" ■■= jpe T'^ + iT^j. e T'^ + iT'^l |Rep| > ^p^^, \lmp\ < <5™, |Imz.| < (5,^| , (4.32) 

depending on positive parameters p*^, > and 6rw > 0. The subscript 'rw' stands for 'random walk' and it will 
always be used for objects related to A^. 

Proposition 4.2. Assume Assumptions \2.1\\T2\\2.3\ and \2.4\ There are positive constants p*^ > OandSrw > 0, determining 
1 23™ above, such that the following properties hold 

1) For small fibers p, i.e., such that {p, v) £ T)'"^', the operator UyMpU-^ is bounded and has a simple eigenvalue frw{p), 
independent ofv, 

sp(C/,XpC/_,) = U (4.33) 

The eigenvalue frwip) elevated above the rest of the spectrum, uniformly in p, i.e., there is a positive gl°^ > such 
that 

Re rip,, < ~giZ < Re fru, (p) < 0, for all {p, v) e T>iZ (4.34) 
The one-dimensional spectral projector UuPrw{j))U-u corresponding to the eigenvalue frwip) satisfies 

sup \\U,PrUp)U-.\\ < C (4.35) 
(p,i/)es'°"' 

2) For large fibers p, i.e., such that (p, 0) e l^^w^' operator UuMpU-u is bounded and its spectrum lies entirely below 
the real axis, i.e., 

sup Resp([/,XpC/_j.) < -g^;^'', for some g'^'j'' > (4.36) 

3) The function frw{p), defined for all p such that (p, 0) G S^-iT' ^ negative real part, Re frwip) < 0, and satisfies 

frwip ^0)^0, and Vj,/™(p))|^^o = (4.37) 

The Hessian Drw '■— i'^p)'^ frwip)\p_Q has real entries and is strictly positive (4.38) 
The spectral projector Prwip — 0) is given by 

Prwip = 0) = \erl){Cll (4.39) 

with 

1 e-f^^ 

= and e^^(fc) = keT^ (4.40) 

The conclusion of Proposition 14.21 is sketched in Figure H) The proof of this proposition is very analogous to 
|7| (which, however, does not consider internal degrees of freedom). For completeness, we reproduce the proof in 
Appendix ICl 

From Proposition 14.21 one can derive that the semigroup e*(^'^'^(^)+-^) exhibits diffusion, decoherence and 
equipartition. This follows by analogous reasoning as in Sections l3.3.1[|33l2l and l3l3.3l but starting from Proposition 
14.21 instead of Theorem 13.31 The matrix Drw is the diffusion constant, the inverse decoherence length has to be 
chosen smaller than 5rw and the function ^^2; is the 'equilibrium state'. We do not derive these properties explicitly 
as they are not necessary for the proof of our main results. 
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spectrum of Mp (real part) 



fiber p 




Figure 4: The spectrum of A4p as a function of the fiber momentum p. Above the irregular black line, the only spectrum consists of the 
isolated eigenvalue frw (p), in every small fiber p. Below the irregular black lines, we have no control. 



4.3.1 Bound on At in position representation 

By virtue of Proposition l4.2l we can write 

{AJp = P™(p)e^'/'-™(f)*+<7(i,p)e-^'3™*, (p,^^)eS)tr (4.41) 

{AJp = i?«3"(p,t)e-^'f™'S (4.42) 
with Prw (p) as defined above, satisfying (|4.351 , and the operators Rl°^ , i?^™ satisfying 



:flow 

(v,i^)&'S\?^ *>0 



sup sn^\\U,RZ{t,p)U-u\\<C (4.43) 



sup sup 1 1 [/.i?^?;^'' {t,p)U-,\\<C (4.44) 

The appearance of the factor is due to the fact that multiplies M in the definition of the semigroup Kt- 

Next, we derive estimates on Ai (see ee.g. the bound ( |4.55t below) starting from ll4.41H4.42t and | |4.43H4.44t . 

without using explicitly the semigroup property of Af. This is important since in the proof of Lemma [6.21 we will 

carry out an analogous derivation for objects which are not semigroups. 

We choose n = (k^, k^) e x such that Re = Re and we calculate, using relation H2.60I I, 

J.^tJ-. = j ApU, {Kt}p+^p with Ap ^L^-^ v ^fi^ (4 45) 

where we are using the analyticity in (p, v), see ( |4.35t and | |4.43H4.44t . Recall that {Ajj, acts on "i^^ ^ ^ 
//^(T"^, ^(.5^)). Our choice for k ensures that Ap and v are purely imaginary. 

Next, we split the integration over p e into small fibers (|p| < P™,) and large fibers (|p| > P™) t'y defining 

{p + Ap|peTMp|<p;„}, /'"s":={p + Ap|peTMp|>p;^} (4.46) 
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Using the relations (|4.411 and l|4.42|l , we obtain 

e ® 
J.^tJ-. - J dpe^"f'-^P^*U,PrUp)U-, + e-^'s™ * y dp[/,i?5,T(p,0C/-. (4.47) 



= =:B2 



y dpc/,i??;f''(p,t)[/-. 



Jhigh 



We establish decay properties of the operators Si. 2, 3 in position representation. For B2 and B3 we proceed as 
follows. Recall the duality (|2. 61112. 62|l . By varying p and ly, and using the bounds (|4.43II4.44|| , we obtain 

\\iB2,3)x^,x^;x'^,x'Jy < Ce-*|('='^+'^'«)"('^^+'^«^le-''l(^^-^«'^(^'^"^«)| (4.48) 

for any 7 < 6rw 

For Bi, we need a better bound, which is attained by exploiting the fact that Prw{p) is a rank-1 operator with a 
kernel of the form (recall the notation of l|2.51 ) 

P™(p)(fc,fc')= |(C™(p))(fc))((^™(p))(fc')|, forsome e ^ (4.49) 

where both ^riu(?'), ^rtud?) are bounded-analytic functions of fc, k', respectively, in a strip of width Srw- This follows 
from boundedness and analyticity of Prw{p) by the same reasoning as in Lemma 13.41 By the definition of Bi in 
gliTll and (l236t . (g^Qj, 

(Bi):,^,^^;:ri,,x^ ^ I dpe^'-''™(p)*e-'§(("'^+^«'"("^^+"^«)) (4.50) 



dfc / dfc'e~"=(^^^^«)+"='(^'^-^^)P™(p)(fc,fc') (4.51) 

By the analyticity of Prw{p){-, •) in both fc and k', we derive, for 7 < (5riu, 

< Ce'''-('^'^)*e-*l('=^+'^«^"('=^+'=^^le-Tl^^-^«le-Tl^'^-^«l (4.52) 
where the function rrw{j, A) is given by 

r-ru;(7: A) := sup Re/™(p). (4.53) 

|Imp|<7,|Rcp|<p*^ 

Note that 

r™(7,A) :=0(A2)0(72), A\0,7\0, (4.54) 

as follows from Re frw{p) < 0. The bound (|4.54t will be used to argue that the exponential blowup in time, given 
]-,y gr^„(7,A)t compensated by the decay e^f by choosing 7 small enough, see Section|5l 

Putting the bounds on -81,2.3 together, we arrive at 

+ C"e~'^^^™*e~^ l^'^'^'^^^^^^'-'^^^'^^^le"''!*^'^^^'^"^^*^'^'^"'^^-'! (4.55) 

for 7 < 5rw and with the rate grw ■= min {g\°w i9rw^)- The bound | |4.55|I is the main result of the present section and 
it will be used in a crucial way in the proofs. The importance of this bound is explained in Section lSTil 
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For completeness, we note that a bound like 

ll(AO.,,.«;.i,xgi <e^^''?^(-^)*e-^l^'--^-le-^l^'«-^-l (4.56) 
can be derived simply from the fact that JkMJk is bounded for complex k, see (|4.6t , since 



is conjugate to (x^ — x^) (4-57) 

is conjugate to [x'^ — x^), (4.58) 



5 Strategy of the proofs 



JX 

^SR ■ — 



In this section, we outline our strategy for proving the results in Section |31 We start by introducing and analyzing 
the space-time reservoir correlation function V'(a:, t). Then we introduce a perturbation expansion for the reduced 
evolution Zt (which involves the reservoir correlation function). Afterwards, we describe and motivate the tem- 
poral cutoff that we will put into the expansion. Finally, a plan of the proof is given. 

5.1 Reservoir correlation function 

A quantity that will play an important role in our analysis is the free reservoir correlation function ip{x, t), which 
we define next. Let 

dq ((/.(<z)e'«-l, ® a, + 0(^e-"?-l, ® a*^ (5.1) 

where 1^ — ^x{X) is the projection on J^, acting as {lxip){x') = 5x,x'V{x) for e P{Z'^, y). The operator /gp is 
the part of the system-reservoir coupling that acts at site x after setting the matrix W £ equal to 1 (recall 

that the matrix W describes the coupling of the internal degrees of freedom to the reservoir). We also define the 
time-evolved interaction term, with the time-evolution given by the free reservoir dynamics 

/sR(i) e"^-/s-Re-"«- (5.2) 
The reservoir correlation function ^ is then defined as 

i,{x,t) pi [/s\(iKsR(0)] , 

= (0%T^e"-0)„ + (0%(l + T;3)e-'*"0)^ 

dwVi(a;)e'"* ^ dse'""'^ (5.3) 

where (0^)((7) :— e^'^'^<j){q) and S''^^ is the d — 1-dimensional hypersphere of unit radius. The 'effective squared 
form factor' was introduced in i2.27} , and the density operator in Section l2.3.3l 

Assumptions l2.2] and l2.3l imply certain properties of the correlation function that will be primary ingredients of 
the proofs. We state these properties as lemmata. In fact, one could treat these properties as the very assumptions 
of our paper, since, in practice. Assumptions 12. 2] and l2.3l will only be used to guarantee these properties. Lemmata 
I5.1l and [5.2l The straightforward proofs of Lemmata l5.1l and l5.2l are postponed to Appendix lAl 

The following lemma states that the free reservoir has exponential decay in t whenever \x\/t is smaller than 
some speed . 

Lemma 5.1 (Exponential decay at 'subluminal' speed). Assume Assumptions \2.2\ and \2.3\ Then there are positive con- 
stants > 0, .gn, > such that 

\x\ 

|V'(a;,i)| < Cexp(— if -j- <v^., for some constant C . (5.4) 
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Property (|5.41 is satisfied if the reservoir is 'relativistic', i.e., if the dispersion law uj{q) of the reservoir particles 
is linear in the momentum \q\, temperature (3^^ is positive and the form factor (j> satisfies the infrared regularity 
condition that k H> |(/)(A;)p|/c| is analytic in a strip around the real axis. The speed v* has to be chosen strictly 
smaller than the propagation speed of the reservoir modes given by the slope of cj. In fact, the decay rate 
vanishes when v* approaches the propagation speed of the reservoir modes. Lemma [5. II does not depend on the 
fact that the dimension d > 4. 

Lemma 5.2 (Time-integrable correlations). Assume Assumptions \2.2\ and \2.3\ Then 

/ d< sup |i/'(a;,i)| < oo (5.5) 

This property is satisfied for non-relativistic reservoirs, with uj{q) cx in d > 3 and for relativistic reservoirs, 
with uj{q) oc \q\, in d > 4, provided that we choose the coupling to be sufficiently regular in the infrared. 

□ 



5.2 The Dyson expansion 

In this section, we set up a convenient notation to handle the Dyson expansion introduced in Lemma |231 
We define the group Uf on ^(^) by 

and the operators Ix.i, with x G Z'^ and / e {L, R} {L, R stand for 'left' and 'right'), as 



~i {U®W)S 
i S(lx®W) 



if l = L 
if l^R 



S e.^{.M;)- 



(5.6) 



(5.7) 



where the operators = lx{X) are projections on a lattice site x e Z'*, as used in Section |5A] 

We write {ti,Xi,li),i = 1, . . . , 2ri to denote 2n triples inRxZ'^x{L,i?} and we assume them to be ordered by the 
time coordinates, i.e., ti < t^+i. We evaluate the Lie-Schwinger series (|2.431 using the properties (|2. 21112.22112. 23t , 
and we arrive at 



E E 



d<i . . . dt2nC(^, iU,X,, k)i2l)Vt{{U,X,,h)i2l) 



0<ti<...<t2„<t 



where tt e P„ are pairings, as in l|2.23b , and we define 



Vt{{ti,Xi,li 



■ ■^X2:l2^t2-ti^Xi,lii^ti 



with Ut as in 



and 



(5.8) 



(5.9) 



(r,s)G7r 



^/j{Xs 


X J* • ~^ s 


-u) 


^ 


= L 




" 3^ ■f • i ^ 


-tr) 


L/r. — ig 


= R 






~tr) 


If — Zy, Is 


= R 


■>P{Xs 






If — R-i Is 


= L 



(5.10) 



with the correlation function ?/; as defined in I l5.3b . We recall the convention r < s for each element of a pairing vr. 
For n = 0, the integral in l|5.8|l is meant to be equal to Ut- In Section [71 we will introduce some combinatorial 
concepts to deal with the pairings it £ Vn that are used in < |5.8l l. For convenience, we will replace the variables 
(tt, {ti,Xi,li)^I^i) G P„ X ([0, t] X Z'' X {L, i?})^" by a single variable a that carries the same information. 
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(tl,Kl,il) it2,X2,l2) (t4,X4,U) (ts.aJs./s) itv,XT,lr) (tg.Xg.lg) (t lo ,Xio , Jlo) (* 12 .2:12 I2) 




(t3,X3,l3) (te,xe,ls) (ts.xs,ls) (tn 



Figure 5: Graphical representation of a term contributing to the RHS of tSlSl with n = {(1,3), (2,4), (5,8), (6, 10), (7, 11), (9, 12)} e Ve- The 
times ti correspond to the position of the points on the horizontal axis. 

Starting from this graphical representation, we can reconstruct the corresponding term in )5.8> - an operator on .^2{-^))- as follows 

• Toeachstraightline between the points (ti, aii, /i) and (tj , aij, ij), one associates theoperators e^''*j ^'''^s^ with itbeing — for = Ij = L 
and + for li = Ij = R. 

• To each point (ti, Xi, li), one associates the operator A^I^,. defined in 15.71 . 

• To each curved line between the points (tr, x,., Ir) and {tg,Xs, h), with r < s, we associate the factor ip'f^ixs — Xr,ts — tr) with ■0* being 
ip or ij), depending on i, , 's/ as prescribed in IS.lOt . 

Rules like these are commonly called "Feynman rules" by physicists. 

5.3 The cut-off model 

In our model, the space-time correlation function t) does not decay exponentially in time, uniformly in space, 
i.e, 

there is no g>0 such that snp \ip{x,t)\ < Ce~^^*^ (5-11) 

xez.'^ 

The impossibility of choosing the form factor </> or any other model parameter such that one has exponential decay 
is a fundamental consequence of local momentum conservation, as explained in Section [Ol 

However, if the correlation function ip{x,t) did decay exponentially, we could set up a perturbation expan- 
sion for Zt around the Markovian limit A(. Such a scheme was implemented in [321 , building on an expansion 
introduced in |31|. 

In the present section, we modify our model by introducing a cutoff time t into the correlation fimction ip{x, t). 
More concretely, we modify the perturbation expansion for Zt by replacing 

%l^{x,t), — > l\t\<ri^{x,t) (5.12) 

The cutoff time t will be chosen as a function of A satisfying 

t(A) ^ oo, t(A)A ^0, as A \ (5.13) 

However, we will take care to keep r explicit in the estimates, until Section [8] where the r-dependence will often 
be hidden in generic constants c(7, A), c'(7, A). With the cut-off in place, the correlation function ^(x, t) decays 
exponentially, uniformly in a; e l/\ i.e., obviously, 

sup l|t|<,|V^(a;,i)| < Ce-^. (5.14) 
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The modified reduced dynamics obtained in this way will be called Z/" . 
That is; 

^t^Y^Y. J dii...dt2„Cr(7r,(t„x„?,)'=i)Vt((i„a;„/,)?^i) (5.15) 

n€Z+7reV„ o<ti <. . .<t2„ <* 

with 

Cr{7T,{U,X,,k)tl):= i n Mt.-tA<r]ci^At^,^^,h)tl) (5.16) 

If r is chosen to be independent of A then one can analyze by the technique deployed in 1321 . It turns out 
that for a A-dependent r, one can still analyze the cutoff model by the same techniques as long as A^r(A) \ as 
A \ 0, which is satisfied by our choice (|5.131 . The analysis of Z^ is outlined in Lemma [6. II in Section l6A] The main 
conclusion of the treatment of the cutoff model is that 



The cutoff reduced dynamics Z^ is 'close' to the semigroup At (5.17) 

This conclusion is partially embodied in Lemma 16.21 For the sake of this explanatory chapter, one can identify ZJ 
with A^. 

The reason why it is useful to treat the cutoff model first, is that we will perform a renormalization step, effect- 
ively replacing the free evolution Ut in the expansion (|5.8b by the cutoff reduced dynamics Z^ . The benefit of such 
a replacement is explained in Section lSTil 



5.4 Exponential decay for the renormalized correlation function 
5.4.1 The joint system-reservoir correlation function 

We recall that the free reservoir correlation function ijj{x, t) does not decay exponentially in t, uniformly in x. This 
was mentioned already in Section I5l3l and it motivated the introduction of the temporal cutoff r. 

In the perturbation expansion for the reduced evolution Zt, the correlation function ?/'(a;, t) models the propaga- 
tion of reservoir modes over a space-time 'distance' (x, t) and it occurs together with terms describing the propaga- 
tion of the particle. Let us look at the lowest-order terms in the expansion of Zt, introduced in Section lS^ above: 



Zt^y 



/ <lt2<lti ^ ip*{x2- xi,t2-ti) Ut-t2lx2,h^t2-tilxi,hl^ti +i^is'^er orders in X 



(5.18) 



0<tl<t2<t 



with ip"^ being ip or ip, as prescribed by the rules in l|5.10|l . It is natural to ask whether the 'joint correlation function' 

■lp*ix2 - Xi,t2 -ti) Ix^^iJAt2-ttIxt,h (5-19) 

has better decay properties than t) by itself. In particular, we ask whether | |5.19l l is exponentially decaying in 
t2 — ti, uniformly in a;2 — xi. This turns out to be the case only if li = I2 since in that case, the question essentially 
amounts to bounding 



\lp{x2 - Xi,t2 - ti)\ X 



,±i(i2-ti)A^£(P) 



ixi,X2] 



(5.20) 



The expression Il5.20b has exponential decay in time because 



For speed 



t2-tl 



greater than some v > 0, we estimate 



< e-(7"-A"ge(7))l*2-ti| 



for < 7 < 



(5.21) 



with 6e,qe{-) as in ||2.12| | and Assumption l2.1l Hence, for fixed v, one can choose 7 so as to make the exponent 
on the RHS of l|5.2H l negative, for A small enough. 
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• For speeds ^t\Jll smaller than w* > 0, the reservoir correlation function (0:2 — a;i, t2 — ^i) decays with rate 

gn, as asserted in Lemma [5 . 1 1 with v* as defined therein. 
When 7^ ^2 in (|5.191 , there is no decay at all fromZ^f^-f^, in other words, the 'matrix element' 

(5-22) 



is obviously not decaying in the variables Xj^ — x'j^ or — x^, since it is a function of Xj^ 



x'j^ and x^ 



■c'j^ only. 



Hence, for li 7^ l\, the joint correlation function (|5.191 has as poor decay properties as the reservoir correlation 
function i\)(x^ t). 

The situation is summarized in the following table 



joint S — R correlation fct. 


\x\/t > V* 


\x\/t < V* 




ll=l2 


h, I2 arbitrary 


No exp. decay 


exp. decay from 


exp. decay from 

Tp^X.t) 



5.4.2 Renormalized joint correlation function 

The bad decay property of the joint correlation function l|5.191 suggests to perform a renormalization step, replacing 
the free propagator Ut by the cutoff reduced dynamics , for which (|5.19l l has exponential decay when li ^ I2. 
The cutoff reduced dynamics was introduced in Section 1531 where we argued that it is well approximated by 
the Markov semigroup Aj. Hence, we replace the group Ut by the semigroup At in ( |5.19t , thus obtaining a 'renor- 
malized joint system-reservoir correlation function'. We then check that the so-defined renormalized correlation 
function has exponential decay in time, uniformly in space: For A small enough, 

with the decay rate grw as in (|4.551 . To verify l|5.23|l , we assume for concreteness that li — L and I2 — R, and we 
estimate by the triangle inequality 

\xr ~Xl\<IWr~x'l\ + 1 K^r + A) - {xr + + \x„ - xl\ (5.24) 

We note that the three terms on the RHS of ||5.24| | correspond (up to factors i) to the three spatial arguments 
multiplying 7 in the first line of l|4.55|l . By (|5.241 , at least one of these terms is larger than ^ — Hence we 
dominate | |4.55l l by replacing that particular term by i — | . Setting all other spatial arguments in | |4.55l l equal 
to zero, we obtain 

\\{^t).,,.^,.',,.'J < Ce'-™(^-^)*e-*l-«-^-l + C'e-(^'3'-»)' (5.25) 

Assuming that \x'j^ — Xj^\ > v*\t2 — ii| and using that rr^(7, A) = 0(7^)0(A^), see (|4.54l l, we choose 7 such that 
the first term of Il5.25b decays exponentially in t2 — ti with a rate or order 1. Hence, at high speed (> v*) Il5.23b is 
satisfied. At low speed (< v*), | |5.23t holds by the exponential decay of ip and the bound || Ai|| < Ce'^''*'^^*, which is 
easily derived from | |4.55l l. 

For li = I2, we can apply the same reasoning, and hence ( |5.23t is proven in general. However, in the case 
h = h, the proof is actually simpler. We can follow the same strategy as used for bounding (|5.20b , but replacing 
the propagation estimate (|2.12l l for Ut by the analogous estimate 114.61 for Ai . Indeed, the exponential decay in the 
case li = I2 was already present without the coupling to the reservoir, as explained in Section [5.4.11 whereas the 
decay in the case li 7^ ^2 is a nontrivial consequence of the decoherence induced by the reservoir. 



renormalized S — R correlation fct. 


\x\/t > V* 


\x\/t < V* 






li, I2 arbitrary 


exp. decay from 
decoherence of Af 


exp. decay from 

e±"'f^s(0,a;) 


exp. decay from 

tp{x,t) 



<g-tA^s.™^ for h,l2e{L,R} 



(5.23) 
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Along the same line, we note that the decay rate in Il5.23b cannot be made greater than O(A^), since the effect 
of the reservoir manifests itself only after a time 0(A^^). This should be contrasted with the decay rate for (|5.20b , 
which can be chosen to be independent of A. 



5.5 The renormalized model 

We have argued in the previous section that it makes sense to evaluate the perturbation expansion (|5.8l l in two steps 
by introducing a cutoff r for the temporal arguments of the correlation function (. The resulting cutoff reduced 
evolution was described in Section 15.31 By reordering the perturbation expansion, we are able to rewrite the 
reduced evolution Zt approximatively as 



where the restriction that t2 — ti > t reflects the fact that the short diagrams have been resummed. Note that 
it is somewhat misleading to call the remainder of the perturbation series 'higher order in A', since r will be A- 
dependent, too. 

The main tools in dealing with the renormalized model are 

1) The exponential decay of the renormalized joint correlation function, as outlined in Section |5l4l This prop- 
erty holds true thanks to the decoherence in the Markov semigroup Aj and the exponential decay for low 
('subluminal') speed of the bare reservoir correlation function. The latter is a consequence of the fact that the 
dispersion law of the reservoir modes is linear (see Lemma [5. 11 1. The necessity of the exponential decay of the 
renormalized joint correlation function for the final analysis will become apparent in Lemma [9.4l 

2) The integrability in time of the correlation function, uniformly in space, as stated in Lemma lS^ This property 
allows us to sum up all subleading diagrams in the renormalized model. This will be made more explicit in 
Section 1221 in particular in Lemma [9. 31 

The most convenient description of the renormalized model will be reached at the end of Section [8] and the be- 
ginning of SectionlH where a representation in the spirit of l|5.26b is discussed. The treatment of the renormalized 
model is contained in Section|9] 

5.6 Plan of the proofs 

In Section m we present the analysis of the cutoff reduced dynamics Z^ and the full reduced dynamics Zt, starting 
from bounds that are obtained in later sections. The main ingredient of this analysis is spectral perturbation theory, 
contained in Appendix |B1 

In Section[71 we introduce Feynman diagrams and we use them to derive convenient expressions for the cutoff 
reduced dynamics Z^ and the full reduced dynamics Zt. We will distinguish between long and short diagrams. 
The cutoff reduced dynamics contains only short diagrams. 

Section [8] contains the analysis of the short diagrams. In particular, we prove the bounds on Z^, which were 
used in Section[6l 

In Section |91 we deal with the long diagrams. In particular, we prove the bounds on Zt from Section |6l At the 
end of the paper, in Section [9341 we collect the most important constants and parameters of our analysis. A flow 
chart of the proofs is presented in Figure |6l 



<ti <t2 <t 

1*2 - til > T 




i:*{x2-xi,t2-ti) Zl_tI^.^j^Zl^_tI^^j^Zl^ +higher orders in A (5.26) 
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Spectral analysis 
M: Proposition |4!2| 



Control of short 
diagrams: Lemma 
EH 



Analysis of the 
cutoff model: 
Lemma 16.21 



Control of long 
diagrams: Lemma 
[63] 




Analysis of the 
cutoff model: 
Lemma 16.21 



Analysis of the full 
model: Theorem 



Spectral perturb- 
ation theory: 
Lemma lB.ll 



Properties of 
tjj{x,t): Lemmata 



Figure 6: A flow chart of the proofs. The arrows mean ' 
of the implication. 



Main results: The- 
orems 13.11 and 13.21 

implies". The arrows pointing to arrows specify the proof 



6 Large time analysis of the reduced evolution Zt and the cutoff reduced 
evolution 

In this section, we analyze the evolution operators Zt and Z[ starting from bounds on their Laplace transforms 

n^z) := J dte-''-Zl (6.1) 

and 

n{z) := j dte~*^Zt. (6.2) 

M+ 

These boimds are proven by diagrammatic expansions in Sections [71 [8| and [3 However, the present section is 
written in such a way that one can ignore these diagrammatic expansions and consider the boimds on Tl'^{z) and 
Ti{z) as an abstract starting point. Our results. Lemma 16.21 and Theorem 13.31 follow from these bounds by an 
application of the inverse Laplace transform and spectral perturbation theory. For convenience, these tools are 
collected in Lemma [B. II in Appendix [B] 



36 



6.1 Analysis of Zl 

Our main tool in the study of TV {z) is Lemma |6. II below. Loosely speaking, the important consequence of this 
lemma is the fact that we can represent the Laplace transform TV , defined in l|6.11 , as 

7^^(z) = (z - (-iad(r) + + A^{z)))-^, (6.3) 

where the operator {z) is "small" wrt. A^A^, in a sense specified by the theorem. Note that if we set A^{z) — 0, 
then the RHS of l|6.3l l is the Laplace transform of the Markov semigroup Af . This is consistent with the claim that 
Zl is 'close to' Af. 

The subscripts 'Id' and 'ex', introduced below, stand for "ladder" and "excitations", respectively. These sub- 
scripts will acquire an intuitive meaning in Section [7] when the diagrammatic representation of the expansion is 
introduced. The (sub)superscript t indicates the dependence on the cutoff t, but sometimes we will also use the 
(sub)superscript c. This will be done for quantities that are designed for the cut-off model but that do not necessar- 
ily change when r is varied. Lemma lOl can be stated for any t, but, as announced, it will be used for a A-dependent 

T. 

Lemma 6.1. For A small enough, there are operators TZl^{z) and in ^(^2(^))/ depending on A and r, satisfying 

the following properties: 

1) For Re z sufficiently large, the integral in (16. It converges absolutely in SS{SSi (^)) and 

TV{z) = (z - (-iad(i/s) + 7^^d(z) + TlU^z^K (6.4) 

2) The operators 7ll^{z),TZl^{z) are analytic in z in the domain Re z > Moreover, there is a positive constant Si > 
such that 

r \\J^'Rl^{z)J-4 = 0(A2)0(A2r), AV \ 0, A \ 

sup < (6.5) 

iimK^,„i<5i3c.>-^^ \ \\j^ni^{z)j^^\\ < X^C 

3) Recall the operator C{z), introduced in Section Wl] It satisfies 

r+oo 

sup IIj-, (7^[d(z)-A2£(z))^_,|| < A'C / dtsup\iP{xA)\ + X^tC (6.6) 

The proof of this lemma is given in Section|8l 

From Lemma lOl one can deduce, by spectral methods, that Z[ inherits some of the properties of the Markovian 
dynamics Aj. Instead of stating explicitly all possible results about Z[, we restrict our attention to Lemma [6. 2i in 
particular, to the bound | |6.7| |. This bound is the analogue of the bound | |4.55l l for the semigroup dynamics At, and 
it will be used heavily in the analysis of Zt in Section[8l 

Lemma 6.2. Let the cutoff reduced evolution Zl he as defined in Section \53[ with the cutoff time t = t(A) satisfying (|5.13t . 
Then there are positive numbers (5c > 0, Ac > and gc> such that, for < | A| < Ac and 7 < Sc, 

for constants c\^,c\ > 0, and with 

r,(7,A) = 0(A2)0(72) + o(A2) A\0,7\0 (6.8) 
where the bound o(A^) is uniform for 7 < 6c- 

The constants 5c > and decay rate gc > are in general smaller than the analogues Srw ai^d grw in the bound 
ll435l l. 
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Proof. We apply Lemma [BH in Appendix B with e := and 

V{t,e) UAZnpU^u (6.9) 

Aiiz,e) := C/,({7^L(z)}p + {7^^,(z)}p-A2i{ad(£(P))}J^7_. (6.10) 

N UAMjj^U-, (6.11) 

B := -[/,{ad(r)}p;7_, ~-ad(y) (6.12) 



and (p, I/) e 2)^;°"' with 



D;,"'^ := n <{ llmpl < mm{Si,5e), Im < i min(5i, 5^) (6.13) 



The set has been defined before Proposition 14. 2[ the bound on p, v involving 8\ ensures that we can convert 
the domain of analyticity in the variable n in Lemma fe.ll into a domain of analyticity in the variables (p, v), via the 
relation | |2.60|I . Similarly, the bound on p, involving 5^ ensures that 

sup ||C/,{ad(e(P))}p[/^,|| < C (6.14) 

as a consequence of the bound on Jf^z,A(e(P))J-K, provided by Assumption 12.11 and eq. | |2.12|| . 

We now check, step by step, the conditions of Lemma [B. II First, the continuity of V{t, e) and the bound l IB.ll 
follow from Lemma |231 and Statement 1) of Lemma |6A] The relation l|B.4|l is Statement 1) of Lemma [O] 
Condition 1) of Lemma FB. II is trivially satisfied since y is a Hermitian matrix on a finite-dimensional space. 

To check Condition 2) of Lemma [BAl we choose gA as gA — "^91°^ ^^i^ we will actually show that the bound l|B.61 , 
which is required in the region Re z > —X^gA, holds in the region Re z > —1/ (2t), as long as A^r is small enough. 
By Cauchy's formula, this implies that 



-Aiiz,X) = — 6 dz' \' ^OiX'r), ioTRcz> -X'gA, (6.15) 
Oz Ztti Jq^ yz — Z J 

with Cz a circle of radius 0(1/t) centered at z. Hence | |B.7| | follows. 

To check l lB.6t , we use the bound l|6.14l l for ad(e(P)). The boundedness of the other terms in Ai{z, A) follows 
immediately from l|6.5|l . 

Condition 3) contains conditions on the spectrum of Mp that are satisfied thanks to Proposition |4]2l It remains to 
check l|B.8|l . By the bound on 7^^^ (z) in l|6.5t , it suffices to check that, for any a e sp(ad(y )); 

l,(ad(r))J«(A2X - (-A2iad(e(P)) +7^[<J(-ia)))J7_«la(ad(r)) = oiX"), asA ^ (6.16) 

This follows by the estimate in (|6.6l l and the relation between A4 and C in I I4.41 . Note that we used that t(A) — > oo 
as A — >■ to get o(A^) from the estimate I l6.6b . 

Hence, we can apply Lemma [B. 1 1 and we obtain a number /r(p, A), a rank-one projector Pt{p, A) and a family 
of operators P^""" {t, p, A) such that 

{Z^p = e/-(f^^)*P,(p,A) +e-(^'3'°")*Pt''-(t,p,A) (6.17) 

for some g'"™ > (which can be chosen arbitrarily close to g^"^ by taking | A| small enough), and such that 

sup \\U,Pr{p, X)U-,\\ < C (6.18) 

sup sup\\U,Rl°^{p,t)U-,\\<C (6.19) 
Refr{p,X)>-X^g'r (6.20) 
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The above reasoning applies to small fibers, since we use the spectral analysis of Proposition 14.21 We now 
establish a simpler result about the cut-off reduced evolution {Z^}^, for large fibers. Let 

S);?*^'' n ||Imp| < mm{Si,S,),lmi^ < imin(5i,5,)| (6.21) 

Although for {p, v) G 'j;)high^ cannot apply Lemma [BAl we can still apply Lemma lB!2l to conclude that, for A 
small enough, the singularities of {7?.^(2^)}p in the domain, say. Re z > — 2A^(7^^^'' lie at a distance o(A^) from spA^p. 
One can then easily prove that {R'^{z)}^ is bounded-analytic in a domain of the form Re z > —}?grw^ + o(A^) and 
hence 

{z;}p===i?f«^p,t)e-(^^«?'''')* (6.22) 
with a rate 5^*^'' > (which can be chosen arbitrarily close to g^w^'' by making A small enough) and 

sup sup||i7^i?^'s''(t,p,A)C/_^|| < C (6.23) 

Finally, we note that one can easily find a constant 5c such that S)^""" and '^^'^sh gj.g Qf ^j-^g fQj-jn | |4.3Hl and (|4.321 
with the parameters 5c instead of 5rw (the parameter p*^, does not need to be readjusted). 

With the information on obtained above, we are able to prove the bound l l6.7t by the same reasoning as we 
employed in the lines following Proposition 32] to derive the bound l|4.55t . 

The function r^- (7, A) in the statement of Lemma [6l2l is determined as 

r,(7,A):= sup Re/,(p,A) (6.24) 

peT<',|p|<7 

and the bound I I6.8I I follows by | |4.54| | and 

A(p,A)-A2/(p)-o(A2), A\0, (6.25) 

which follows from (|B.13|I in Lemma [BA] The decay rate gc is chosen as.gc mm{gl°^ , g'^^^^) . This concludes the 
proof of Lemma |621 

□ 

We close this section with two remarks which are however not necessary for an understanding of the further 
stages of the proofs. 

Remark 6.3. As apparent from the bound 1 16.61 , one cannot take r = const, since in that case, this bound becomes O(A^) 
instead of o{\^). This would mean that there is a difference ofO{\^) between Zl and kt, whereas the important terms in Aj 
are themselves ofO()?). This is however not an essential point: as one can see from the classification of diagrams in Section^ 
one could easily modify the definition of the cutoff model such that Zl is close to At even at t = const. This can be achieved 
by performing the cutoff on the non-ladder diagrams only, which is a notion introduced in Section^ The true reason why t 
must diverge to 00 when A \ will become clear in the proof of Lemma W^ in Section \93\ 

Remark 6.4. One is tempted to say that any claim that is made about Zt in Section\3\could be stated for Z[ as well. While 
this is correct for Theorem [l3[ it fails for Proposition \3.5\ The reason is that the identity f{p = 0, A) = follows from the fact 
that Zt conserves the trace of density matrices, as it is the reduced dynamics of a unitary evolution. This is not true for Z[, 
and hence we cannot prove (or even expect) that /^(p = 0, A) =0. 

6.2 Spectral analysis of Zt 

In this section, we state Lemma 16.51 the r = 00 analogue of Lemma 16.11 This Lemma leads to our main result, 
Theorem l3.31 via reasoning that is almost identical to the one that led from Lemma [6. II to Lemma [6.21 

Essentially (and analogously to Lemma |6A] |, Lemma r6.5l states that the Laplace transform Tl{z), defined in l|6.2b , 
is of the form 

n{z) = {z- (-iad(r) + X^M + A{z)))-^ (6.26) 

where A{z) is 'small' w.r.t. X^M. 
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Lemma 6.5. There is an operator TZc^{z) £ SS{^^2 i'^s))/ depending on A and satisfying the following properties, for A small 
enough: 

1) For Re z sufficiently large, the integral in (|6.2| | converges absolutely in ^(^2(^)) and 

n{z) = (7^^(z)-l - 7^e.(z))-l (6.27) 
where TV{z) was introduced in \6.1\ and t = t(A) was defined in (|5.131 . 

2) There are positive constants Sex, gex such that the operator TZex{z) is analytic in z in the domain Re z > —X'^gex and 

sup ||J«7^ex(^)^-«|| = o(A2), A\0. (6.28) 

The proof of Lemma 16.51 is contained in Section |8l Starting from Lemma 16.51 we can prove our main result, 
Theorem l3.3i by the spectral analysis outlined in Appendix |Bl 

Proof of Theorem 1331 We apply Lemma lB.lj with e := A^ and 

V{t,e) := UAZt}pU-, (6.29) 

M^, e) := ({7^e.(z)}, + + {7^^,(^)}p - A^I {ad(£(P))} J C/_. (6.30) 

N := U,{M}pU-, (6.31) 

B := -t/,{ad(y)}^C/_, ~-ad(y) (6.32) 

and 

(p,!/) e Dj,™ n ||Imp| < min((5e.,<5e),Imi/ < i min(,Se., <5e)| (6.33) 

Hence, the only difference with the relations (|6.9H6.10l6.11H6.12|l is that we have added the term {'Tiex{z)}p in 
| |6.30|I , we consider Zt instead of in l|6.29b , and we replace 5i by 6ex in (|6.33l l. This means that we can copy the 
proof of Lemma [6.21 except that, in addition, we have to check the bounds ||B.61 and (|B.7|l for the term TZf.x{z)- We 
choose gA '■= ^/2gex- Then the bound 1IB.6I 1 follows from ll6.28l l, and l|B.7) l follows since, by the Cauchy integral 
formula and (|6.281 , 

sup II— f/,{7^e.(z)}^C/_,|| = |Rcz-(-A2g,,)f'o(A2) ^o(|A|"), A\0 (6.34) 

where we use the same argument as in (|6.15l l, but with a circle radius equal to ^ |Re z — {—X'^gex)\ ■ This application 
of the Cauchy integral formula is the reason for the factor i into the definition of gA- Lemma iB.ll yields the 
function f(p, A), the rank-one projector P{p, A) and the operator R''°^{t,p, A) required in the small fiber statements 
of Theorem l3.3l 
For 

{p,i^) e Dj?;^'' n ||Imp| < miniSex, Se),lmiy < i min(5e., , (6.35) 

we can again apply Lemma [B. 21 to derive the large fiber statements of Theorem l3.3l As in the proof of Lemma [6. 2i 
we can again choose parameters 5,p* such that domains Ji'o™^ as defined in (|3.15H3.16l l, are included in the 

domains for {p, u) specified by l|6.33|l and (|6.35l l. □ 



7 Feynman Diagrams 

In this section, we introduce the expansion of the reduced evolution Zt and the cutoff reduced evolution Z^ in 
amplitudes labelled by Feynman diagrams. These expansions will be the main tool in the proofs of Lemmata 16. II 
and 16.51 We start by introducing a notation for the Dyson expansion of Zt which is more convenient than that of 
Section |5l2l 
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7.1 Diagrams a 

Consider a pair of elements in / x Z"^' x {L, R] with / c IR.+ a closed interval whose elements should be thought 
of as times. The smaller time of the pair is called u and the larger time is called v, and we require that u ^ v, i.e. 
u < V. The set of pairs satisfying this constraint is called S}. We define S" as the set of n pairs of elements in 
I X y. {L, R] such that no two times coincide. That is, each a E Yj"^ consists of n pairs whose time-coordinates 
are parametrized by (u^, vi), for i = 1, . . . ,n, and with the convention that Ui < Vi and Ui < u^+i. The elements a 
are called diagrams. As announced in Section I5l2l there is a one-to-one mapping between, on the one hand, a set of 
triples {ti{a), Xi{<T), ^i(cr))i"i with ti < ti^i and ti € /, together with a pairing tt E Vn, and, on the other hand, a 
diagram cr G S" as defined above. 

To construct this mapping, proceed as follows: Choose from the pairing tt the pair (r, s) for which r = 1 and set 
ui = tr,vi = tg. The pair {{tr,Xr,lr), {ts,Xs,ls)) becomes the first pair in the diagram a. Then choose the pair 
(r', s') E TT such that 

r' = min{{l,2,...,2n}\{r,s}} (7.1) 

Set U2 = tr', V2 — tgi. The pair ((ir', Xr', Ir'), {ts', Xs', Is')) becomes the second pair of cr. Repeat this until one has n 
pairs, each time picking the pair whose r is the smallest of the remaining integers. The mapping is easily visualized 
in a picture, see Fig. [Zl 

We also use the notation t{(T),x{(j),l{a) to denote the 'coordinates' of the diagram a. Here, t{a),x{a),l{(T) are 
2n-tuples of elements in /, Z''. {L, R}, respectively, and such that the i'th component of these 2n-tuples constitutes 
the z'th triple {U{a), Xi{(T),kia)). 

Note that the time-coordinates t = ti{a), . . . , t2n{(j) can also be defined as the ordered set of times containing 
the elements {ui, Vi,i — 1, . . . , n}. Evidently, the triples (ti(cr), Xi{(T), h{cr))f"i do not fix a diagram uniquely since 
the combinatorial structure that is encoded in tt is missing. That combinatorial structure is now encoded in the 
way the time coordinates t{a) are partitioned into pairs {ui, Vi), see also Figure[7l 

We drop the superscript n to denote the union over all n > 1, i.e. 

S/ := U S/ (7.2) 

n>l 

and we write \a\ — n to denote that ct G Sp. 
We define the domain of a diagram as 

n 

DomCT := y [u„ v,], for cr G (7.3) 

1=1 

We call a diagram cr G S/ irreducible (notation: ir) whenever its domain is a connected set (hence an interval). 
In other words, cr is irreducible whenever there are no two (sub)diagrams cri , cr2 G S/ such that 

cr = cri U 0-2 , and Domcri n Domcr2 = (7-4) 

where the union refers to a union of pairs of elements in / x Z'' x {L, R}. For any cr G S/ that is not irreducible, we 
can thus find a unique (up to the order) sequence of (sub) diagrams ai, . . . such that 

cri , . . . , cr^ are irreducible and cr = cti U . . . U cr^ (7.5) 

We fix the order of CTi, . . . , a^n by requiring that maxDomcri < miiiDomcri+i and we call the sequence (crj, . . . , cr„j) 
obtained in this way the decomposition of a into irreducible components. 

We let E" (ir) c S" stand for the set of irreducible diagrams cr (with n pairs) that satisfy Domcr = /, that is, 
ui = <i(cr) ~ mini and max; ti(cr) = maxi Vi — max/. 

A diagram cr G S/ (ir) is called minimally irreducible in the interval / whenever it has the following property: 
For any subdiagram a' c a, the diagram a\a' does not belong to S/(ir). Intuitively, this means that either the 
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diagram g' contains a boundary point of / as one of its time-coordinates, or the diagram cr \ cr' is not irreducible. 
The set of minimally irreducible diagrams (with n pairs) is denoted by S" (mir). See pictures[7]and|8]for a graphical 
representation of the diagrams. Since, up to now, most definitions depend solely on the time-coordinates, we only 
indicate the time-coordinates in the pictures. In the terminology introduced below, we draw equivalence classes of 
diagrams [cr] rather than the diagrams cr themselves. 




Ul Vi U2 Us U4 V2 Va U5 V3 V5 

Figure 7: A diagram (T g S/ with \a\ = 5. Its time coordinates are shown explicitly. Note that the parametrization by Ui, Vi encodes the 
combinatorial structure (the way the times are connected by pairings), whereas the ti are ordered. We consistenly draw the long pairings (see 
later) above the time-axis and the short ones below. 




Figure 8: The left figure shows an irreducible diagram cr in the interval / = [/_ , /+] with \a\ =7. This diagram is not minimally irreducible. 
The right figure shows a minimally irreducible subdiagram. In this case, there is only one such minimally irreducible subdiagram, but this 
need not always be the case. 



A diagram cr in S/ for which each pair of time coordinates (w, v) satisfies jw — w| > r, or |u — u| < r, is called 
long, or short, respectively. The set of all long/small diagrams with n pairs is denoted by S"(> r) / (< r). Note 
that S"(> r) U I]y(< r) is strictly smaller than E" whenever n > I. 

In addition to the sets I]"(ir), S"(mir), S"(> t), we will sometimes use more than one specification (adj) to 
denote a subset of E/ or S", and we will drop the superscript n to denote the union over all |cr|, as in (|7.3b , for 
example, 

I]?(<T,ir), I],(>r,mir) (7.6) 

are the sets of short irreducible diagrams with |cr| = n and long minimally irreducible diagrams, respectively. 
On the set Ey, we define the "Lebesgue measure" dcr by 

j AaF{a) := j dwi . . . du„ j dwi . . . dw„ ^ i^(cr) (7.7) 

/_<ui<. ..<«„</+ M,<«i a;((j),i(o-) 

where I — [/_,/+]. 

Since S" (ir) is a zero-measure subset of E", the definition of the measure da on E" (ir) has to modified in an 
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obvious way: For all continuous (in the time coordinates t{a)) functions F on , we set 

j daF{a) = j d<7S{maxt{a) - I+)6{mmt{<7) - I^)F{a) 



(7.8) 



S?(ir) 



where the Dirac distributions (5(/+ — •) and S{I^ — •) are a priori ambiguous since I-,I+ are the boundary points 
of the interval /. They are defined as 



Si- -!+):= lim S{- - s), d{- - I-) := lim S{- - s) 



We extend the definition of the measure do- also to S/ and the various I]/(adj) by setting 

J daF{a):=Y, J d(7F„(a) 



(7.9) 



(7.10) 



Sj(adj) 



- Sy(adj) 



where F„ is the restriction to S"(adj) of a function F on S/(adj). 

We will often encounter functions of <t that are independent of the coordinates x{a),l{<T) and that must be 
integrated only over t{(T) and summed over \(t\. To deal elegantly with such situations, we let [a] stand for an 
equivalence class of diagrams that is obtained by dropping the x, ^-coordinates. That is 



\CT\ = \a'\ 

u.,{(j) = Ui{a'),Vi(a) = Vi{a'), for alii = 1, . . . , |cr| 



(7.11) 



The set of such equivalence classes is denoted by IItE/ (the symbol can be thought of as a projection onto the 
time coordinates) and we naturally extend the definition to IItS/ (adj) where adj can again stand for ir, mir, > r, < r 
. The integration over equivalence classes of diagrams is defined as above in (|7.7b and (|7.8l l, but with 

omitted, i.e., such that for all functions F on S/(adj): 



d[<j]F{[a] 



j dCTF(cr), with F{[a]) = ^ F[a) 



(7.12) 



nTS^adj) 



S/(adj) 



Lemma [7.11 contains the main application of this construction. It is in fact a simple — L°°-bound. 
Lemma 7.1. Let F and G he positive, continuous functions on S/. Then 



j daF{a)G{a) < J d[a] 



S/(adj) 



nTSi(adj) 



sup G{cr) 



sup 



(7.13) 



where it is understood that the sum and sup over x{a),l{(j) are performed while keeping \(t\ and t{a) fixed. 



In (|7.131 , the sum/sup, over x{a),l(a) is in fact a shorthand notation for the sum/sup over all a' such that 
[cr'] = [a] for a given a. Hence, supj,(o-)^;(o-) G{a) is a function of [a] only, as required. The supremum sup((o.) is over 
/2kl^ with \a\ fixed. Hence, the second factor on the RHS of ll7.13l l is in fact a function of\a\ only. 

Proof. We start from the explicit expressions in i7.7i or (|7.8b , and we use a — L°° inequality: first for the sum 
over x{a),l{(j) and then for the integration over Ui,Vi. □ 
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7.1.1 Representation of the reduced evolution Zt 

Recall the operators Vt((ti, defined in l|5.9|l . Since, by the above discussion, there is a one-to-one corres- 

pondence between a diagram a and (tt, (i^, Xi, ) where tt e Pio-i and ti < U+i, we can write Vt(cr) instead of 
VtiiU,x„k)f^l) andC(CT) instead ofC{Tr,{U,x^,k)f^l),i.e. 

Vtia) Vt((t.(a),x,(a),Z,(a))t';l) (7.14) 

and 

{ijj{x' — x,v — u) I = I' = L 

f{x'-x,v-u) l^l'^R (7-^5) 

i^{x' -x,v-ti) l = L,l' = R 
ijix'-x,v~u) l = R,l' = L 

As a slight generalization of the operators Vt(o'), we also define V/(o') for a closed interval / [/_ , J+J by 

V/(cr) ■.= Ui^-t2„^x2„,i2„ ■ ■ ■^X2,i2^t2-tilxuh^ti-i-, for (T such that DoniCT C / (7.16) 

The only difference with Vt (cr) is in the time-arguments of Ut at the beginning and the end of the expression. 
With this new notation, Vt(cr) = V[o,t](cr). Next, we state the representation of the reduced evolution Zt as an 
integral over diagrams 

Zt=llt+ f df7C((T)V[o,t]((T) (7.17) 

Similarly, the cutoff dynamics Z^ is represented as 

Zl=Ut+ j d(TC(fT)V[o,t](a) (7.18) 

S[o,t](<T) 

Formulas l|7.17l l and | |7.18| | are immediate consequences of l |5.8t and I I5.15I I, respectively 

We use the notion of irreducible diagrams a to decompose the operators V[o tj (cr) into products and to derive a 

new representation, | |7.23| |, for Zt and Z^ 

Let (fJi, . . . ,(Tp) be the decomposition of a diagram cr G Sjo tj into irreducible components. Define the times 

si, . . . , S2p to be the boundaries of the domains of the irreducible components, i.e., [s2i-i, S2i\ — DomcTi, for i = 

1, . . . Then 

V/(cr) = Ul^-S2p "^^[s2p-l,S2p]('^p)^S2p-l-S2p-2 ■ • -^83-82 , s^] (fJl ) Z^^i , (7.19) 

as can be checked from l l7.15H7.16l l. Here, the essential observation is that all time coordinates of (T; are smaller than 
those of (Ti+i. We introduce 

2't:= I d<(a)V[o,ti((T), Z^'":= J d<(a)V[o,t] (a) (7.20) 

S[o,t](ii-) S(0.t](<-r,ir) 

and we remark that the definitions in | |7.20| | allow for a shift of time on the RHS, that is 

J daC{a)Vi{a), for any I ^ [s, s + 1], seR (7.21) 

and similarly for Z^'". By this time-translation invariance, the factorization property 1 17.191 and the factorization 
property of the correlation function in l |7.15|l , i.e., 

p 

C(aiU...Uap) = nC(^.), (7.22) 
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we can rewrite the expression 117.1 7|l as 



Zt^ j . . . ds„, (Ut^s„A-s^-^ ■ ■ ■Us,-s.Z%_,^Us,) . (7.23) 

™62Z+ o<si<...<s„<t 

where the term on the RHS corresponding to m = is understood to be equal to Ut- The idea behind (17.231 1 is that, 
instead of summing over all diagrams, we sum over all sequences of irreducible diagrams. An analogous formula 
holds with Zt and Z" replaced by Z[ and Z^'". 

7.2 Ladder diagrams and excitations 

We are ready to identify the operators TZl.^ {z) and 7?.f^(z), whose existence was postulated in Lemma [6lT] and the 
operator TZex [z), which was postulated in Lemma [631 

The Laplace transform, 'R.{z),oi Zt has been introduced in l|6.2l l. We calculate 'R{z) starting from l|7.23|l 

7^(z) = / die-*^Zt (7.24) 



(z + iad(i/s))-i (7.25) 



= ^ (z + iad(i/s))-i [ dte-'^Z, 

m>0 L 

= (z + iad{Hs) - niriz))-\ with7^ir(z) := / die~*^Zf (7.26) 

The second equality follows by J^^ Ate~*^^Ut — (2;+iad(iJs))^^ for Re z > 0. The third equality follows by surmning 
the geometric series. 

An identical computation yields 

TViz) = {z + iad(iJs) - ni{z))~^, with ^[^{z) / dt e-*^Z['''' (7.27) 

The definition of TZl^{z) and TZl^{z) relies on the following splitting of Tllj.{z) 



nUz) := j Ate~'^ j C(^)lk|=iV[o,t](cT) (7.28) 

K+ S[o.tj(<r,ir) 

-KA-A ■■= j die-*^ j CMl|.|>2V[o,t](a) (7.29) 

R+ S[o.tj(<r,ir) 

The subscripts refer to "ladder'- and 'excitation'-diagrams. The name 'ladder' originates from the graphical repres- 
entation of diagrams whose irreducible components consist of one pair (it is standard in condensed matter theory). 
Since obviously TZl^{z) + Tll^{z) = 7?.f,.(z), the relation (|7.27l l implies Statement (1) of Lemma [6lT] 

In the model without cutoff, we do not disentangle ladder and excitation diagrams, since every diagram that 
contains a long pairing, is considered an excitation. We can thus define 

n,x{z):=nUz)~TZ[,{z) (7.30) 

We will come up with a more constructive representation of TZex (z) in formula (17.331 1. 

7.2.1 The reduced evolution as a a double integral over long and short diagrams 

We develop a new representation of Z" by fixing the long diagrams, i.e., those in I][o,t](> t), and integrating the 
short ones. 
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We define the conditional cutojf dynamics, Ct{cri), depending on a long diagram ai e ^[o,t] (> t), as follows: 



Ct(cr;) = l(T,eS[o,tj(>Tar)V[o,t](cri) + J dcrC(cr)V[o,t](o-U cr;) 



(7.31) 



S[0,t)(< r) 
(TiUcT & E[o,t]{ir) 



In words, Ct{ai) contains contributions of short diagrams a G I]i(< r) such that ai U a is irreducible in the interval 
[0,t]. Hence, if ai is itself irreducible in the interval [0,t], then there is a term without any short diagrams; this is 
the first term in (|7.311 . In general, ai need not be irreducible. Note that the constraint on a (in the domain of the 
integral) in the second term of | |7.31| | depends crucially on the nature of cr/. In particular, if Domcr; does not contain 
the boundary points or t, then a has to contain or t, and this introduces one or two delta functions into the 
constraint on cr. To relate Ct((T/) to Z", we must explicitly add those cr; that contain one or both of the times and 
t. This is visible in the following formula, which follows from (|7.31t and the definition of Z" in | |7.20|I . 

Zl' - = I daiCiaOCtiai) [1 + S{h{ai)] [l + 5(t2|.,|(^0 - *)] (7-32) 

S[o,t](>r) 

We must subtract Z/^'" on the LHS, since all contributions to the RHS involve at least one long diagram. The 
5-fimctions on the RHS are defined as in l|7.9b . 

The following formula is an obvious consequence of | |7.30l l and | |7.32| |: 



TZeAz) = J dte''' J daiaai)Ct{cTi) [1 + Sihiai)] [1 + Sit2\,^\{ai) ~ t)] (7.33) 

M+ E[o,t](>r) 

All of Section |8] will be devoted to proving good bounds on TZex (z), as claimed in Lemma [6.51 



7.3 Decomposition of the conditional cutoff dynamics Ct{cri) 

Our next step is to decompose the conditional cutoff dynamics Ct{(Ti), as defined in | |7.31| l, into components. Since 
Ct{(Ji) is defined as an integral over short diagrams cr, we can achieve this by classifying the short diagrams a that 
contribute to this integral. The idea is to look at the irreducible components of cr whose domain contains one or 
more of the time-coordinates of cr; (In our final formula, ll7.45l l, these domains correspond to the intervals [s\,, s\]). 
The irreducible components whose domain does not contain any of the time coordinates of cr; can be resummed 
right away, and they do not play a role in our classification (this corresponds to the operators Z^ in | |7.45| |). We 
outline the abstract decomposition procedure in Section [7.3.11 and we present an example (with figures) in Section 
[73:21 



7.3.1 Vertices and vertex partitions 

Consider a long diagram cr/ e S[o,(](> t) with |cr/| = n and time-coordinates t(cr/) = (ti, . . . , i2n)- With this 
diagram, we will associate different vertex partitions Z. First, we define vertices. A vertex [ is determined by a label, 
hare or dressed, and a vertex set S{\), given by 

S[\.) — {tj,tj+i . . . , tj+„i_i}, for some l<j<j + ni — l<2n (7.34) 

Hence, the vertex set is a a subset of the times {ii(cr;), . . . , t2n{<^i)}- Moreover, a vertex I with \S{1)\ > 1 is always 
dressed. Hence, if [ is bare then S{1) is necessarily a singleton, i.e., S{1) — {tj} for some j. 

A vertex partition £ compatible with cr/ (Notation: £ ^ cr/) is a collection of vertices Ii , . . . , [,„ such that 

• The vertex sets S{li), . . . , S{lp) form a partition of {ti(cr/), . . . , i2n(o'/)}. By convention, we always number the 
vertices in a vertex partition such that the elements of S'(Ia;) are smaller than those of S{lk+i)- The number, p, 
of vertices in a vertex partition is called the cardinality of the vertex partition and is denoted by |£|. 
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• Any two consecutive times tj,tj+i such that [tj,tj+i] ^ Domcr/, belong to the vertex set S{lk) of one of the 
vertices Ik- Such a vertex Ik is necessarily dressed since its vertex set contains at least two elements. 

• Ifti = 0, then S{li) = {ti} and [i is bare. If ti > 0, then S{li) 3 ti and li is dressed. 

• If hn — t, then 5([,„) — {t2,i} and [,„ is bare. If t2n < t, then S{lm) 3 t2n and [,„ is dressed. 
The idea is to split 

Ct(aO= ^ Ct(a,,£) (7.35) 

where the sum is over all £ compatible with cri and Ct [ai , £) contains the contributions of all short pairings a that 
match the vertex partition £; 

{V dressed Ik : 3!irr. component (Tj C a such that S{lk) C DomCTj 
and S ( Ik' ) n DoniCTj = for all k' ^ k (7 36) 

Vbare Ik : S'(Ife) n Domcr = 

For the sake of completeness, we define the operators Ct{<7i, £), below, but, in Section [7.3.31 we will provide 
a more constructive expression for them. First, assume that the vertex partition £ contains at least one dressed 
vertex. Then Ct {cri , £) is defined by restricting the second integral in (|7.311 to those a that match £; 

Ct((T;,£):= j dcrC((T)V[o,t] (crUcr;):^ matches £ (7-37) 

S[0,t](< -r) 
o-i U(T G S[o,t](ir) 

Next, we assume that the vertex partition £ contains only bare vertices. If <Ji is irreducible in the interval [0, t], i.e., 
(T; e S[o_t] (> T, ir), then the vertex partition with only bare vertices is compatible with cr;. If cr; ^ ^[o.i] (> ir)/ then 
this vertex partition is not compatible with cr; . Hence, we assume that cr/ G '^[o.t\ (> ir) and we define 

Ct(fT/,£) :=V[o.t](^/)+ j daC(fT)V[o,t](aUaO X (7.38) 

S[o.tl«r) 

The second term is the same as in | |7.37|| (but specialized to the partition with only bare vertices) and the first term 
is a contribution without any short diagrams. This first term equals the first term (on the RHS) of l|7.31t . 

In Section 17.3. 2i we give examples of vertex partitions that are intended to render the above concepts more 
intuitive. 



7.3.2 Examples of vertex partitions 

We choose a long diagram ai e (> r) which consists of three pairs such that the time coordinates [ui, Vi)'-^^ 
are ordered as 

tl t2 ts ti t5 te ,„ OQ-, 

^ < U\ < U2 < V\ < U'i < V2 < V-i — t 

Hence, cr; is irreducible in the interval [ui, t], but not in the interval [0, t], at least not if <i ^ 0. 

Below we display three diagrams a S I][o.t](< t) satisfying the condition gi\J a ^ S[o,t](ir). To assign to each 
of those diagrams a vertex partition, we proceed as follows. Starting on the left, we look at the time-coordinates 
i(cr;) and we check whether these times are 'bridged' by a short pairing, i.e., whether they belong to the domain 
of a short diagram. If this is the case then such a time belongs to the vertex set of a dressed vertex. The vertex set 
of this vertex is the set of all time-coordinates that are connected to this point by short pairings. If this is not the 
case, i.e., if a time-coordinate of cr; is not 'bridged' by any short pairing, than such a point constitutes a bare vertex, 
whose vertex set is just this one point. 
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t2 ta ^4 *5 *6 

Figure 9: A long diagram with time coordinates as in )7.39) 



Actually, for the first time-coordinate (in our case ui), this is particularly simple. Either the first time-coordinate 
is not equal to 0, in which case it has to be 'connected' by short pairings to (indeed, if this were not the case, then 
CT; U (T cannot be irreducible in [0, t]), or the first time-coordinate is equal to 0, in which case it cannot be connected 
by short pairings to the second coordinate, because then the first time-coordinate of the short diagram would have 
to be as well, which is a zero measure event (for this reason, we have excluded this case in the definition of the 
diagrams in Section (TAJ. In our example ui 7^ 0, and one checks that, in all three choices of a, there are short 
diagrams connecting ui and 0. 

Let us determine the vertices in the three displayed figures 



vertex partition 1 



vertex partition 2 



vertex partition 3 



h 


{tl} 


dressed 








[2 


{t2} 


bare 






{h} 


dressed 


[3 


{h} 


bare 




[2 


{t2} 


bare 


k 


{U} 


dressed 




h 


{h, ^5} 


dressed 


[5 




bare 




k 


{^6} 


bare 


[6 




bare 









(1 


{tl} 


dressed 


(2 


{t2M 


dressed 


I3 


{t,} 


bare 


(4 


{is} 


bare 


(5 


{is} 


bare 



In the example displayed above, it is also very easy to determine which vertex partitions £ are compatible with 
ai (£ ~ ai). Apart from the fact that the vertex sets S{\.k) of the vertices in £ have to form a partition of {ti, . . . , to], 
we need that Ii is dressed and [|£| (the last vertex in the partition) is bare. 

To each vertex ii^ in the above examples, we can associate time coordinates s\ and s| as the boundary times 
of the domains of irreducible diagrams bridging the times in the vertex. Eventually, we intend to fix a vertex 
partition and associated time coordinates s* and and to integrate over all short diagrams that are irreducible in 
the interval [s*, s^]. This integration gives rise to the vertex operators, see eq. UA2\ . To illustrate this, we zoom in 
on a part of a long diagram, shown in Figure [TTJ A formal definition is given in the next section. 



7.3.3 Abstract definition of the vertex operator 

Let £ be a vertex partition compatible with cr/, with vertices [k,k = 1. . . . , |£|. In what follows, we focus on one 
particular vertex \k which we assume first to be dressed. The vertex is assumed to have a vertex set S{\k) = 
{ij, . . . , tj+m-i}- This means in particular that the time-coordinate tj_i belongs to the vertex set of the vertex 
[fc-i (unless i = 1) and the time-coordinate i^+m belongs to the vertex set of the vertex [^+1 (unless j+m— 1 = 2|o-; |). 
We fix an initial time s\ and final time s\. such that 

tj^l < S\ < tj < tj+,n-l < 4 < tj+m (7.40) 

where it is understood that = if j = 1 and tj+m = tifj + m- l = 2\ai\. The vertex operator B{lk, s]^, s^) is 
defined by summing the contributions of all a G Ej^,^ j (< r, ir) 



48 




Il I2 [3 [4 (5 

Figure 10: The picture shows three different choices of short diagrams a e '^{f3,t\ (< ''")• Recall that short diagrams are drawn below the 
horizontal (time) axis. In each picture, we show the resulting vertex partition by listing the vertices Ii, b, • ■ ■• The dressed vertices are denoted 
by a horizontal bar whose endpoints represent the vertex time-coordinates s\,s\. The bare vertices are denoted by a short vertical line whose 
position represents the (dummy) vertex time coordinates s^, = = tj ■ The time-coordinates of the bare vertices are not shown since they 
coincide with time-coordinates of long pairings. For example, in the bottom picture, (1 , 12 are dressed and I3, (4, (5 are bare. 



To write a formula for the vertex operator B{ik-, Sfc, we need to relabel the time-coordinates of ai and a e 

Consider the m triples {U{ai), Xi{ai), h{ai)), for i = j, . . . ,j + m - 1, i.e., a subset of the 2\ai \ triples determined 
by the long diagram 0-;, and the 2\a\ triples U{a),Xt{a), li{a) with i — I, . . . ,2\a\ determined by a G Ej^,^ ,5t^](< r, ir). 

We now define the triples {t", x'-, by time-ordering (i.e. such that t'- < t"_^_^) of the imion of triples 

{U{a),Xi{a),k{a))f:l and {ti{ai),Xi{ai),k{ai))itT-' ■ (7.41) 
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B{l,s\s^)= J da — 



Short diagrams a 



Figure 11: A part of a long diagram ct; G ^[o.t] (> ''") is shown, suggesting a dressed vertex I with vertex set S(() = {tj, . . . , <j+4}. The end 
points of the pairings that are 'floating' in the air are immaterial to this vertex, as long as they land on the time-axis outside the interval [s', s']. 
The vertex operator B(l, s', s') is obtained by integrating all diagrams in Sj^i ^/j (< r, ir). 



The vertex operator B{lk, s^) is then defined as follows 

B{lk,sl.si) J daC(<T)V[,. (7-42) 

S, i ,t,(<r,ir) 

where the dependence of the integrand on a is imphcit in the above definition of the triples (t", x'^ , I'/). The double 
primes in the coordinates {t^, x'l . I'D are supposed to render the comparison with later formulas easier. 

We now treat the simple case in which the vertex \k is bare. In that case, there is a j such that S{\.k) — {tj} and 
the vertex operator is simply defined as 

4 = 4=*. (7.43) 
Hence, in this case, the vertex time-coordinates s)., s\. are dummy coordinates, see also Figure [TOl 

7.3.4 The operator C* (ct/ , £) as an integral over time-coordinates of vertex operators 

We are ready to give a constructive formula for Ct((T/, £), as annoimced in Section r7.3.1l First, we define the integ- 
ration measure over the vertex time-coordinates s\, s\.; 

Vivi Jl ds\<ls{ 

fc = 
Ifc dressed 

To understand this formula, we observe that only non-dummy vertex time coordinates need to be integrated 
over. A dummy vertex time coordinate is a time coordinate whose value is a-priori fixed by ai and the vertex par- 
tition £. The non-dummy times are the time coordinates of the dressed vertices, except at the temporal boundaries 
0, t, where such a time coordinate is also a dummy coordinate. The terms between { j-brackets in formula (|7.44)l 
take care of this . Finally, the formula for Ct [cri , £) is 

Ct{ai,£)^ f VsW BH\^\,s\s,^,s\^^)z:. ...z;; ,^([2,4,4)2;. ,;S(li,4,4) (7.45) 

< sj, < 4 < t 
s\, < sJ,, for k' > k 

where the indices k, k' correspond to vertices 1^, Ik'- only the time-coordinates of dressed vertices are integrated 
over, even though all vertices appear on the RHS. This formula can be checked from the definition l |7.37t and the 
explicit expressions for the vertex operators B{-;-,-) above. The cutoff reduced dynamics in | |7.45l l appears by 
summing the small diagrams between the vertices, using formula (|7.181 . 



5{s[) li dressed 1 J - 1) t|£| dressed 

1 [i bare J 1 1 li^i bare 



(7.44) 
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8 The sum over ''small" diagrams 



In this section, we establish two results. First, we analyze the cutoff-dynamics Z^. . The main bound is stated in 
Lemma 18.31 and a proof of Lemma 16.11 (concerning the Laplace transform of Zl^ is outlined immediately after 
Lemma [8.31 Second, we resum the small subdiagrams within a general irreducible diagram: Recall that the con- 
ditional cutoff dynamics Ctipi) is defined as the sum over all irreducible diagrams in [0,t] containing the long 
diagram crj. In Lemma 18.61 we obtain a description of Ctipi) that does not involve any small diagrams. In this 
sense, we have performed a blocking procedure, getting rid of information on time-scales smaller than r. 

Since this section uses parameters and constants that were introduced earlier in the paper, we encourage the 
reader to consult the overview tables in Section|9]4l 



8.1 Generic constants 

In Sections H] and m we will state bounds that will depend in a crucial way on the parameters A, 7 and r. The 
parameter 7 is a momentum-like variable used to bound matrix elements in position representation, see below in 
Section l8l3l It appeared first in Section l4.3.11 To simplify the presentation, we introduce the following notation and 
conventions. 

• We write 0(7) for functions of 7 > with the property that 0(7) is decreasing as 7 00, and 0(7) is finite, 
except, possibly, at 7 = 0. It is understood that 0(7) is independent of A. 

• We write 0(7, A) for functions of 7 > and A G K that have the asymptotics 

c(7,A)=o(7")0(A2)+o(A2), 7^0, A^O (8.1) 

• We write c'(7, A) for functions of 7 > and A e M that have the asymptotics 

c'(7,A)=o(7°)0(A2) + c(7)o(A2), 7 ^ 0, A ^ (8.2) 

• The cutoff time r = t(A) is treated as an implicit function of A, satisfying l |5.13l l. In particular, 0(7, A) and 
c'(7, A) can depend on t. 



8.2 Bounds in the sense of matrix elements 

In Section 1231 , we introduced the kernel notation Ax^^Xj^-x'^,x'^, for operators A on S2(/^(Z'^, y)); Ax^,Xj^-x'^,x'^ is 
an element of S§{i3§2 such that 

{S,AS') = ^ {S{xL,XR),Ax^,x^;x'^,x'^S'{x'L,x'f.))^^(^^) (8.3) 

First, we introduce a notion that allows US to bound operators by their 'matrix element r r' r' ■ 

Definition 8.1. Let A and A be operators on 62(^^(2'^, y)) and B2{P{'L'^)), respectively We say that A dominates 
A in 'the sense of matrix elements', denoted by 

A < A, (8.4) 

7n.e. 

iff 

\\Ax^,x^-x'^,x'J[SB{Sg2{y)) < Ax^,x^-x'^,x'^ (8.5) 

Note that, if A is an operator on S2(/^(Z'^)), the inequality A < A literally means that the absolute values of 

m.e. 

the matrix elements of A are smaller than the matrix elements of A. We will need the following implication 

A< A Mil < Mil (8.6) 
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Indeed, for any S e i^2(^^(Z'') (E) -V) ^ x Z'*, ^2(^)), we construct 

S{xl,xb.) := ||5(xi,,a;j^,)||<^2(^) (8.7) 
such that ||5'|lp(zdxZ'i) = l|5'|lp(zdxzrf„'^2(^)) ^"'l 

< (8.8) 

from which Il8.6b follows. 



8.3 Bounding operators 

We introduce operators on ^2('^(^^)) that will be used as upper bounds 'in the sense of matrix elements', as 
defined above. These bounding operators will depend on the coupling constant A, the conjugation parameter 
7 > and the cutoff time r = r(A). Let the function r^(7, A) and the constants c^, be as defined in Lemma [6.21 
and, in addition, let 

r,{j, A) 2A2g,(27), for 27 < S,, (8.9) 
with ge(-) and Se as in Assumption l2.1l We define 

^ ^'^Q-^'^g<!tQ-i\ix'L+^'R)-(^L+XR)\^--f\{xL-XR)-{x'^-x'j^)\ (8.11) 

{Li^)x^,x^:x'^,x'^ e''-('''^)*e-*K'^^+"'«)-(='^+'^«)le-^|('=^-'^«'"('^'^-~'='«'l (8.12) 

In order for definitions l|8.111l8lT2l l to make sense, A and 7 > have to be sufficiently small, such that the functions 
re (7, A) and (7, A) are well-defined. In particular, we need conditions on A and 7 such that Lemma I6l2l applies. 

The operators I^j, Z^'' inherit their notation from the operators they are designed to bound, as we have 
the following inequalities, for A, 7 small enough: 

(8.13) 
(8.14) 
(8.15) 

m.e. 

The first inequality is obvious from the definition of X^^; in | |5.7|| and the fact that || VK|| .®(,y) < 1. Indeed, Xj^ ; can be 
obtained from I^^i by replacing T4^ by 1. The second inequality is the result of Lemma [6.21 and the third inequality 
follows from the bounds following Assumption l2.1l 

We start by stating obvious rules to multiply the operators Z^'' and . 

Lemma 8.1. For A, 7 small enough, the following bounds hold (with c{j) and 0(7, A) as defined in Section \8?]} 

• For all sequences of times si, . . . , s„ with t — ^i' 

UI..MIUI < [c(7)]"-^e=(^-^)*ZY| (8.16) 





< 


^x,l 




m.e. 






< 






m.e. 




Ut 


< 





Zli^ ...Z::^ZIC < [c(7)]"-ie=(^'^)*Z;'' (8.17) 
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• For all times s < t, 

z;lIu: < c(7)e*e^('^-^)*ir* (8.18) 

m.e. 

UlJ:^^ < c(7)e'^e^(^-^)*i:'^ (8.19) 

m.e. 

Proof. Inequalities l l8.16|l and (|8.171 are immediate consequences of the fact that 

Q-r\^-o=i\e-j\^-^2\ < J2 e-*l^l, forany7 > (8.20) 

To derive inequalities I l8.18|l and ll8.19l l, we use jS.lOi and we dominate exponential factors e'^^^^'* on the RHS by 
6 2^, using that rA^ — > as A \ 0. 

□ 

Lemma [8.21 below, shows how the boimds of Lemma [8. II are used to integrate over diagrams. This lemma will 
be used repeatedly in the next sections, and, since it is a crucial step, we treat the following simple example in 
detail: We attempt to bound the expression 

/ dhdt, J2 C{<^)^l^-J^2,iM-tJ^^,hK-s^ (8.21) 

S <tl<t2<S — 

in "the sense of matrix elements", with cr being the diagram in Sj^^i consisting of the ordered pair {{ti,xi,li),{t2,X2,l2 
We proceed as follows: 

1) We bound ({a) by sup^^ IC(o')|. Note that the latter expression is a function of t2 — ti only. 

2) Since the only dependence on xi,X2,li,l2 is in the operators ixi,ii, we perform the sum J2x — 1/ for 
i = l,2. 

3) Since the operators Xx^ have disappeared, we can bound 

l^y,fil-tK-s^ < [c(7)]^e^('''^)l^'-^'l (8.22) 

m.e. 

using Lemma [8. II 

Thus 

< .e^^'^'^)!^'-^'! / dtidt2 [c(7)]' sup |C(a)| (8.23) 



S^<tl<t2<s' 



Note that sup^^^^ ;^ ;^ \Ci<y)\ = sup^ ^2 — ^i) I because |a-| = 1. The short derivation above can be considered to 
be an application of Lemma [7. II as we illustrate by writing 

^x,,x,;x',,x',- J daG{a)F{a), withG(a) = CM, (8-24) 



and hence (|8.23l l follows from Lemma [7.11 after applying (|8.22l l. 
Lemma [8.2l is a generalization of the bound (|8.231 l above. 

Lemma 8.2. Fix an interval I = [s\s^ and a set of m triples (t-, a;^, Z-)™ i such that t[ <E I and t[ < For any 

a e S/(ir), we 
union of triples 



a e S/(ir), we define the set ofn := ni + 2|cr| triples {t'l, a;-', Zf )™^^''^' hy time-ordering (i.e., such that t'l < t"^i) of the 



(t'^xlJ'^Zi, and Ua),x.,{a)M^))f:l (8-25) 
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Then 



J da |C(a)| . . M^_,.i.^,^U^^,.i.^,^ (8.26) 



< 

m.e. 



S/(ir) 
/ 

e 

\ nrSf(ir 



/ d[a] [c(7)]^l^l sup |C(a)|| X U},, I,,,dAl_,/I,^,,Ul_^ 



Moreover, f/ie statement remain true if one replaces Z^''^ on the LHS and -> Z^' ^ on the RHS of ( |8.26t . 

Proof. The proof is a copy of the proof of the the bound (|8.23b . The steps are 

1) Dominate |C(cr)| by sup^(^)_j(„) |C(cr)| 

2) Sumover x{a),l{a)hy using J2xi.h^xiA, =1 

3) Multiply the operators or Z^'^^, using the bound | |8.16l l or ( |8.17t . 

4) Interpret the remaining sum over \a\ and integration over t{a) as an integration over equivalence classes [a]. 

□ 

8.4 Bound on short pairings and proof of Lemma [611] 

We recall that the crucial result in Lemma [Ql see Statement 2 therein) is the bound 



(8.27) 



S[o.t](<T,ir) 



uniformly for Re z > — and for |Im k| small enough. 

In the first step of the proof of | |8.27| |, we sum over the ^(c), {(cr)- coordinates of the diagrams in | |8.27| |. The 
strategy for doing this has been outlined in Sections 18. 2l and[ 8. 31 

Lemma 8.3. For A, 7 smalll enough, 

I dal|,|>2C(a)V[o,t](a) < e^(^'^)* z:/7 / d[a] c(7)l"l 1|.|>2 sup |C(a)| (8.28) 

S[o,f](<T,ir) nTS[o,t](<'r,ir) 

Proof. In the definition of V/((t), see e.g. (|7.16l l, we bound Ix.i by Ix,i and Ut by . Then, we use the bound (|8.261 
with m = to obtain (|8.28l l. Note that, since rn = 0, the set of triples {f- , x'- , l'/)'^"^^"^ is equal to the set of triples 
{ti{a),Xi{a),k{a))f^l. Note also that we use (|8.26b with !]/(< T,ir) instead of ^/(ir), and with the restriction to 
I (7 1 > 2. However, this does not change the validity of (|8.26l l, as one easily checks. □ 

To appreciate how Lemma [8.3l relates to the bound l|8.27l l, we note already that the bound (|8.281 remains true if 
one puts left and right hand sides between J7k • J^-k for purely imaginary k (for general k the matrix elements can 
become negative, which is not allowed by our definition of < ). 

In the second step of the proof, we estimate the Laplace transform of the integral over equivalence classes [a] 
appearing on the RHS of | |8.28|I . This estimate uses three important facts 

1) The correlation fimctions in (|8.28t decay exponentially with rate 1/r, due to the cutoff. 

2) The diagrams are restricted to | it | > 2, they are therefore subleading with respect to a diagram with | ct | = 1 . 

3) We allow the estimate to depend on 7 in a non-uniform way. Indeed, 7 will be fixed in the last step of the 
argument. 
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Concretely, we show that, for < a < ^ and for A small enough (depending on 7) 

die"' / dHl|,|>2 ([c(7)]'l'^l sup |C(a)| 1 -0(A2)0(AV)c(7), A \ 0, AV \ (8.29) 

M+ S(o.tj(<r,ir) ^ ^ 

To verify I l8.29b , we set 

k{t) := A2c(7)l|t|<,sup|V'(x,t)| (8.30) 
and we calculate, by exploiting the cutoff r in the definition of fc( ), 

||e^*fc||i < A2c(7), ||te(^+ll^-|li)'fc||i = rO(A2)c(7) (8.31) 

The norm II • II 1 refers to the variable i, i.e., II II 1 = dt\h{t)\. Hence ( I8.29t follows from the bound (|D.4l l in Lemma 
ID.ll in Appendix IDI after using that tO(A^) < C, as A \ 0, and choosing A small enough. 

In the third step of the proof, we fix 7. By using the explicit form ||8.12| | and the relation (|2.60b , we check that 



<f.cij.^)t^ for any |ImKi^^|<7. (8.32) 



Next, we make use of the following general fact that can be easily checked (e.g., by the Cauchy-Schwarz inequality): 
If, for some 7 > and C < 00, 



II {JkAJ-^,)^^ ,^^.^,^^^,^ II < C, uniformly for k^j^ s.t. llmK^^^I < 7, (8.33) 

then 

\\A\\ < c(7) (8.34) 

where the norms refer to the operator norm on ^{^2(1'^ C^"^, as in Definition 1 8. II 

Hence, from | |8.32| | we get 

sup ||^-«W7J«|| <c(7)e^(^'^)*. (8.35) 

|ImK^^|<7/2 



By the first equality in ( |8.27t and Lemma [8.31 

J.nMJ-. < I Ate-'^^-^jM^J.^e^'^'''^' [ d[a] l|,|>2[c(7)]'l"l sup \(:{<y)\ (8.36) 

"i-e-jR+ J ^((t),U,(t) 

S[o,t](<T,ir) 

We combine l|836l l and (|835t with l l8l29ll , setting 

a = max(-Rez,0) + c(7, A) (8.37) 

for A small enough such that 0(7, A) < ^ and such that (|8.29b applies. At this point, the parameter 7 has been fixed 
and this choice determines the maximal value of |Im k.\. This concludes the proof of the bound in (|8.27l l. The other 
statements of Lemma [O] are proven below. 

Proof of Lemma [01 The claim about TZJ^{z) (Statement 2) follows by a drastically simplified version of the above 
argument for TZl^{z). 

To establish the convergence claim in Statement 1) of Lemma [6.11 it suffices, by l|7.27l l, to check that ||Z(^'"^|| < e^* 
for some constant C. This has been established in the proof of Statement 2), above, since 7?.[^(z) + 'R-l^{z) is the 
Laplace transform of Zj^'". The identity (|6.4|l was established in Section [7!2l 
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To check Statement 3), we employ expression l|C.Hl for C{z) and Il7.28b for TlJ^{z). The latter differs from C{z) in 
that it is the Laplace transform of a quanitity with a cutoff att = T and in the fact that it includes the propagator 

^ gi(ad(Y)+A^ad(e))t whereas C{z) includes only e'f^'i^^))*. We observe that 



\j^{niaiz)-\^ciz))j. 



< 4A^ 
+ 



dt sup 1-0(2;, i) I 

X 

dt sup \ip{x, t)\ 



-iad(Y)t 



J- 



iad(y)t / iA^ad(e(P))t 



I J- 



(8.38) 
(8.39) 



where the factors '4' originate from the sum over li, I2 and we use that Rez > 0. In the first term on the RHS, 
II j7Ke"*"^*'^''*i/-K|| = 1 since Y commutes with the position operator X. The second term is bounded by 



4A2 / ds sup|V(x,s)| X sup (Xh J«,ad(e(P))e'^'^'^(^(-f'»*JL^ 
where we have used Lemma|5]2]and the bound (|2.12l l. 



<tX^C 



(8.40) 
□ 



8.5 Bound on the vertex operators B{1, s\ s^) 

In this section, we prove a bound on the 'dressed vertex operators', which were introduced in Section [7331 Since 
such 'dressed vertex operators' contain an irreducible short diagram in the interval [s\ s^], we obtain a boimd that 
is exponentially decaying in — s'|. In l|8.411 , this exponential decay resides in the function w( ) and it is made 
explicit through the calculation in | |8.45l l. 

The proof of the next lemma parallels the proof of Lemma [8 . 3 1 above . Consider m triples {t[, x'^, Z^)™ ^ and let [ 
be a (dressed) vertex with vertex set S{1) = {t[, . . . , t'^^}- Let s\ s' be vertex time-coordinates associated to [, i.e., 
such that s' < t[ and > t'm- 

Lemma 8.4. For A, 7 small enough, the following hound holds: 

BH,s\s^) < w(s'~s')U}_,i^,^,,Jl^, Ui Z^^jjUi (8.41) 

m.e. ^ ^ m — L Z L i 

where 

m;(/-s') :=e^'^'l"'-*'l / d[(7] Cl"! sup \C{<t)\ (8.42) 

The RHS of | |8.42|| indeed depends only on — s', since the correlation function ({a) depends only on differences of the 
time-coordinates of a. The function w{-) depends on the coupling strength A via the correlation function C(o')/ see l|7.15|l . 

Proof. Starting from the definition of the vertex operator B{{, s', s') given in l|7.42|l . we bound the operators Xx,iMt 
by Ix,i , and we apply Lemma[8]2]to obtain 

B{[,s\s') < e"(^'^'l^'-^'l / d[(T] [c(7)]2|-l sup |C(a)| 

nTS[,i^tj(<r,ir) 

^Lt'J^L,i'Ji^t'_, ■ ■ ■i^l-t'M^vUl_^. (8.43) 

From the definition of in | |8.12|| , we see that 

Qii < e^(T'^)*^Y7^ for 72 < 71 (8.44) 

We dominate the RHS of Il8.43|l by fixing 7/2 = 71 and applying (|8.44l l for any 72 < 71. This yields (|8.411 , with 
the constant C in | |8.42| | given by fixing 7 = 71 in 0(7). One sees that the maximal value we can choose for 71 is 
7i = |;(5e, with 5^ as in Assumption 12. II □ 
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For later use, we note here that, for A sufficiently small and with wit) defined in l|8.42|l : 

/ dt|i|w(i)c^ < tC dtw{t)e- 

< 0(AV2), A\0,AV\0 (8.45) 
where the second inequality follows by the bound l|D.3l l in Lemma [D. 11 with 

k{t) := X^C sup \^p{x,t)\lt<r and a - (8.46) 

X T 

for A such that A^C" < 1 /t with C" as in the exponent of (|8.421 . 
8.6 Bound on the conditional cutoff dynamics Ct{ai) 

In this section, we state bounds on Cf (cr; , £) and Ct (cr; ), defined in Sections l7.2.1l and |7.3.11 respectively. Our bounds 
will follow in a straightforward way from Lemma ISTil and formula (|7.45l l, which we repeat here for convenience 

Ctic7i,£)= [ Ps'Psf S(l|£|,s|£|,s[£|)z;f_^.^ ...z:,_^,B{l2,sl4)z:,_^,Bih,s[,s{) (8.47) 

J 1 32 21 

< sj, < 4 < t 
< s] , for k' > k 

By inserting the bound from Lemma I8l4l in (|8.47l l, we obtain a bound on Ct {<ti , £) depending on the vertex time- 
coordinates s\ s^. In the next bound. Lemma 1531 we simply integrate out these coordinates. To describe the result, 
it is convenient to introduce some taylor-made notation. Let the times {ti, . . . , t2n) be the time-coordinates of ai. 
We will now specify the effective dynamics between each of those times, depending on the vertex partition £. 

• If the times ti and ti+i belong to the vertex set of the same vertex, then 

K+uU ■■= Gl^.-u ' with ^7 (8.48) 

• If the times ti and i^+i belong to different vertices, then 

The idea of this distinction is clear: Within a dressed vertex, we get additional decay from the short diagrams; this 
is the origin of the exponential decay " in Q2 ■ Between the vertices, we encounter the cutoff reduced evolution 
Zl , as already visible in l |8.47|l . Moreover, we get an additional small factor for each dressed vertex. To make this 
explicit, we define 

I ■^^ I dressed '■— #{dressed [fc} {— number of dressed vertices in the vertex partition £) (8.50) 

Lemma 8.5. Let the operators '^Q. ^.^^ he defined as above, depending on the diagram ai and the vertex partition £. Then, 
for A, 7 small enough, 

Ct{<7i,£) < [(|A|r)2c(7)] G7-t,J..„,i.M.,t,^-. ■ ■ ■ nl-tX..M^l-t±.,MQl (8.51) 

m.e. 

Note that between the times and ti, we always (for each vertex partition) put Q'^_^ . This is because either i i = 0, 
in which case Q^_^ = 1, or ti belongs to a dressed vertex whose initial time coordinate, s\, is fixed to be s\ = 0. The 
same remark applies between the times i„ and t. 
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Figure 12: Consider the long diagram ct; G S[o_t] (>t) with \a\ = 4 shown above. In the upper figure, we show a short diagram a such that 
(T; U (T is irreducible in [0,t]. The corresponding vertex partition £ = {li, . . . , Ib} is indicated by vertical lines for the bare vertices [1,(4, 15 and 
horizontal bars for the dressed vertices b, I3, fe- In the picture below we suggest the representation that emerges after applying Lemma 1831 
There are no vertex time coordinates any more. The time-coordinates of the long diagrams correspond to operators X. The intervals between 
time-coordinates of the long diagrams correspond to operators Z'^''^ or gT. The intervals corresponding to Q"' are those which in the upper 
picture belong entirely to the domain of a short diagram. 



Proof. Theproof starts from the representation of Ct(cr;, £) in (|8.471 and the bound for the vertex operators B{lk, s),,, 
given in Lemma [8!4l Then we integrate out the s'j. , s^-coordinates for the dressed vertices . The main tool in doing 
so is the fast decay of the function w( ), as follows from | |8.45| |. 

We consider a simple example. Take ti = and = t and let |£| = 1, i.e. there is one vertex I. It follows that I 
is dressed and S{Vj = {t2, . . . , i2n-i}- hi this case, formula | |7.45| | reads 

Ct{ai,&)= J dsWl,,^,i„z;_^,Bil,.s\.s^)Zl^^I,,^i, (8.52) 

< s' <t2 
t2n-\ < <t 

and the boimd in Lemma |8^ is 

Bil,s\s') < w{s' -s')W,^^,^^_^xI.,,^_,,i,^_Ml^_^_,^^_^.^ (8.53) 



where the operator is defined as the 'interior part' of the vertex operator. The sole property of that is relevant 
for the present argument is that 

JC < e'=('''^)*^* (8.54) 

m.e. 

as follows from the definition of and the bound ll8l6l l. From | |833|| , l(834l l and (|8l8ll819t , we obtain 

Ct(CT,,£) < e~ ""'^''''" / ds'd/u;(/-s')e^(c(7))^ (8.55) 

m.e. J 

< s' < t2 

t 

-^X2n ,l2n ■^t2n~t2„-l-^ 2 - 1 1 "^aj i , i i 
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where we have used the decomposition — s' = (s^ — i2n-i) + {hn-i — + {h — s') and we have chosen A, 7 
small enough such that 0(7, A) < 1/ (6t) in (|8.54l l. By a change of integration variables, we find that 



/ 



dsM/w(/ - s')e^^ < / dt |t|u;(t)e* (8.56) 



0<S' <t2 R+ 

t2n-l < < t 

and we note that this bound remains valid if, in the integration domain on the LHS, we replaced by a smaller 
number, or by a larger number. Hence, by the bound | |8.45l l, we obtain 

CM,£) < [c(7)]^(|A|r)2e-^^^^^I.,„,,„i,;':_,^_^^^i;-Z,^i.,,, (8.57) 

m.e. 

where the constant that originates from the RHS of the bound i8A5i has been absorbed in 0(7). The bound | |8.57| | is 
indeed I |8.51| | for our special choice of £ in which |£|dressed — 1- To obtain the general bound, one repeats the above 
calculation for each dressed vertex. These calculations can be performed completely independently of each other, 
as is visible from the remark below | |8.56|I . 

□ 

In Lemma [8.51 the bound depends on £ through H'', see (|8.481 and (|8.491 . The next step is to sum over £,. First, 
we weaken our bound in (|8.51t to be valid for all £, such that the sum over £ amounts to counting all possible 
£ ~ (7;. By "weakening the bound", we mean that we bound some of the operators by Zp'' . This can always be 
done, since, for A small enough, 

57 < Zl'^ (8.58) 

m.e. 

with Ql as in ll8.48l l (in fact, Q2 is smaller than the second term of Zl''^ , see (|8.11l l). Let cti, . . . , am be the decom- 
position of (Ji into irreducible components and let S2i-i, S2i be the boundaries of the domain of ai. These times 
Si should not be confused with the vertex time-coordinates , £ that were employed in an earlier stage of our 
analysis. In particular, the times S2i-i, S2i, i = 1, . . . m, are a subset of the times t^, i = 1, . . . , 2n. The central remark 
is that 

For any i, the times S2i, S2i+i belong to the same vertex for all vertex partitions 2, ^ ai. 

Indeed, since the interval [s2i, S2i+i] is not in the domain of cr;, it must be in the domain of any short diagram 
contributing to Ct{(Ji), or, in other words, any vertex partition £ ^ ct/ must contain a vertex whose vertex set 
contains both S2j, S2i+i- Consequently, the operators "HJ^. ^^i+i I l8.51| | are always (i.e., for each compatible vertex 
partition) equal to ^J^. s^i+i' ^^'^ ^^^^ replace them. However, we replace all other T-Q. i.e. those with 
the property that the times tj , ij+i are in the domain of the same irreducible component of cr/, by 2/^'^^ ■ 
This procedure is illustrated in Figure [131 

After this replacement, the operator part of the resulting expression is independent of £, and we can perform 
the sum over £ ^ ct/ by estimating 

E [(|A|r)2c(7)]"""^'^" < (|A|t)2-(-')c(7)I'^'I, for |A|r < 1 (8.59) 

with 

v{ai) := min |£|dressed (8.60) 

To obtain | |8.59|I , one uses that 

#{£- CT,} < 42|'^'l-i (8.61) 

Indeed, 2^1°^' is the number of ways to partition the time-coordinates into vertex sets. The extra factor 2 for each 
vertex takes into account the choice bare/ dressed. 
We have thus arrived at the following lemma 
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Z^'~< Z^-'< Z^'~< Q~< Z^-'< Z^'~< Z^'-i Q-i 



Figure 13: Consider the long diagram cr; £ S[o_t] (> t) with |o-[ = 4 shown above. It has two irreducible components with domains [si , 52] 
and [53,34]. In the upper figures, two different vertex partitions (compatible with ai ) are shown together with their respective bounds, obtained 
in Lemma 1831 These bounds are represented by the operators Q^^ and .E^'^, as in Figure [121 except for the fact that we omit the operators X 
corresponding to the time coordinates of o-j . In the lower figure, we show the (weaker) bound that gives rise to Lemma |831 To establish this 
weaker bound, we replace the Q"< that are 'bridged' by the long diagram by .E^-^. 



Lemma 8.6. Let S2i-2, S2i the boundaries of the domain of Ui, the i'th irreducible component ofai. Then, for A, 7 small 
enough, 

Ct{oi) < (|A|r)2"(-')^7_,^,„£^(a™)^J^^_^_,^^_/~^(a™_i)...f^(ai)g,\ (8.62) 

ni.e. 

where, for an irreducible diagram a with \a\ = p, 

£^{(j) [c(7)]''''^x2p(^),i2p(^)4''J(<x)_t2p-i(<T) • • • ^^t2ja)^t,(a)^M'y)M{<y) (8-63) 
with v{ai) as defined above. 

Note that v{<7i) is actually the number of factors G2 ™ the expression (|8.62t for which u ^ {u can be zero only 
for the rightmost and leftmost GZ)- Or, alternatively, 

v{<7i) = #{irreducible components ino-;} — 1 + lt2„^t + ^ti^o (8.64) 



8.7 Bounds on TZpjfz) in terms of S^'(a) 



To realize why the bound Il8.62b in Lemma[8]6]is useful, we recall that our aim is to calculate TZex (z), given by (see 



daC(a)Ct(a) [I + 5{t^{ai)] [l + 5{t^\„^\{ai) - t)] 



(8.65) 



S[0.t](>T) 



We calculate TZex {z) by replacing the integral over diagrams by an integral over sequences of irreducible diagrams, 
as we did in (17. 231 , i.e.. 



n 



— 0<si <...<S2„<i 



S[0,t](>T) 

Using the bound | |8.621 l, we obtain, with the shorthand z.,. := Re 2;, 

Tiex{z) < ^(l + 7^g(z,))7^£~(z,)(7^g(z,)7^£^(z,))"(l + 7^g(z,)) 



(8.66) 



il+ng{Zr))n^{Zr) (1 - ng{Zr)n^{Zr)) ' (1 +7^g(z,)) 



(8.67) 
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where, for Re z large enough. 



7^.(z) := (t|A|)2 / dte-*^^7 (8.68) 



7^£(z) := / dte^'^ / AaC,{a)e^{a) (8.69) 

S[o,t](>T,ir) 

Since ^7 is kriown explicitly, the only task that remains is to study TZg{z). This study is undertaken in Section|9l 

9 The renormalized model 

In this section, we prove Lemma [631 thereby concluding the proof of our main result, Theorem |3.3l We briefly recall 
the logic of our proof. As announced in Section |5.4.1[ we analyze Z" and Zt through a renormalized perturbation 
series, where the short diagrams have already been resummed. However, we do not study the Laplace transform 
of irreducible diagrams (defined in H7.20I I) 

TZ^iz) = J dtc-'''Zl' ^ J dte-*' J daC{a)Vt{a), (9.1) 

K+ R+ S[Q,tj(ir) 

directly, but rather the Laplace transform of irreducible renormalized diagrams (defined in Lemma [8 . 6 1 and Section 



= J dte-''' J d<7C{<7)£''{a). (9.2) 

K+ S[o,tj(>r,ir) 

Although the quantities ( 19. 2t and | |9.1t are not equal, we will argue below (in the proof of Lemma 16.51 starting 
from Lemma |9.1| l that good bounds on Tl£{z) yield good bounds on TZ"{z) and hence also on TZ{z). The reason 
that the expression ( 19. It itself cannot be bounded by an integral over long irreducible diagrams, is the fact that an 
irreducible diagram in the interval [0, t] does not necessarily contain an irreducible long subdiagram in the interval 
[0,t]. Indeed, Lemma 18.61 decomposes the domain of an irreducible diagram into domains of long, irreducible 
subdiagrams and intermediate intervals. These remaining intervals give rise to operators , which are easily 
dealt with, as we will see below, in the proof of Lemma l63l since they originate from short diagrams and therefore 
have good decay properties. 

Nevertheless, we clearly see the similarity between ( 19. It and i9.2i . To highlight this similarity, we write the 
inverse Laplace transform of 7?.^ (z): For r > large enough, we have 

^ J dze*^7^^(z)= J dac(7)l^kM^..„.b„i;:-t.„_,---^...;.it"2A,h (9-3) 

r+iR S[o,tj(>T,ir) 

where t,x,l are the coordinates of cr and, since a is irreducible in [0,t], ti = Q and t2n = t. The inverse Laplace 
transform of TZ"{z), i.e. Z", is 

J daC(a)I,,„,,„Z^,^„_,^_^...I,„z,Z^,^_,^I,„z, (9.4) 

S[o,t](ii') 

where t, x, Z have the same meaning as above. Thus, the perturbation series in (19.31 1 is indeed a renormalized version 
of (Hill. The dia grams are constrained to be long, and the short diagrams have been absorbed into the 'dressed free 
propagator' Zj'''. This point of view has also been stressed in Section l531 Observe, however, that Z^'"' depends on 
the positive parameter 7, whereas there is no such dependence in ( 19. 11 1. 
The following lemma is our main result on 7?.^(z). 
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Lemma 9.1. Recall that TZgiz) depends on 7, because £''{■) does. One can choose 7 such that there are positive constants 
gex > and Sex > 0, such that 

sup |lJ;7^£-(z)J_„|l =o(A2), flsA\0 (9.5) 

|ImK^ ^ I <(5eK ,Rc — A2(7g^ 

The main tools in the proof of Lemma [9. II will be the exponential decay of the 'renormalized correlation func- 
tion', which follows from the bounds on stated in Lemma 16.21 and the strategy for integrating over diagrams 
presented in Lemma [8.21 With Lemma |9?1] at hand, the proof of Lemma |631 is immediate. 

Proof of Lemma\6^'We only need to prove Statement 2) since Statement 1) will follow by a remark analogous to that 
in the proof of Statement 1) of Lemma fe.ll Clearly, for A small enough, 

sup \\J.ng{z)J^4<0{X), asA\0 (9.6) 

This follows from the properties oiU^ , see e.g. the proof of Lemma [6. 11 and the definition of Q] , see ll8.48ll . Next, 
we remark that 

WJ.nMJ-d < U + JnT^g{Zr)J-n\\^\\J.ng{Zr)J-^\\\\ (l - J«7^g (z,) || (9.7) 

with Zr ~ Re z. This follows from the bound l|8.67b , the fact that J^J-k — 1/ and the implication (|8.6b (which allows 
to pass from ' < ' to an inequality between norms). 

m.e. 

Hence, Statement 2) follows by plugging the bounds of Lemma [93] and l |9.6t into the the RHS of | |9.7t . □ 
9.1 Bound on the renormalized correlation function 

In this section, we prove Lemma [9.21 which establishes (as its first claim) the property (|5.23l l with Af replaced by 
ZJ^ . Indeed, in Section |6Al we argued that Zl is very close to A^, and this was made explicit in Lemma [6l21 Let 

( e-(i/2)ffRt \x\/t<v* 
hit) := X^ch sup < (9.8) 

a^ez" i \i>{x,t)\ \x\/t > V* 

with the velocity v* and decay rate gji as in Lemma [5. II and the constant Ch chosen such that 

sup \ip{x, t)\ < h{t), fort>T (9.9) 

X 

Lemma [S.ll ensures that such a choice is possible. 

Lemma 9.2. There are positive constants 6r > and gr > such that, for all j < 5r, X small enough, and n = (kl, k^) 
satisfying \lmKj^j^\ < ^, 

X^mxl-xi^,t)\ X \\{jar''J-.)x^,x,-x'„x'J < h{t)e-''"ar\ for hM € {L,R} (9.10) 

and 

{J.Zl^'J-.)x,,x,;x',,x'j\ < Ce^^'''^' (9.11) 

This lemma is derived from the bound (|6.7l l in Lemma [6.21 in a way that is completely analogous to the proof 
of Il5.23b starting from ll4.55b , as outlined in Section [5. 4. II The only difference is that in Lemma [9.21 we allow for a 
small blowup in space given by the multiplication operator J7k. 

For future use, we also define 

Kit) := lit\>rhit) (9.12) 

and we note that 

\\hr\\i:= [ hrit)=oiX^), asA\0 (9.13) 
since \\h\\i < oo and t(A) — t- oo as A \ 0. 
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9.2 Sum over non-minimally irreducible diagrams 

In a first step towards performing the integral in (|9.2b , we reduce the integral over irreducible diagrams to an 
integral over minimally irreducible diagrams. Indeed, since any diagram that is irreducible in / has a minimally 
irreducible (in /) subdiagram, we have, for any positive function F, 



J dcrF{a) < J da F{(j) + j Aa' F{<7 U a') 

S/(ir,>r) Sj(mir,>r) \ ^l{>r) ) 



(9.14) 



The first term between brackets on the RHS corresponds to the minimally irreducible diagrams on the LHS. The 
second term contains the integration over 'additional' diagrams a' . The integration over these diagrams is uncon- 
strained since u U ct' is irreducible in / for any a' , provided that a is irreducible in /. This is also explained and 
used in Appendix O see (|D.5l l and (|D.6l l. 

Lemma [9.31 shows that such an integration over unconstrained long diagrams yields a factor exp{c'(7, A)|/|}, 
with the generic constant c'(7, A) as introduced in Section lHTTl 

Lemma 9.3. For A, 7 small enough, 

[ da£''{a)\Cia)\ < e"'^''^^^' [ da \Ci<j)\£i (a) (9.15) 

J me. J 

S[o,t](>T,ir) S[o.t] (>r,mir) 

Proof. By formula l|9.14l l (applied in the case where F{a) is a matrix element of the operator \({a)\£'' (a)), we have 
that 

J dcrf'^(cr)|C(a)| < J da r {a)\Cia)\ + J da J da' T (a U a')\C{a U a')\ (9.16) 

S[o,t] (>-r,ir) S[o,tj (>T,mir) Ejo.tj (>T,mir) i;[o,t](>-r) 

First, we bound 

J da' r{a[Ja')\Cia(Ja')\ (9.17) 

S[o.t](>r) 

with a fixed. To perform the integral over a' in l|9.17l l, we recall that £''{■) consists of products of the operators 
Ixidi and Z/^'J^_j. . Hence, by Lemma |8^ with replaced by Z^''', we can sum over the x, {-coordinates of a' and 
multiply the Zl^'^^_^. operators using the boimd l|8.17l l. This yields 

ei3 < £;^(a)|C(a)|e^(^'^)* / d[a']c(7)l^'l sup |C(a')|. (9.18) 

■m.e. J x(cr'),!_(cr') 

nTSt(>r) 

The integral on the RHS of | |9.18| | is estimated as 

d[a']c(7)l'"'l sup \C{a')\ < e"^'')"''-"^* - 1 (9.19) 



a^(cr'),/_(cr') 
nTSt(>r) 

with hr as defined in | |9.12|| . This follows from the bound l|D.8|l (integral over unconstrained diagrams) in Appendix 
[Dl To bound the first term in | |9.16| |, we dominate 

r{a) < e'='('''^)*f^(cr) (9.20) 

m.e. 

The lemma follows by inserting the bounds i9.2Qi and (I9.18ll9.19l in (|9.16t and using that c(7)||/z^||i = cJ{j, A) since 

WhrWl ^ 0{X'). □ 



63 



9.3 Sum over minimally irreducible diagrams 

In this section, we perform the integral 

J da\aa)\£'<ia) (9.21) 

St(>r,mir) 

that appears on the RHS of the bound in Lemma |93l (upon replacing 7 — > 7/2), and we prove that it is exponentially 
decaying in time with decay rate O(A^), for well-chosen 7 and A small enough, depending on 7. It is in this place 
that we use the decay property of the renormalized correlation function that was stated in Lemma |921 

The key idea is the following. If <t is a long diagram with |cr| = 1 consisting of the two triples (ti, xi, then 

\aa)\£\a)=c{^)\^*{x^^x,M-ti)\ X,,,,, J^^Z, I,,,,, (9.22) 

In this case, we can obviously use Lemma 19.21 to deduce exponential decay in t2 ^ of Il9.22b , uniformly in 
xi, a;2j ^1: '2- In general, there is of course more than one pairing in an irreducible diagram and so one has to 
'split' the decay coming from Ix2d2^t2-ti^xidi between the different pairings, thus weakening the decay by a 
factor which can be as high as 

However, since we are considering minimally irreducible pairings, there are at most two pairings bridging any 
given time t', see Figure [141 Hence, one can attempt to split the decay from Ir^2,i2^t2-ti'^xi,h in half. This can be 
done and it is described in Lemma [9.41 

For CT e S[o_t], we define the fimction 

kl 

H,{a)^Y[hr{v,-Uj) (9.23) 

i=i 

with hr as in | |9.12| |. Note that Hr{a) depends only on the equivalence class [a], and hence we can write H.r{[a]) := 

Hr{(j). 

Lemma 9.4. Let the positive constants gr and Sr be defined as in Lemma \92\ and choose 7 < Sr- Let k — {kl, k^) such that 
jlniK^^I < ^ / 8 and fix a long irreducible diagram class [a] e nTS[o.t](> T,mir). Then 



J2 \a'T)\jj''ia)j-^ 



< c(7)l"le"(^'^)*e-^^'9'-*i/,([r7]) (9.24) 



Note that the operator between \\ ■ \\ depends on the equivalence class [a] only, due to the sum over x{a),l{cr). 

Proofi For concreteness, we assume that |<t| = n is even, the argument for \<t\ odd is analogous. We can find cti 
and a2 such that cri U (T2 = f, Wi\ ~ \(T2\ = n/2 and ai and (T2 are ladder diagrams, i.e. their decompositions into 
irreducible components consists of singletons. To be more concrete, the time-pairs of ai are 

(ti^ts), (t4, t^), (^8, til) ■ ■ ■ (^2n-4, hn-l) (9-25) 

and those of ct2 are 

(^2,^5), (^6,^9), (^10,^13) • ■ • (i2n-2, hn) (9-26) 

The possibility of making such a decomposition is a consequence of the structure of minimally irreducible dia- 
grams, as illustrated in Figure [141 

We now estimate the LHS of | |9.24| | in two ways. In our first estimate, we take the supremum over the x,l- 
coordinates of fT2 and we keep those of ui. In the second estimate, the roles of <ti and <T2 are reversed. 
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Fig ure 14; The decomposition of a minimally irreducible diagram cr into two 'ladder diagrams' (j\ and u'l. In the upper figure, one can easily 
check that any point on the (horizontal) time-axis is bridged by at most two pairings. 



We estimate 

E ICMI^~"M (9-27) 

x{ai).l(pi) a;((Ti),i(c7i) 

< [c(7)]l"le^(^-^)*[c(7)]l^^l ( sup |C(a2)|l 

m-e. \x(a.2).l(a^) j 

E K(o'l)|-2t2j^-.f2„_jIa;2n-l,'2,.-l^t2„^-l-t2n-4 " ' -^^ih^t-i-t^xx 
x{rTi),l(ai) 

The equality follows from the definition of £''{cr)- To obtain the inequality on the third line, we perform the sum 
over x{a'2) , 1{(J2) by the same procedure that was used to obtain l|9.18t , i.e., by using Lemma 18.21 Since |C(o'i)l 
factorizes into a function of the pairs in ai, the operator part in the last line of i9.27} is a product of two types of 
terms, namely; 

E - a;4», i4i+3 - i4i)l ^'a:4i+3,i4i+3^Ui+3-t4.^^4.,Z4, (9-28) 

la, hi+s 

for i = 1, . . . ,n/2 — 1 (and an analogous term where we replace Ai ^ 1 and Ai + 3 ^ 3, corresponding to the first 
pair in cti) and 

^^4'it4.-x (9-29) 

for i = 1, . . . ,n/2. 

We note that Lemma [9.2| provides a bound on the matrix elements of these expressions. In particular, we use 
| |9l0l l to bound and | |9ll] l to bound i929\i . We obtain, for |Im k^^] < ^, 



(J.(|929i^-.),^,,^^,,y^ < e=(^-^)(*--*--) (9.31) 

By the relation stated in | |8.34| | and the line following it, we can convert these boimds on the kernels into bounds 
on the operator norms, yielding, for |Ini < ^7 

\\J.^MJ-k\\ < c(7);i,(t2.+i-t2»-i)e-^'»'-(*-+^-*-) (9.32) 
\\J^^MJ-n\\ < c(7)e^('''^)(*«-*^'-^) (9.33) 
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and hence, by multiplying these bounds for the operators appearing in (|9.27b and using that 



n/2-l 

hrits - ti) hr{W+z - Ui) < Hr{ai), sup |C(0'2)| < Hr{<J2), 

£(o'2), 1(0-2) 



(see | |9.9| |), we arrive at 



\JAWMJ-^\\<ch)^"^Hr{a)e 



-A^g^lDomo-il c(7,A)t 



(9.34) 



(9.35) 



The claim of the lemma now follows by applying the same boimd with the roles of cti and a2 swapped, taking the 
geometric mean of the two bounds and noting that 



[0,t\< |Domcri| + |Domcr2| 



(9.36) 
□ 



Next, we use Lemmata 19. 3l and 19. 41 to prove Lemma [9.11 By these two lemmas, the integral over renormalized 
irreducible diagrams is reduced to an integral over minimally irreducible equivalence classes [a] . Each equivalence 
class [a] essentially contributes c{'^y'^^Hr{cr) to the integral. Since Hr{(y) is not exponentially decaying in Domcr, 
the Laplace transform of HT{a) cannot be continued to negative Re z. However, the factor c^^^ha^* in Lemma [9.41 
enables us to do such a continuation since the factors c'(7, A), 0(7, A) from Lemmata l9.3l and l9.4l can be made smaller 
than A^ by first choosing 7 small enough, and then adjusting A. 

Proof of Lemma \9J\ We choose 7 small enough, as required in the conditions of Lemma ta l9 . 3l and l9 .41 and we estimate, 

for |ImK^ o| < 7/8, 



da\aa)\jJ-'{a)J. 



S[o,t] (>-r,ir) 



< e' 



c'(7,A)t 



S[o,t](>i",mir) 



'.37) 



^ \aa)\jJHa)J^^ 



nrSjo.tj (>T,mir) 

nTS[o_t](>T,mir) 



The first inequality is Lemma 19.31 the second inequality uses the definition of the measure da, and the third in- 
equality follows from Lemma [9.41 

We will now estimate the Laplace transform of the integral in the last line of l|9.37b . To prove Lemma [9Al we fix 
dex '■= gffr and we show that one can choose 7 such that, for A small enough and Kez > —X^Qex, 



j d[(j]c(7)l'"le-3^'f'-*i/^(cr) =o(A2), asA\0 (9.38) 



nrS[o,t] (>T,mir) 



Of course, the choice of 7 will have to depend on the specific value of c'(7, A). To show l|9.38b , we first choose 7 
such that, for A small enough, we have 

c'(7,A)-iAV<-^AV (9.39) 



withe' (7, A) asinl|938ll. 

Consequently, we can dominate the factor e-^^e'-'^''-i-^)-\^ 3r)t by 1 jji ||9.38t . Next, we note that, if a is minimally 
irreducible in the interval [0, t], then 

kl 

^\v,~u,\<2t (9.40) 

i=l 
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where {ui,Vi) are the pairs of time-coordinates associated to cr. This follows from the observation that each point 
in the interval [0, t] is bridged by at most two pairings of cr, see also Figure [ill Consequently, we find that 

kl 

il^(cr)e-3^'f'-* < Yl K{vj - u^Oe-s^'s-l"^"-"^"! (9.41) 
i=i 

We estimate the LHS of (|938l l, with e-^*e('='('>''^)"^^'^''-)* replaced by 1, by invoking l|a2)l in Lemma IdH with 

fc(0 := c(7)c-^^'9'-*/irW and a := (9.42) 
hideed, by using — o(A^) and the exponential decay e"^'*'^^'-*, we obtain that 

II fell 1 = c(7)o(A2) (9.43) 

||ifc||i = c(7)o(|An, asA\0 (9.44) 
Therefore, the bound (p!2)l yields l|938t . □ 

9.4 List of important parameters 

We list some constants and functions that we use throughout the paper. We start with the different decay rates; 
in the third column we indicate where the constant appears first. By "full model", we mean: "the model without 
cutoff". 



9K 


bare reservoir correlation fct. (for subluminal speed) 


Lemma 15. 11 




bare reservoir correlation fct. in the cut-off model 


Section |53j 




renormalized joint S — R correlation function 


Lemma 19.21 


A Qrw 


Markov semigroup 


Proposition 14.21 




cut-off model 


Lemma 16.21 


A^gex 


excitations in the full model 


Lemma 19. 11 




full model 


Theorem 13. 21 



Additionally, the rates grw, dc 9 come with a superscript low, high indicating that the gap refers to small, large 
fibers p, respectively. 

The following constants restrict the values of complex deformation parameters, in particular the parameter k in 
Jk, as defined in l2.59l 



6e 


particle dispersion law 


Assumption 12.11 




reservoir dispersion law 


Assumption 12.31 


6 


full model 


Theorem |33j 


^rw 


full Markov semigroup 


Proposition 14.21 




excitations in the full model 


Lemma 16.51 


5r 


renormalized S — R correlation fct. 


Lemma 19.21 



The following functions of 7, A appear as blowup-rates in exponential bounds. 



rrw{j, A) 


Markov semigroup 


Lemma 18. 11 






Sectionl8.3l 






Lemma 16.21 



In the final Sections [8] and |9l these rates are represented by the generic constants 0(7, A), c' (7, A), introduced in 
Sectionim 
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A Appendix: The reservoir correlation function 

In this appendix, we study the reservoir correlation function ip{x, t) and we prove Lemmata 15. ll and lS. 21 Recall the 
definition of the "effective squared form factor" if; in (|2.27l l. It is related to ^'(a;, t) by, see II5.3I I, 

'<p{x,t) = j j dsVi(a;)e'*"e'"'"^ (A.l) 

R gd-i 

From this expression, one imderstands that ^'(a;, t) cannot have exponential decay in t, uniformly in x. One also 
sees that, for x fixed, there is exponential decay provided that ) is analytic in a strip around M. 
Let q{-) be the Fourier transform of if), then 

%l){x,t)= i dsq{t + s-x) (A.2) 



By Assumption |2.3[ there is a (5r > such that q{t) decays as Ce '''^1*1 . Choosing v* = i, we obtain, for |a;|/|t| < v*, 

<e-»«l*IC, with.gR (A.3) 

which proves Lemma f5.ll From now on, we assume that > v* . 

We remark that jA.li can be rewritten, after an explicit calculation, as 



d-3 



^p{x,t) = J dr]qit + r]\x\)a{r]), a{r]) := Volume{S''-^) {l - r]^) ^ (A.4) 

By Assumption 12. 3[ in particular the condition t/'(0) = and the analyticity of ip, we deduce that is analytic 

in a strip around M, as well. Its Fourier transform, Q, is an exponentially decaying C^-function (since ^^^^ G L^) 
whose derivative equals q. 

Hence, by partial integration, the fact that a{ri) = a{jf) | = for d > 3, and the change of variables C, = \x\ri, 
we obtain 

. + ^«'(^), Q' - ^ (A.5) 

Here, Q' and a' stand for the derivatives of Q and a. We evaluate this integral by splitting it into the regions 

- |a:| < C < -|a^l + 1, -\x\ + l<C<\x\-l, |a;|-l<C<|x| (A.6) 

In the second region, we dominate the integral l|A.5|l by ||Q||i x II j^jryrTjo'l |§| )l|oo (we assume here that |a;| > 1, 
otherwise the decay in t has been proven above). In the first and third region, we dominate the integral l|A.51 by 
IIQIloo X II |-j^p^a'(|§|) 111. Using the explicit form of a' and the fact that Ixl/lt I > w*, we conclude 

sup|V'(x,t)| <C(l + |^|)"^/^ ford>4 (A.7) 

X 

which implies Lemma I5l2l Obviously, dispersive estimates l|A.7|l can be derived in much greater generality, see e.g. 

GHl. 



B Appendix: Spectral perturbation theory 

Let e e K be a small parameter and consider a continuous function R+ 3 t ^ V{t, e), taking values in a Banach 
space, and such that 

supe"*"||F(t,e)|| < oo, for some to > 0. (B.l) 
t>o 
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The Laplace transform 

A{z,e):= [ dte-*'V{t,€) (B.2) 
is well-defined for Rez > m and it follows (by the inverse Laplace transform) that 

V{t, e) = / dz e^*A(2, e), with := m' + IM for any m' > m (B.3) 

2771 J 

where the integral is in the sense of improper Riemann integrals. 

We will state asstimptions that allow to continue A{z, e) downwards in the complex plane, i.e., to Re 2; < m and 
to obtain boimds on V{t, e). 

Lemma B.l. For Re z large enough, let 

A{z,e):={z-iB- A^{z,e))-^ (B.4) 

and assume the following conditions. 

1) B is hounded and its spectrum consists of finitely many semisimple eigenvalues on the real axis, that is 

B= J2 ^lf(^)' (B.5) 

where lb{B) is the spectral projection corresponding to the eigenvalue b. For concreteness, we assume that e spB. 

2) For e small enough, the operator-valued function z Ai {z, e) is analytic in the domain Rez > —eg a and 

sup \\Ai{z,e)\\ = 0(e), (B.6) 

Rez>-£gA 

sup \\^Ai{z,e)\\ = o{\e\% e\0 (B.7) 

Rez>-egA '^^ 

3) There are bounded operators N)„for b e spB, acting on the spectral subspaces Ranlfe(B) and such that, for all b e spB, 

eNi, - h{B)Ai{ib,e)lb{B) = o(e), e \ 0. (B.8) 

Consider the operator 

N:= (B Nb, with[B,N]=0, (B.9) 
and assume that N has a simple eigenvalue /jv such that 

spN = {/jv} U fijv und supRefiiv < —Qn (B.IO) 
for some gap qn > 0. We also require that 

RefN>-gN, RefN>-gA (B.ll) 

The eigenvalue fN is necessarily an eigenvalue of Nb for some b e spB. For concreteness (and to match with our 

applications), we assume that it is an eigenvalue of Nq 

Then, there is an cq > such that, for \e\ < eo, there is a number /(e), a rank-one operator -P(e), bounded operators 
R{t, e) and a decay rate g > 0, such that 

V{t, e) = P(e)e^('^)* + R{t, e)e-^s* (B.12) 

with 

f{e)-efN = o(e) (B.13) 

||P(e)-l/„(iV)|| = o(|e|°) (B.14) 

sup ||i?(t,e)|| = 0(|e|°), as |e|\0 (B.15) 
teM+ 
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with 1 {N) the spectral projection of N associated to the eigenvalue f^. The decay rate g can be chosen arbitrarily close to 
mm{gN,gA} by making eo small enough. In particular, one can choose g and eo such that Re /(e) > —eg for all \e\ < eg. 

If in addition N and A\ depend analytically on a parameter a in a complex domain V c C, such that l lB.6ll - l|B.7|l -l lB.8ll - 
||B.101 - I|B.11|| hold uniformly in a e V, then (|B.12|l holds with f, P and R analytic in a and the estimates (|B.13l l- (|B.14|l - l|B.151 
are satisfied uniformly in a eV. 

Lemma IB. II follows in a straightforward way from spectral perturbation theory of discrete spectra. For com- 
pleteness, we give a proof below, using freely some well-known results that can be found in, e.g., [26J. 

Lemma B.2. The singular points of A{z,e) in the domain Rez > —egA He within a distance of o{e) of the spectrum of 
iB + eN (provided that there are any singular points at all). 

Proof. Standard perturbation theory implies that the spectrum of the operator 

iB + Ai{z,t), iovYlez>-egA (B.16) 

lies at a distance 0(e) from the spectrum of \B. Here and in what follows, the estimates in powers of e are uniform 
for Re z > —eg a- Let 1^ = lb{B) be the spectral projections of B on the eigenvalue h. As long as e is small enough, 
there is an invertible operator U = U{e, z) satisfying — l|j = 0(e) and such that the projections 

Ifc C/lgC/-\ bespB (B.17) 

are spectral projections of the operator l iB.161 associated to the spectral patch originating from the eigenvalue b at 
e = 0. It follows that the spectral problem for l lB.16t is equivalent to the spectral problem for 

^ U-Ht {iB + A, [z, e)) hU - + + Ae.,fc(z, e)) (B.18) 



(B.19) 



where 

A,,A^,e) := lgC/-i(iS)C/lO-i6lO, {0{e^)) 

+ elgC/-i7VfcC/lO - eNb, {0{e^)) 

+ l'',U-'{A,{ib,e)~eNt)Ull (o(e)) 

+ l<',U-HAdz,e)-A,iib,e))Ull (|z - i6|o(|e|0)) 

The estimates in powers of e are obtained by using U — 1 ~ 0(e), the property li,U ~ J71° and the bounds I IB.61 1- 
| |B.7|I -| |B.8|I . When z is chosen at a distance 0(e) from ib, then all terms in I IB.191 1 are o(e). The claim of Lemma [B.2I 
now follows by simple perturbation theory applied to the RHS of l|B.18| l. 

□ 

Lemma B.3. The function A{z) has exactly one singularity at a distance o{e)from ef^. This singularity is called f = /(e). 
The corresponding residue P = P{e) is a rank-one operator satisfying 

||/'-lAr(/iv)|| =o(|e|°), e\0 (B.20) 

Proof. By Lemma IB.2I there can be at most one singularity. We prove below that there is at least one. By the 
reasoning in the proof of Lemma [B.2l and the fact that the eigenvector corresponding to Jm belongs to Ranl^^Q (see 
condition 3) of Lemma fB.lb , it suffices to study the singularities of the function 

z^ {z-eNo + A,,fl{z,e)y^ (B.21) 

Let the contour = [e) be a circle with center e/w and radius er for some r > 0. Clearly, for r small enough, the 
entire spectrum of eA^'o lies outside the contour F-', except for the eigenvalue c/at. The contour integral of {z — eN)~^ 
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along equals the spectral projection corresponding to /at. We estimate 

/ dz [{z - eNo - A,,.o{z, e))-' - {z - eNo)-'] (B.22) 
Jrf 

dz [z - eNQ)-^A,,^o{z, e)(z - tN^ - A,,^q{z, e))'^ (B.23) 



rf 



= (f dz{€-^c{r)fo{€), as e \ with c(r) sup ||(2 - iVo)"^||, (B.24) 

Jrf \z-f^\=r 

The last estimate holds in norm and it follows from the bound ||yle2,^o(2, = o(e), see (|B.19ll . The expression llB.241 
is 0(1), as e \ 0, since the circumference of the contour is 27rre. From the fact that the contour integral of l lB.211 
does not vanish, we conclude that A{z) has at least one singularity inside F-^. 

The claim about the residue is most easily seen in an abstract setting. Let F{z) be a Banach-space valued analytic 
fimction in some open domain containing 0, and such that e spF(O) is an isolated eigenvalue. We have the Taylor 
expansion 

Fiz)^y^—F„, F„ :=i^(")(0), OespFo (B.25) 
n! 

n>0 

If — 1|| is small enough, then also F^^Fq has as an isolated eigenvalue. We denote the corresponding spectral 
projection by lo{F^^Fo) and we calculate 

Res{F{z)-^) ^ Res(Fo + zFi)-^ ^ (Res(Ff ^i^o + z)"^) F^^ = lo(i^f Vo)Ff \ (B.26) 

The last expression is clearly a rank-one operator. In the case at hand, F^^ = l + o(|e|'^), ase\0, which yields 

We proceed to the proof of Lemma IBIT] 

First, we choose the rate g such that fN<9< ininl^yi, gjv} and we fix the contours F-^ and F^ (see also Figure 
K 

• The contour F-'^ is as described in Lemma [B.3[ with r < \g~ /nI- In particular, for small e, it encircles the point 
/ but no other singular points of A(z). 

• The contour F^ is given by F^ := — e.g + iM. 

By Lemma [B. 2 1 we know that for small e, there are no singularities of A{z) in the region Re z > —eg except for 
the point z = f. Hence, we can deform contours as follows 

Vit,e) = dze^*A(z,e) (B.27) 

= — /dze^*A(z,e) + — /dze^*A(z,e) (B.28) 

27ri J ZTTl J 

rf 

The first term in (IB.28I I yields e*-^ P. The second term of (IB.28I I is split as follows 

j dze^*A(z,e) (B.29) 

= j dze^*(z-iB-eiV)-i (B.30) 

+ j dze^\z-iB-eN)-^{Ai{z,e)-eN)A{z,e) (B.31) 
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Rez 



T 




^9N 




Fi gure 15; The (rotated) complex plane. The black dots indicate the spectrum of iB + eN (which need not be discrete). The upper dot is the 
eigenvalue efpf. In the picture, we have assumed that the spectrum of B consists of 3 semisimple eigenvalues: 0, bi, The gray patches 
contain the possible singularities of the function A{z) above the irregular black line. These singularities lie at o{e) from the spectrum of iB + tN. 
Below the irregular black line, i.e., in the region Kc z < —eg a, we have no control since A^z, e) ceases to be analytic in that region (hence we 
have also not drawn a patch around fe- 1 ). The integration contours F"* , r_> and are drawn as dashed lines. 
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The term l|B.30t equals 

since the contour can be closed in the lower half-plane to enclose the spectrum of iB + eN minus the eigenvalue 
efN, i.e., the set eil^. 

The integrand of lIB.Slb decays as \z\^^ for z oo , since for a bounded operator M 

=0(^), |z|^oo (B.33) 

Using that Ai{z, e) = 0{e), it is now easy to establish that the integral in (|B.3Hl is 0(1), as e \ 0. One extracts e*^''^ 
from the integration | |B.31| I to get the boimd 0(e~^^*). Together with | |B.32t , this proves Lemma IBJ] 

C Appendix: Construction and analysis of the Lindblad generator M. 

The operator A4 was introduced at the beginning of Section ID We provide a more explicit construction and we 
prove Propositions 14. Il and l4. 21 

C.l Construction of M 

First, we note that by using the notions introduced in Section |5^ the operator £{z), defined in Section |4J1 can be 
expressed as 

Ciz) = J dtc-'^' J2 i^*{x2-xi,t)X,,,i,e~'^^^''^'X,,,i, (C.l) 

where ip'^ equals V' or tp, depending on li, I2, according to the rules in l|5.10l l. In words, A^£(z) contains the terms 
of order in the Lie-Schwinger series of Lemma [2.51 
Next, we define some auxiliary objects. 



T := Im WaW* j di V(0, i)e'"* (C.2) 

aesp(ad(y)) "^^^ 

"^^P) - ( f dte^^'i^ix^-xj^t)] X {l,^^Wa)p{l.^^War (C.3) 

s^.i^GZ'' aesp(ad(y)) ^ * ^ 

The operator T = T* e was already referred to in Section ID From the above expression and the 

definition of Wa in ( |2.30t , we check immediately that [Y, T] = 0. Further, we can rewrite | |C.3l l as 

5'(/9) = 27r Y J ds^p{a)V{s,a)pV*{s,a) (C.4) 

aesp(ad(Y)) g<i_i 

with Vi(-) as in l|2!27t and ^i^, and 

y(s,a) := ^ e'"^-n,®M^„ (C.5) 

The expression | |C.4| | is essentially the Kraus decomposition of "J*, see fll/ and hence it shows that is a completely 
positive map. This means in particular that for p > 0, 1 = Tr and hence, by using (IC.3I I, 

=5]Tr^.((*p)(x,2:)) < dtmO,mwr^ IIpIIi (C.6) 
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where the finiteness of the factor between brackets on the RHS is implied by Lemma 15.21 Since any trace class 
operator can be written as a linear combination of four positive trace class operators, it follows that "if is bounded 
on ^i(J^) and by a similar calculation, one can check that ^ is also bounded on 3§2{J^)- 
We are now ready to verify that 

M{p) = -i[eiP) + T,p]+^ip)-^i^*il)p + p^*il)). (C.7) 

Indeed, this is checked most conveniently starting from (|4.41 and employing ||C.1| |. The terms with h ^ h give rise 
to '^{p), while the terms with li = I2 give rise to — i[T,/9] and — + p'^*{l)). Moreover, by the bounded- 
ness and complete positivity of 'f and the representation | |C.7| |, it follows that A4 is of Lindblad type, see e.g. [IJ. 
Starting from | |C.4| |, one can derive the momentum space representation of Ai given in Section l4!2l For example, 
by expressing V{s,a) in momentum representation, one obtains 

^{p){kL,ku)^27T J dsi^{a)WaPiki^+sa,kj, + sa)W: (C.8) 

aesp(ad(y)) gd-l 

which gives rise to the first term of l|4.16|l . 
C.1.1 Proof of Proposition O 

By the integrability in time of the correlation function iIj{x, t), as stated in Lemma [5.2[ the expression dC.lt implies 
immediately that C{z) can be continued continuously to z G M. This proves | |4.5l l. The boundedness of A4 on 
^2i'^) and ^i{J^) follows from the boundedness of 5*, which was explained above. The complete positivity 
of the map 5* and the canonical form jC.T} imply that is a Lindblad generator, see e.g. [IJ. Consequently, 
— iad(y) + X^Ai is also a Lindblad generator and the semigroup At is positivity-preserving and trace-preserving. 
To check ||4.6|I , we note that 

JnMJ-^ -M^-i [J^ad{e{P))J-^ ~ ad(e(P))] (C.9) 
and hence ( 14.61 1 follows immediately from Assumption l2.1l 

C.2 Spectral analysis and proof of Proposition 14.21 

The claims of Proposition l4.2l require a spectral analysis which we present now. We recall the decomposition 

M= dpMp^ dp © Mp^a (C.IO) 

Jjii JTi aesp(ad(y)) 

and we keep in mind that Proposition 142] treats Mp as an operator on the Hilbert space ^ ^ L^(T'^, ^2(-5^))- 
C.2.1 Explicit representation of Mp,o 

By exploiting the nondegeneracy condition in Assumption l2.4[ we can identify A4p.a=Q, for each p with an operator 
on L'^{T^ X spy). This was explained in Section l42l We introduce explicitly gam, loss and kinetic operators operators; 
G, L and Kp, acting on (jd x spy), by 

G^{k,e) ■■= J2 [ dfcV(fc',e';fc,e)(^(fc',e') (C.ll) 
e'espY •^T'^ 

Lip{k,e) - / dk'r{k,e;k',e')(p{k,e) (C.12) 

e'espY J'^" 

Kpip{k,e) := i(e(fc + |)-e(fc-|))^(fc,e), ip e L'' (T"" x spY) (C.13) 
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The kinetic operator Kp models the free flight of the particle between collisions. The operators L and Kp act by 
multiplication in the variables k, e. The expression for Mp^, given in given in Section l4!2l (in particular in (|4.151 ), 
can be rewritten as 

Mp,o ^G + L + Kp (C.14) 

We define the similarity transformation 

A:=ei'^^Ae-i'^^, for A e ■^{L^iT'^ x spY)) (C.15) 

where we have slightly abused the notation by writing Y to denote a multiplication operator on spY, i.e., Yip{k, e) = 
e(p{k, e). Since L and Kp act by multiplication, we have L = L and Kp = Kp. The usefulness of this similarity trans- 
formation resides in the fact that G, and hence also Ma.or are self-adjoint on L^(T'^ x spF). 

C.2.2 Explicit representation of Mp.a^a 

To write an explicit expression for Aip^a^o/ we first define the operators (acting on 

adJT) - l„(ad(r))ad(T)l,(ad(r)), a G sp(ad(r)) (C.16) 

which satisfy ado (T) = and ad(T) — ®aada (T) since [T, F] = 0. Due to the nondegeneracy condition in Assump- 
tion l2.4i both la(ad(F)) and ada(T) are rank-one operators and we identify ada(T) with a number Tq, such that 
ada(T) = Tala(ad(Y')). In fact, as already remarked in Section l42l the operator Mp^a^o itself acts as the rank-one 
operator la(ad(y)) on and hence we identify it with an operator on i^(T'^) (which is also called Mp^a^o 

here): 

Mp,a^ov{k) = -i{Ta + e{k + |) - e{k - - i {j{k, e) + e')) V & L^T") (C.17) 

where j{k, e) are the escape rates introduced in l|4.27|l and e, e' are determined by a = e — e'. To check | |C.17| |, one 
starts from | |4.15| | and one uses 

• The fact that r{k', e'; k, e) vanishes for e' = e, as remarked following (|4.171 . 

• The definition of the matrices Wa in | |2.30|I and the escape rates •) in l|4.27|l . 

• The definition ip{k) = (e, ^ik)e')j^ in l l423ll . 

In particular, the last term on the RHS of | |C.17|| appears because 

/ dk'r,,.,{k,k')(e',{WaW:ak)+ak)WaW:)e) =(j(fc,e)+j-(fc,e'))^(fc) (C.18) 
Hence, Pip.ai^Q acts by multiplication in the variable k. 
C.2.3 Analysis of TWq o 

We already established that Alo,o is a bounded Markov generator on L^(T'^ x spy). This implies that 

Resp£iA^o,o < 0. (C.19) 

The operator A^o,o is not longer a Markov generator, but its spectrum is identical to A^o.O/ since e*^'^^ is bounded 
and invertible. The loss operator L is a multiplication operator and its spectrum is found to be (see l|4.27|l ) 

sp(L) = - { j(e, fc) I e e spy, k e T''} (C.20) 

It is important to note that the escape rates j(e, fc) are bounded away from 0; this is a consequence of Assumption 
I2.4l and more concretely, of the fact that for each e, there is a e' such that l|2.34b holds. Hence, we have 

sp(L) < (C.21) 
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Next, we argue that G is a compact operator on L"^. Indeed, for fixed e, e', the kernel r(fc', e'; k, e) depends only on 
Afc = k — k' and, hence, its Fourier transform acts on ^(2"^) by multiplication with the fimction 

Z'^3x^ I d(Afc)e''^'=-^r(0,e'; A/c,e) (C.22) 

From the explicit expression for r{k\ e'; fc, e), one checks that the function l|C.22l l decays at infinity if the dimension 
d> 1 (recall that d> Ahy Assumption l2.2|l . Hence G is compact. 

Given the compactness of G, Weyl's theorem ensures that the self-adjoint operators A^o.o and L have the same 
essential spectrum. 

By inspection, we check that TMq.o has an eigenvalue 0, corresponding to the eigenvector (p^''(fc,e) = e~'^'^. 
Note that the corresponding right eigenvector of Alo.o is the (unnormalized) Gibbs state (p'^''{k, e) = e^^'^ and the 
corresponding left eigenvector is the constant function, since indeed 

^■^^ = e^^V^^ = e-^^'lT.xspy (C.23) 

Sinceanyeigenvalueof jMo,o oni^ has to be an eigenvalue of 7Mo,o on (Note that i^(T''xspy) c L^(T'^ x spF)), 
the relation llC.19l l implies that there are no eigenvalues with strictly positive real part. 

We now exploit a Perron-Frobenius type of argument to argue that the eigenvalue is simple and that it is 
the only eigenvalue on the real axis. See e.g. Theorem 13.3.6 in |8| for a version of the Perron-Frobenius theorem 
that establishes this in our case, provided that the semigroup e*^" " is irreducible, i.e., that for any nonnegative 
fimctions ip G L^{T'^ x spY), the inclusion 

Supp (e*^" V) C Supp((^), (Supp stands for 'support') (C.24) 

implies that either Supp(iy9) = T'' x spY or (p = 0. This irreducibility criterion is easily checked starting from 
Assumption 12 .41 in particular its rephrasing in terms of a connected graph. Theorem 13.3.6 yields that the eigen- 
value 1 of e*-^" " is simple, which implies that the eigenvalue of A^o,o is simple. To exclude purely imaginary 
eigenvalues ib of A4a,o, we apply this theorem for t such that c'''* = 1. 



C.2.4 Analysis of Mp^o and Mp^ 

We investigate the spectrum of Xip.o as follows. By the same reasoning as in Section |C.2.31 any spectrum with real 
part greater than (the negative number) supspL consists of eigenvalues of finite multiplicity. Assume that A4p has 
an eigenvalue nip with (right) eigenvector) ipp. Then 

Remp{ipp,ipp) = Re{ipp,Mpflipp) (C.25) 
Re{ipp,Kpipp) +Re{ipp,Mo,oV'p) (C.26) 
The first term in llC.261 vanishes because the multiplication operator Kp is purely imaginary. The second term 
can only become positive if = (^^"J, with ip'^'^ the eigenvector of A^o.o corresponding to the eigenvalue 0. This 
means that either the eigenvalue rup has strictly negative real part, or the vector (f^'^ is an eigenvector of Mq.o with 
eigenvalue 0. In the latter case, (p'^'^ must also be an eigenvector of Kp with eigenvalue 0, which can only hold if 
e{k + §) — e(fc — |) = for all k. This is however excluded by the condition (|2.11l l in Assumption l2.1l 

We conclude that for all p G T^X {0}, we have Rc spAip < 0. By compactness of T'* and the lower semicontinuity 
of the spectrum, we deduce hence that 

sup Re spA^p.o — c{Iq) < 0, for any neighborhood Iq of (C.27) 
For a 7^ 0, the operator Mp^a^o is a multiplication operator in k and 

RespA^p.Q < — i inf j(fc, e) < 0, independently of p (C.28) 
as follows by | |C.17| | and the fact that j{k, e) is bounded away from 0. 
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C.2.5 Proof of PropositionH^l 



We summarize the results of Sections lC.2.3l and lC.2.4l For a 7^ 0, the real part of the spectrum of the operators Aip^a 
is strictly negative, uniformly in p, see (|C.28l l. The real part of the spectrum of Aip^ is strictly negative, uniformly 
in p except for a neighborhood of e T'*. 

The operator A4o,q has a simple eigenvalue at with corresponding eigenvector ip'^'^, as defined in Section |C.2.3[ 
The rest of the spectrum of Aio,Q is separated from the eigenvalue by a gap. 

Since A4p — ®a-A^p,a/ and using the uniform bound ||C.28|| , we obtain immediately that the operator A4q has a 
simple eigenvalue at with corresponding eigenvector ^'^'^ 

f 9 := <^«9 © ffi . . . ® 0, (C.29) 

separated from the rest of the spectrum of Mo by a gap. By the analyticity in k, see l|4.61 , and the correspondance 
between k and {p, v), as stated in l|2.60t , we can apply analytic perturbation theory in p to the family of operators 
Mp. We conclude that for p in a neighborhood of 0, the operator Mp has a simple eigenvalue, which we call 
frwip), that is separated by a gap from the rest of the spectrum. We also obtain that the corresponding eigenvector 
is analytic in p and v. 

In this way we have derived all claims of Proposition 1121 except for the symmetry Vpjrwip) = and the strict 
postive-definiteness of the matrix {Vp)'^frw{p)- These two claims will be proven in Section lC.2.6l We note that the 
function frw{p), which we defined above as the simple and isolated eigenvalue of M.p with maximal real part, is 
also a simple and isolated eigenvalue o{ M.pfi with maximal real part. 

C.2.6 Strict positivity of the diffusion constant 

By the remark at the end of Section lC. 2.51 and the fact that spA^p.o — spA^j,,o/ we view frw{p) as the eigenvalue of 
Mpfi that reduces to for p — Q. 

We recall that Mp^o — A^o,o + Kp and we define the operator-valued vector V :— ^pKp\^^Q (note that y is in 
fact a vector of operators). The first order shift of the eigenvalue is given by 

To check that l|C.30|l indeed vanishes, we use that (p'^^ is symmetric under the transformation k ^ —k (in fact, it 
is independent of k) while V is anti-symmetric under k ^ —k (this follows from the symmetry e{k) = e{—k) in 
Assumption 12.11 1 . 

The second order shift is then given by 

where the first term on the RHS of llC.311 is well-defined since V(p'^'^ is orthogonal to the 0-spectral subspace of 
A^o.O/ by I IC.3OI 1. The second term vanishes because {\7p)^Kp = 0, as can again be checked explicitly. 
Let V e R'' and V^, -.^ v ■ V (recall that y is a vector). Then, by l lOTl l, 

^ ■ = -(^^^^^^'''^^-^o.oK^^') (C.32) 

Upon using the spectral theorem and the gap for the self-adjoint operator A^o.O/ we see that the RHS of the last 
expression is positive and it can only vanish if 

which is however excluded by Assumption 12.11 The strict positive-definiteness of the diffusion constant Drw is 
hence proven. 



-/3e 



dk\v ■Ve{k)\'' 



(C.33) 
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D Appendix: Combinatorics 



In this appendix, we show how to integrate over irreducible equivalence classes of diagrams. In other words, 
we assume that the x, /-coordinates have already been summed over (or a supremum over them has been taken) 
and we carry out the remaining integration over the time-coordinates t and the diagram size \a\. We first define a 
function of diagrams, K{a), that depends only on the equivalence class [a]. Let khe a positive fimction on IR+ and 
put 

kl 

K{a) :^Ylkiv,~u,) (D.l) 

i=l 

where {ui,Vi) are the pairs of times in the diagram a. In the applications, the function k will be (a multiple of) 
sup^ \ip{x, t)\, sometimes restricted to t < t or i > r. 

Lemma D.l. Let a > and assume that ||te°*fc|| i = J^^ dt te"*fc(i) < 1, then 

H+ nTS[o,tj(mir) 

If in addition, ||te°*/c||i < 1 with a a+ \\k\\i, then 

dte- / dMi^(a) < 2|le^'fc|li^-^^ (D.3) 

nrS[o_t] (ir) 

die- / d[a]l|.|>,if(a) < 2||e'^^fc||i /^g!fc||^ (D-4) 

E+ nrS[o,tj(ir) 

Proof. First, we note that for each irreducible diagram a e E[o^t](ir), we can find a subdiagram a' c a such that 
cr' is minimally irreducible in [0, t], i.e., a' G Sjo.f] (mir). Note that the choice of subdiagram a' is not necessarily 
unique. Conversely, given a minimally irreducible diagram cr' G ^lo,t] (mir), we can add any diagram a" E ^lo,t] to 
cr', thereby creating a new irreducible diagram <j := <j' U a" E S[o,t] (ir)- By these considerations, we easily deduce 

J d[a] l|,|>2if(r7) < I J d[a']l|,,|>2if(a')j M + J d[<j"]K{a")] (D.5) 

+ [ / d[a']l|.,|=ii^(a')V 1 d{a"\K{a")\ (D.6) 

Thel+in||D3J covers the case in which the diagram a was itself minimally irreducible, and hence no diagrams a" 
are added to a' . In llD.61 , one always has to add at least one pair to cr', since |cr| > 2 but |cr'| — 1. In fact, the equival- 
ence classes in the inequality could be dropped, i.e., one can omit the projections IIt and replace d[cr], d[cr'], d[cr"] 
by dcr, dcr', dcr", respectively. 

We recall that if a diagram a with |cr| = 1 is irreducible (or minimally irreducible) in the interval /, then its 
time-coordinates are fixed to be the boundaries of /; i.e., there is only one equivalence class of such diagrams. 
Hence 

y d[a]l|„|=ii^(a) - y d[a]l|,|=ii^(a) = fc(t) (D.7) 

nTS[o,tj(ir) nrS[o,t](mir) 

The unconstrained integral over all (equivalence classes of) diagrams, that appears in l ID.SI l and 1ID.6I 1, can be 
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performed as follows 



J dad[a]K{<7) = ^ J dui . . . dun J dvi . . . dw„ ( Y[ Hvi ~ Ui) j 

nTS[o_t] 0<ni<...<«„<t Vi>Ui J 

E / dui...d«„(iifciiir = ES(ii^iii)" = ^*"'"^-i 



< 



Next, we perform the integral over (equivalence classes of) minimally irreducible diagrams. For a G ^[o.t] (niir) 
with \a\ = n > 1, the relative order of the times Ui, vt is fixed as follows: 

— Ui < U2 < Vi < U'i < V2 < Ua < ■ ■ ■ < < Un ^ I'n-l ^ Vn = t (D.9) 

We have hence 



/ die"* / daX(a)l|^|^^ - / dvik{vi - ui)^''^''^-'''^ / d?/2 / dz^s ■ ■ • 

air) 

poo PVn — 2 rOO 

du„_2 / dv„_2 . . . / du„_i / dD„_ie°'^^"-i"''"-^^A:(i;„_i - m„_i) 
du„ / di;„e"('^"-'^'-i)fcK-w„). (D.IO) 



i;„_2 Ju„_i 



First we extend the domain of integration of u„ from [vn^2,Vn-i] to (— oo, and we estimate the integrals over 
the variables w„ and w„ by 



dun / d^;„e'^(''"-''"-^)/c(i>„ - Un) < ||te'^*/c||i (D.ll) 



Next, we perform the integration over Un-i,Vn-i in the same way, we continue the procedure until only the vari- 
able vi is left (note that wi = is fixed). The integral over vi gives Ije"*/;;!! i. This yields the bound 

LHS of jPlOl l <\\e''*k\\ix\\te''*k\\'l-^ (D.12) 

We are ready to evaluate the Laplace transform of | |D.5| I-| |D.6| I. Using dP.SI I, we bound 

^1+ J d[cr]X((7)^ < e*ll'=ll\ ^ J d[cr]X(cr)^ < e*ll'''lli - 1 < t||/c||ie*ll''^lli (D.13) 

Combining this with ||D.7| | and ||D.12| | , and summing over n>2, we obtain 

Jdte'^' J dMl|„|>2if(a)<||e"*fc||i^^H21|^ + ||A;||i||te'^*fc||i (D.14) 

nTS[o,t](ir) 

where the two terms on the RHS correspond to $D.5\ and llD.6b , respectively. This ends the proof of (|D.4|l . The 
bound in llD.31 follows by adding ||c°*A;||i, which is the contribution of \a\ = 1 (see l|D.7|l ). to l|D.4) l. The bound l|D.2)l 
is proven by summing l|D.12|l over n > 1. 

□ 
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